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Summary 


Objectives:  The  objective  of  this  research  was  to  develop  input  configurations  for  the  Draper 
Semianalytic  Satellite  Theory  (DSST)  which  provide  a  balance  between  computational  speed  and 
ephemeris  accuracy.  This  research  follows  VTA’s  efforts  in  1994  to  port  the  software  system 
containing  DSST  (Draper  Research  and  Development  Goddard  Trajectory  Determination  System, 
R&D  GTDS)  from  a  UNIX-based  platform  to  a  Windows/DOS  environment  (486  and  Pentium). 
These  input  decks,  which  are  required  to  properly  run  the  software,  have  been  developed  for 
distribution  along  with  a  complete  mathematical  description  of  theory  (the  description  was 
developed  by  the  Naval  Postgraduate  School  with  assistance  from  Draper  Laboratory).  Optimal 
input  decks  were  developed  for  low,  medium,  and  high  altitude  circular,  geosynchronous, 
Molniya,  and  near  earth  eccentric  orbits. 

Approach:  In  order  to  determine  the  proper  configurations  for  the  theory,  various  batch, 
differential  correction  fits  to  simulated  observational  data  were  performed.  The  results  of  these 
various  fits  were  then  compared  in  terms  of  accuracy  and  computational  efficiency  to  operational 
orbit  determination  capabilities.  Selection  of  appropriate  input  deck  configurations  was  based  on 
attaining  speeds  comparable  to  operational  algorithms  while  maintaining  at  least  a  5-10  fold 
increase  in  ephemeris  accuracy. 

Alternative  Technologies:  Various  alternative  technologies  for  orbit  determination  (OD)  are 
available;  these  OD  technologies  can  be  compared  based  on  the  orbit  propagation  technique  used. 
Classically,  orbit  propagation  has  been  categorized  as  either  numeric  (very  high  accuracy  with 
slow  run  times)  or  analytic  (very  fast  with  severely  limited  accuracy).  DSST  represents  a  hybrid 
approach  which  combines  the  advantageous  aspects  of  both  numeric  and  analytic  theories 
(accuracies  approaching  that  of  numeric  methods  with  speeds  comparable  to  analytic  techniques). 
DSST,  developed  by  Dr.  Paul  Cefola  and  his  colleagues  at  the  Charles  Stark  Draper  Laboratory,  is 
the  state  of  the  art  in  American  semianalytic  techniques. 

Technical  Challenges:  VTA  has  had  to  combat  (1)  a  limited  knowledge  base  concerning  the 
cutting  edge  astrodynamic  techniques  of  DSST  and  (2)  a  lack  of  knowledge  of  how  to  properly 
run  the  software  with  potential  users. 

Payoffs'  DSST  will  be  used  by  various  research  and  operational  based  institutions.  Currently, 
DSST  is  slated  for  operational  use  with  the  RADARSAT  mission  (work  accomplished  by  Draper 
Laboratory),  the  Air  Force  Maui  Optical  Station  (AMOS),  and  the  STARFIRE  facility.  VTA  has 
been  using  R&D  GTDS  to  support  operational  orbit  determination  experiments  for  the  past  eight 
months.  Specifically,  dark  pass  satellite  illumination  has  been  accomplished.  Kaman  Sciences 
has  been  using  the  input  decks  for  their  Space  Command  supported  Ephemeris  Theory  Accuracy 
Study.  The  Naval  Postgraduate  School  has  been  using  results  generated  with  the  input  decks  to 
increase  Navy  Space  Command’s  understanding  of  semianalytic  techniques. 

Results:  For  the  first  time,  input  decks  to  balance  computational  speed  and  ephemeris  accuracy 
have  been  documented.  This  gives  the  user  with  limited  knowledge  concerning  the  intricacies  of 
DSST  a  better  chance  at  properly  running  the  software. 
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1.0  Introduction 


Various  studies  in  the  past  twenty  years  have  demonstrated  the  accuracy,  computational  efficiency,  and  flexibility 
of  DSST  for  a  variety  of  orbit  types1'24.  These  encouraging  results  have  led  to  a  proliferation  of  the  theory  to  various 
academic  and  research  based  institutions.  In  most  cases,  DSST  has  been  distributed  as  part  of  the  larger  scale  R&D 
GTDS  orbit  determination  system,  which  also  includes  the  Cowell32,  SGP25,  SGP4/  SDP425,26,27,  HANDE28, 
SALT29,30,31,  Brouwer  Lyddane32,  and  Vmti32  orbit  propagation  theories  (the  addition  of  PPT2  to  R&D  GTDS  is 
nearing  completion  as  of  the  date  of  publication  of  this  document).  In  addition  to  these  orbit  propagation  theories, 
R&D  GTDS  includes  estimation  (batch  &  filter),  early  orbit  determination,  data  simulation,  and  error  analysis 
programs32.  These  programs  give  R&D  GTDS  the  flexibility  to  compare  the  performance  of  various  orbit 
propagation  theories. 

Recent  efforts  have  unified  DSST’s  theoretical  documentation  from  a  multitude  of  AJAA/AAS  papers,  journal 
articles,  contract  reports,  laboratory  memorandums,  and  working  notes  to  a  single  document1.  The  knowledge  of  how 
to  properly  run  the  software,  however,  has  not  been  as  easily  grasped  by  the  community.  The  addition  of  the 
Semianalytic  Input  Processor  to  the  software  has  provided  a  tool  which  simplifies  the  configuration  of  DSST’s 
SPG33.  This  document  sheds  light  on  how  to  best  use  the  tool. 

Specifically,  this  paper  contains  a  set  of  DSST  input  decks  which  balance  accuracy  and  computational  efficiency 
for  a  variety  of  orbit  types.  These  configurations  span  the  orbit  classes  established  by  Kaman  Sciences  in  the 
Ephemeris  Theory  Accuracy  Study  (ETAS)5.  These  configurations  can  be  inserted  into  the  code  for  automatic  set-up 
of  the  theory,  or  used  in  a  standalone  fashion  by  the  traditional  R&D  GTDS  user.  The  information  in  this  paper,  when 
coupled  with  the  description  of  the  theory1,  provides  current  and  potential  users  with  a  concise  package  which 
supports  R&D  GTDS. 

One  final  note  before  each  specific  orbit  class  is  described  in  detail  is  that  two  separate  machines  were  used  for 
the  generation  of  results  herein  (both  a  486  and  a  Pentium  586  machine  were  used).  Therefore,  from  case  to  case, 
timing  statistics  may  not  seem  consistent;  however,  within  one  particular  case,  the  results  are  consistent.  In  addition, 
most  cases  used  a  spacecraft  cross-sectional  area  of  1  m2  and  a  mass  of  100  kg.  Even  though  these  values  may  not 
exactly  match  the  parameters  for  the  actual,  “real-world”  spacecraft,  they  were  chosen  to  represent  a  standard  area  to 
mass  ratio.  Furthermore,  appendices  are  attached  which  provide  the  information  necessaiy  to  replicate  test  cases  with 
the  optimal  input  decks. 

2.0  Low  Altitude  Circular  Case  (LACC) 

An  orbit  very  close  to  decay  was  chosen  to  fulfill  the  low  altitude  circular  test  case.  In  the  ETAS  study,  the  LACC 
is  characterized  by  a  SMAH  range  of  0  -  575  km  and  an  eccentricity  range  from  0.0  to  0.05.  For  this  test  case, 
multiple  runs  were  made  in  which  varying  configurations  of  DSST  were  fit  to  a  Cowell  truth  model  over  a  90  minute 
arc.  A  90  minute  DSST  EG  with  the  DC  solve-for  parameters  followed  the  fit  span.  Accuracy  metrics  were  derived 
from  comparisons  between  the  DSST  and  Cowell  trajectories  over  both  the  fit  (DC  RMS)  and  predict  spans  (from  the 
R&D  GTDS  Ephemeris  Comparison  Program).  ORB  1  output  was  produced  every  450  seconds. 

The  osculating  elements  required  for  the  Cowell  truth  EG  were  built  from  a  one  day  DSST  run  with  the 
following  inputs: 

Table  1.  DSST  Specifics  to  Build  Osculating  Elements,  LACC 


Epoch  Date  YYMMDD  HHMMSS.S 

820223  000000.0 

End  Date  YYMMDD  HHMMSS.S 

820224  000000.0 

Semimajor  Axis  (mean) 

6635.0814  km 

Eccentricity  (mean) 

0.010201164 
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Table  1.  DSST  Specifics  to  Build  Osculating  Elements,  LACC  (Continued) 


Inclination  (mean) 

64.9567  deg 

Longitude  of  Ascending  Node  (mean) 

228.6393  deg 

Argument  of  Perigee  (mean) 

271.2229  deg 

Mean  Anomaly  (mean) 

88.164558  deg 

Input  Frame 

True  of  Date 

Integration  and  Output  Frame 

Earth  Equator  and 

Equinox  of  1950 
(Mean  of  1950) 

Integrator 

4th  Order  Runge-Kutta 

DSST  Step  Size 

43200.0  sec 

Gravity  Field 

8x8  GEM10B 

Lunar  /  Solar  Point  Masses 

Yes 

Drag  Model 

Jacchia-Roberts 

CD  =  2.0 

Solar  Radiation  Pressure 

No 

Global  Short  Periodic  Select 

Low  Altitude,  Improved 
Accuracy  Option  (SPSHPER  =  3) 

Average  Spacecraft  Cross-Sectional  Area 

1.0  m2 

Spacecraft  Mass 

100  kg 

The  osculating  elements  which  resulted  from  this  one  day  DSST  run  were  then  used  as  inputs  for  a  six  day 
Cowell  EG.  The  specifics  for  this  Cowell  EG  are  given  in  Table  2: 

Table  2.  Cowell  Truth  EG  Inputs,  LACC 


Epoch  Date  YYMMDD  HHMMSS.S 

820224  000000.0 

End  Date  YYMMDD  HHMMSS.S 

820302  000000.0 

Input,  Integration,  and  Output  Frame 

Earth  Equator  and 

Equinox  of  1950 
(Mean  of  1950) 

Integrator 

12th  Order  Summed  Cowell/ 
Adams  Predict-Partially  Correct 

Step  Size 

60.0  sec 

Gravity  Field 

21x21  GEM10B 

Lunar  /  Solar  Point  Masses 

Yes 

Table  2.  Cowell  Truth  EG  Inputs,  LACC  (Continued) 


Drag  Model 

Jacchia-Roberts 

Cd  =  2.0 

Solar  Radiation  Pressure 

No 

Average  Spacecraft  Cross-Sectional  Area 

1.0  m2 

Spacecraft  Mass 

100  kg 

It  can  be  noted  in  this  table  a  six  day  Cowell  EG  was  created,  even  thought  the  fit  and  predict  spans  totalled  180 
minutes  (in  other  words,  more  ORB1  data  was  generated  than  was  required).  This  test  procedure  was  established  to 
provide  enough  data  if  subsequent  testing  with  other  fit  and  predict  spans  was  desired.  The  output  of  this  EG 
indicated  the  satellite  would  re-enter  the  atmosphere  three  days  after  the  Cowell  epoch  (this  particular  test  case 
reflects  an  orbit  extremely  close  to  decay).  Due  to  the  short  remaining  orbital  lifetime  for  this  satellite,  it  was 
determined  that  fit  and  predict  spans  on  the  order  of  the  satellite’s  period  should  be  used.  This  test  protocol  was 
deemed  analogous  to  operational  procedures  for  a  satellite  this  close  to  re-entering  the  atmosphere. 

The  same  initial  guess  of  the  state  vector  and  observation  noise  standard  deviation  were  used  for  all  of  the  DSST 
DC  runs.  These  values  are  highlighted  in  Table  3: 

Table  3.  DSST  DC  Specifics,  LACC 


Epoch  Date  YYMMDD  HHMMSS.S 

820224  000000.0 

Fit  Span  Begin  YYMMDD  HHMMSS.S 

820224  000000.0 

Fit  Span  End  YYMMDD  HHMMSS.S 

820224  013000.0 

Semimajor  Axis  (mean) 

6628.457  km 

Eccentricity  (mean) 

0.0089 

Inclination  (mean) 

64.84  deg 

Longitude  of  Ascending  Node  (mean) 

224.51  deg 

Argument  of  Perigee  (mean) 

271.89  deg 

Mean  Anomaly  (mean) 

115.16  deg 

DSST  Input,  Integration,  and  Output  Frame 

Earth  Equator  and 

Equinox  of  1950 
(Mean  of  1950) 

DSST  Step  Size 

43200.0  sec 

Geopotential  Model 

GEM10B 

Position  Standard  Deviation 

100  meters 

Velocity  Standard  Deviation 

10  cm/sec 

In  the  DSST  DC  runs,  vaiying  configurations  corresponding  to  different  perturbation  modeling  were  used  to 
determine  the  optimal  balance  of  accuracy  and  computational  speed.  Results  of  the  DC’s  and  subsequent  EG’s  for  the 
LACC  are  given  in  Table  4;  however,  it  is  first  desirable  to  describe  some  of  the  notation  used  in  this  table  (and  other 
tables  in  this  paper): 
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•  GRAV  refers  to  the  geopotential  modeling  in  the  AOG. 

•  RESONPRD  is  a  GTDS  keyword  for  the  input  processor  to  override  the  default  boundary  for  the  modeling  of 
tesseral  resonance  in  the  AOG  (if  the  period  of  the  resonance  is  less  than  10  days  in  the  default  configuration, 
the  resonance  is  modelled  in  the  SPG  as  a  tesseral  linear  combination  high  frequency  term);  by  specifying  a 
value  with  this  keyword,  a  new  boundary  can  be  set  to  force  the  modeling  of  resonance  into  the  AOG  (i.e.,  for 
shallow  resonance). 

•  RESNM  is  a  GTDS  keyword  for  the  input  processor  to  include  resonant  harmonics  in  the  AOG.  Typically,  this 
card  is  used  to  augment  a  small  “base”  gravity  field  (i.e.,  4x4)  specified  with  MAXDEGEQ  and  MAXORDEQ 
cards  with  resonant  harmonics  beyond  the  “base”  field  (for  example,  the  16th  order).  It  should  be  noted  the 
RESNM  card  eliminates  the  need  for  a  RESONPRD  card;  all  harmonics  specified  with  the  RESNM  card  are 
forced  into  the  AOG  (even  shallow  resonance).  For  this  particular  case,  an  8x8  base  field  was  established  and 
the  RESNM  card  was  used  to  include  resonant  effects  at  the  16th  and  17th  orders  (an  AUTOFORC  card  must 
be  used  with  RESNM). 

•  TESS  refers  to  the  tesseral  linear  combination  high  frequency  short  periodic  terms. 

•  SRP  is  an  abbreviation  for  solar  radiation  pressure. 

•  J-R  is  an  abbreviation  for  the  Jacchia-Roberts  atmospheric  density  model. 

•  H-P  is  an  abbreviation  for  the  Harris-Priester  atmospheric  density  model. 

•  L-S  is  an  abbreviation  for  lunar  /  solar  third  body  point  mass  effects. 

•  “trunc”  refers  to  truncated  modeling;  modeling  descriptions  corresponding  to  the  R&D  GTDS  Semianalytic 
Satellite  Theory  input  processor  are  specified  where  needed.  Additional  details  can  be  found  in  the  R&D 
GTDS  Semianalytic  Theory  Input  Processor  keyword  description  booklet33. 

•  Total  Position  RMS  is  derived  over  the  90  minute  predict  span  only  (from  the  R&D  GTDS  Ephemeris 
Comparison  Program). 

•  R&D  GTDS  default  configuration  settings  are  used  unless  specifically  noted. 

•  “Iszak”  refers  to  Iszak’s  J2  height  correction  (a  second  order  atmospheric  drag  effect  included  in  AOG 
parti als).  All  cases  which  attempt  to  solve  for  CD  include  first  order  drag  effects  in  the  AOG  partials  (as  well  as 
Iszak’s  corrections  if  indicated;  if  Cq  is  not  solved  for,  then  all  drag  partials  in  the  AOG  are  shut  off).  All  cases 
include  J2  partials  analytically  in  the  AOG.  No  other  AOG  or  SPG  partials  were  included  in  any  of  the  runs. 

•  Most  cases  have  an  SSTESTFL  card  of  1  3;  this  setting  includes  all  desired  AOG  partials  in  the  A  matrix  (J2 
partials  are  computed  analytically;  all  other  partial  derivatives  are  computed  via  finite  differences).  In  addition, 
this  setting  includes  desired  AOG  partial  derivatives  for  non-conservative  perturbations  in  the  D  matrix  via 
finite  differences  (these  partials  were  also  included  in  the  A  matrix).  Cases  noted  by  **  have  an  SSTESTFL 
card  of  1  2;  this  setting  puts  only  J2  AOG  partials  into  the  A  matrix  (computed  analytically).  All  other  desired 
AOG  partials  for  non-conservative  perturbations  are  included  in  the  D  matrix  via  finite  differences  (with  this 
setting,  third  body  and  J22  AOG  partial  derivatives,  if  desired,  would  not  be  included  anywhere;  if  they  are 
desired,  the  SSTESTFL  setting  of  1  3  should  be  used). 
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Table  4.  DSST  DC  and  Subsequent  EG  Results,  LACC 


GRAV 

(AOG) 

Zonal 

SP 

M- 

daily 

SP 

TESS 

SP 

h2 

SP 

V 

m-daily 

SP 

Drag 

SRP 

3rd 

Body 

Run 

Time 

DC  RMS 

Total 

Position 

RMS 

(Predict 

Span) 

21x21 

RESON 

PRD  = 

1  day 

Yes 

Yes 

Yes 

Yes 

Yes 

J-R 

AOG 

SPG 

Solve 

CD 

No 

L-S 

AOG 

SPG 

I 

DCDIVERG 

-ES 

21x21 

T  2 

J2 

RESON 

PRD  = 

1  day 

Yes 

Yes 

Yes 

Yes 

Yes 

J-R 

AOG 

SPG 

Iszak 

Solve 

CD 

No 

L-S 

AOG 

SPG 

54.62 

sec 

3.01 

meters 

(6  its) 

33.960 

meters 

21x21 

h2 

Yes 

Yes 

Yes 

Yes 

Yes 

J-R 

AOG 

SPG 

Iszak 

Solve 

cD 

No 

L-S 

AOG 

SPG 

54.33 

sec 

2.66 

meters 

(6  its) 

34.034 

meters 

21x21 

h2 

Yes 

Yes 

Yes 

Yes 

Yes 

J-R 

AOG 

SPG 

Iszak 

Solve 

CD 

No 

No 

53.83 

sec 

2.81 

meters 

(6  its) 

34.168 

meters 

Yes 

Yes 

Yes 

Yes 

Yes 

J-R 

AOG 

Iszak 

Solve 

CD 

No 

No 

45.36 

sec 

77.12 

meters 

(4  its) 

163.31 

meters 

a 

Yes 

Yes 

Yes 

Yes 

Yes 

No 

No 

No 

46.26 

sec 

207.81 

meters 

(5  its) 

5238.7 

meters 
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Table  4.  DSST  DC  and  Subsequent  EG  Results,  LACC  (Continued) 


GRAV 

(AOG) 

Zonal 

SP 

M- 

daily 

SP 

TESS 

SP 

h2 

SP 

V 

m -daily 

SP 

Drag 

SRP 

3rd 

Body 

Run 

Time 

DC  RMS 

Total 

Position 

RMS 

(Predict 

Span) 

m 

Yes 

Yes 

No 

Yes 

Yes 

J-R 

AOG 

SPG 

Iszak 

Solve 

CD 

No 

No 

31.36 

sec 

50.59 

meters 

(5  its) 

216.51 

meters 

21x21 

T2 

J2 

Yes 

No 

No 

Yes 

No 

J-R 

AOG 

SPG 

Iszak 

Solve 

CD 

No 

No 

17.19 

sec 

32.71 

meters 

(5  its) 

468.24 

meters 

Yes 

No 

Yes 

Yes 

No 

J-R 

AOG 

SPG 

Iszak 

Solve 

CD 

No 

No 

35.53 

sec 

50.16 

meters 

(5  its) 

309.45 

meters 

21x21 

T  2 

J2 

Yes 

Yes 

Yes 

Yes 

No 

J-R 

AOG 

SPG 

Iszak 

Solve 

CD 

No 

No 

39.16 

sec 

3.62 

meters 

(6  its) 

31.292 

meters 

12x12 

h2 

Yes 

Yes 

Yes 

Yes 

No 

J-R 

AOG 

SPG 

Iszak 

Solve 

CD 

No 

No 

21.64 

sec 

8.39 

meters 

(6  its) 

87.282 

meters 

Table  4.  DSST  DC  and  Subsequent  EG  Results,  LACC  (Continued) 


GRAV 

(AOG) 

Zonal 

SP 

M- 

daily 

SP 

TESS 

SP 

j22 

SP 

J2/ 

m-daily 

SP 

Drag 

SRP 

3rd 

Body 

Run 

Time 

DC  RMS 

Total 

Position 

RMS 

(Predict 

Span) 

12x12 

T  ^ 

J2 

Yes 

Yes 

trunc 

8842 

-1212 

Yes 

No 

J-R 

AOG 

SPG 

Iszak 

Solve 

CD 

No 

No 

19.28 

sec 

10.17 

meters 

(6  its) 

58.072 

meters 

12x12 

h2 

Yes 

Yes 

trunc 

4442 

-8  8 

Yes 

No 

J-R 

AOG 

SPG 

Iszak 

Solve 

CD 

No 

No 

16.81 

sec 

22.16 

meters 

(5  its) 

198.45 

meters 

8x8 

h2 

Yes 

Yes 

Yes 

Yes 

No 

J-R 

AOG 

SPG 

Iszak 

Solve 

CD 

No 

No 

17.56 

sec 

18.30 

meters 

(5  its) 

58.201 

meters 

8x8 

h2 

Yes 

Yes 

trunc 

6642 

-1010 

Yes 

No 

J-R 

AOG 

SPG 

Iszak 

Solve 

CD 

No 

No 

16.80 

sec 

17.40 

meters 

(5  its) 

65.710 

meters 

8x8 

J22 
*  * 

(opt. 

deck) 

Yes 

Yes 

trunc 

6642 

-10  10 

Yes 

No 

J-R 

AOG 

SPG 

Iszak 

Solve 

CD 

No 

No 

16.31 

sec 

17.40 

meters 

(6  its) 

65.720 

meters 
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Table  4.  DSST  DC  and  Subsequent  EG  Results,  LACC  (Continued) 


GRAV 

(AOG) 

Zonal 

SP 

M- 

daily 

SP 

TESS 

SP 

h2 

SP 

V 

m -daily 

SP 

Drag 

SRP 

3rd 

Body 

Run 

Time 

DC  RMS 

Total 

Position 

RMS 

(Predict 

Span) 

8x8 

I  2 

J2 

| 

Yes 

trunc 

6642 

-1010 

Yes 

No 

J-R 

AOG 

SPG 

Iszak 

Solve 

CD 

No 

No 

16.91 

sec 

18.25 

meters 

(5  its) 

60.651 

meters 

8x8 

h2 

1 

trunc 

442 

trunc 

6642 

-10  10 

Yes 

No 

J-R 

AOG 

SPG 

Iszak 

Solve 

CD 

No 

No 

16.48 

sec 

23.15 

meters 

(5  its) 

291.89 

meters 

8x8 

h2 

Yes 

Yes 

trunc 

4442 

-88 

Yes 

No 

J-R 

AOG 

SPG 

Iszak 

Solve 

CD 

No 

No 

16.42 

sec 

30.89 

meters 

(5  its) 

123.48 

meters 

8x8 

h2 

Yes 

Yes 

trunc 

2242 

-6  6 

Yes 

No 

J-R 

AOG 

SPG 

Iszak 

Solve 

CD 

No 

No 

15.87 

sec 

34.75 

meters 

(5  its) 

211.94 

meters 

8x8 

h2 

Yes 

Yes 

trunc 

2242 

-6  6 

Yes 

No 

J-R 

AOG 

SPG 

Solve 

cD 

No 

No 

15.98 

sec 

34.76 

meters 

(5  its) 

212.44 

meters 

8x8 

h2 

16th 

order 

RESNM 

Yes 

Yes 

trunc 

2242 

-6  6 

Yes 

No 

J-R 

AOG 

SPG 

Iszak 

Solve 

CD 

No 

No 

15.98 

sec 

35.28 

meters 

(5  its) 

209.94 

meters 
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Table  4.  DSST  DC  and  Subsequent  EG  Results,  LACC  (Continued) 


GRAV 

(AOG) 

Zonal 

SP 

M- 

daily 

SP 

TESS 

SP 

h2 

SP 

J2/ 

m-daily 

SP 

Drag 

SRP 

3rd 

Body 

Run 

Time 

DC  RMS 

Total 

Position 

RMS 

(Predict 

Span) 

1 

Yes 

Yes 

trunc 

2242 

-6  6 

Yes 

No 

J-R 

AOG 

SPG 

Iszak 

Solve 

CD 

No 

No 

16.20 

sec 

34.90 

meters 

(5  its) 

214.43 

meters 

6x6 

T  ^ 

J2 

Yes 

Yes 

Yes 

Yes 

No 

J-R 

AOG 

SPG 

Iszak 

Solve 

CD 

No 

No 

16.64 

sec 

15.22 

meters 

(5  its) 

143.55 

meters 

6x6 

J22 

Yes 

Yes 

I 

trunc 

4442 

-8  8 

Yes 

No 

J-R 

AOG 

SPG 

Iszak 

Solve 

CD 

No 

No 

16.14 

sec 

27.83 

meters 

(5  its) 

86.261 

meters 

6x6 

J22 

Yes 

Yes 

trunc 

2242 

-6  6 

Yes 

No 

J-R 

AOG 

SPG 

Iszak 

Solve 

CD 

No 

No 

15.60 

sec 

36.55 

meters 

(5  its) 

364.52 

meters 

1 

Yes 

Yes 

Yes 

Yes 

No 

J-R 

AOG 

SPG 

Iszak 

Solve 

CD 

No 

No 

16.09 

sec 

22.46 

meters 

(5  its) 

332.57 

meters 

10 


Table  4.  DSST  DC  and  Subsequent  EG  Results,  L ACC  (Continued) 


GRAV 

(AOG) 

Zonal 

SP 

M- 

daily 

SP 

TESS 

SP 

h* 

SP 

hi 

m-daily 

SP 

Drag 

SRP 

3rd 

Body 

Run 

Time 

DC  RMS 

Total 

Position 

RMS 

(Predict 

Span) 

Spi 

ice  Command’s  SGP  (solve  for  drag  term  ndot/2) 

13.29 

sec 

1182.50 

meters 

(4  its) 

12605.0 

meters 

Space  Command’s  SGP4  (solve  for  drag  term  ndot/2) 

13.35 

sec 

138.64 

meters 

(5  its) 

505.19 

meters 

The  results  in  Table  4  enforce  some  important  concepts  for  low  earth  orbits.  First,  large  effects  stem  from 
atmospheric  drag,  tesseral  m-dailies,  and  tesseral  linear  combination  short  periodics.  In  fact,  attempting  to  even 
slightly  truncate  the  drag  or  m-daily  modeling  greatly  impacts  the  accuracy.  On  the  contrary,  these  results  indicate  the 
tesseral  linear  combination  short  periodic  modeling  can  be  truncated  to  vastly  improve  the  computation  time  with  an 
acceptable  degradation  of  accuracy  (for  most  applications). 

Specifically,  the  following  errors  can  be  noted  for  a  one  revolution  predict: 

•  removing  drag  entirely  results  in  -  5  km  error 

•  drag  short  periodics- 130  m  error 

•  tesseral  m-dailies  -  300  m  error 

•  removing  all  tesseral  linear  combination  short  periodic  terms  —  200  m  error 

These  results  also  indicate  the  J2  /  m-daily  model  serves  to  provide  accuracy  closure  with  Cowell,  but  at  an 
expense  of  computation  time.  For  reasons  of  simplicity  and  computational  time  savings,  this  model  can  be  neglected 
in  an  optimized  DSST  configuration. 

Third  body  effects  were  determined  to  not  significantly  impact  these  orbits.  Specifically,  the  inclusion  of  these 
effects  does  not  dramatically  improve  accuracy  or  greatly  increase  computer  run  time.  For  reasons  of  simplicity,  their 
effects  can  also  be  neglected  in  an  optimized  DSST  configuration  for  this  orbit  type. 

Several  comments  can  be  made  about  the  short  nature  of  these  runs  (90  minute  fit  and  predict  spans).  Short  fit 
spans  can  lead  to  “observability”  problems  when  the  DC  attempts  to  solve  for  drag  terms  (or  any  other  desired  solve 
for  parameter,  i.e.  station  biases).  These  “observability”  problems  result  because  the  DC  can  not  sample  or  observe  a 
large  arc  of  the  satellite’s  trajectory,  which  makes  it  difficult  to  discern  or  separate  the  various  contributions  to  the 
satellite’s  motion.  For  analytic  theories,  these  observability  problems  are  compounded  by  their  limited  physical  force 
models.  To  deal  with  observability  problems  or  a  lack  of  physical  force  models,  longer  fit  spans  are  usually 
employed;  however,  longer  fit  spans  may  not  be  practical  or  even  possible  for  decay  orbits.  For  these  reasons,  only 
complete  perturbation  theories  should  be  used  to  process  decay  orbits.  It  should  also  be  mentioned  the  short  nature  of 
these  runs  negated  the  benefit  of  special  shallow  resonance  processing  for  the  truncated  geopotential  cases  (i.e., 
RESONPRD  and  RESNM). 

As  the  Approach  section  of  the  Summary  on  page  1  states,  the  selection  of  appropriate  input  deck  configurations 
was  based  on  attaining  speeds  comparable  to  operational  algorithms  while  maintaining  at  least  a  5-10  fold  increase  in 
ephemeris  accuracy.  The  results  for  this  particular  test  case  and  test  methodology  indicate  that  several  of  the  DSST 
configurations  very  nearly  meet  this  criteria  (the  timing  statistics  fall  short  of  the  desired  criteria  by  20-30  percent). 
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Testing  was  done  to  try  to  simulate  real-world  uncertainties  with  atmospheric  drag.  In  these  tests,  a  Jacchia- 
Roberts  drag  model  was  used  to  generate  a  Cowell  truth  EG.  Then,  a  Hams-Priester  DSST  configuration  was  fit  over 
a  90  minute  arc  to  the  Jacchia-Roberts  Cowell  truth  trajectory  (in  these  cases,  either  the  default  Harris-Priester  table 
with  F#=l  50  or  the  “hottest”  table  with  F#=275  in  R&D  GTDS  was  used).  Again,  a  90  minute  predict  span  followed 
the  fit  span,  with  accuracy  metrics  derived  over  both  the  fit  (DC  RMS)  and  predict  spans  (R&D  GTDS  Ephemeris 
Comparison  Program).  The  results  of  this  testing  are  presented  in  Table  5: 

Table  5.  Simulation  of  Real  World  Drag  Uncertainties,  LACC 


GRAV 

(AOG) 

Zonal 

SP 

M- 

daily 

SP 

TESS 

SP 

h2 

SP 

V 

m-daily 

SP 

Drag 

SRP 

3rd 

Body 

Run 

Time 

DC  RMS 

Total 

Position 

RMS 

(Predict 

Span) 

8x8 

Yes 

Yes 

trunc 

Yes 

No 

H-P 

No 

No 

15.87 

16.29 

150.97 

h2 

6642 

AOG 

sec 

meters 

meters 

(opt. 

-1010 

SPG 

(6  its) 

deck) 

Iszak 

** 

Solve 

CD 

F#=275 

8x8 

Yes 

Yes 

trunc 

Yes 

No 

H-P 

No 

No 

15.73 

24.68 

311.96 

I  2 

J2 

6642 

AOG 

sec 

meters 

meters 

(opt. 

-10  10 

SPG 

(6  its) 

deck) 

v 

Iszak 

" 

** 

Solve 

CD 

F#-150 

Note  the  configurations  listed  in  Table  5  are  designated  as  “opt.  deck”  even  though  they  use  the  Hams-Priester 
atmospheric  model  (the  optimal  configuration  uses  the  Jacchia-roberts  atmospheric  model).  This  designation  was 
given  because  all  geopotential,  partial  derivative,  solve-for,  and  force  model  options  for  the  AOG  and  SPG  match  the 
optimal  configuration  (even  though  the  particular  atmospheric  model  may  be  different). 

The  setup  for  the  LACC  which  provides  an  optimal  balance  between  accuracy  and  computational  speed  can  now 
be  given  (opt.  deck): 

•  8x8  geopotential  (AOG) 

•  J22  (AOG  and  SPG) 

•  Zonal  short  periodics  with  default  settings 

•  M-daily  short  periodics  with  default  settings 

•  Truncated  tesseral  linear  combination  high  frequency  short  periodics  (SPTESSLC  =  6642  -10  10) 

•  Jacchia-Roberts  drag  (AOG  and  SPG) 

•  Solve  CD 

•  SSTESTFL  1  20  0.0 

•  SSTAPGFL  1  00  1.0  6.00.0 

Sample  input  decks  for  the  LACC  are  given  in  the  appendices. 
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3.0  Medium  Altitude  Circular  Case  (MACC) 

A  Landsat  type  orbit  was  chosen  to  fulfill  the  medium  altitude  circular  orbit  test  case.  In  the  ETAS  study,  the 
MACC  is  characterized  by  a  SMAH  range  of  575  -  1000  km  and  an  eccentricity  range  from  0.0  to  0.05.  For  this  test 
case,  multiple  runs  were  made  in  which  varying  configurations  of  DSST  were  fit  to  a  Cowell  truth  model  over  a  three 
day  arc.  A  three  day  DSST  EG  with  the  DC  solve-for  parameters  followed  the  fit  span.  Accuracy  metrics  were 
derived  from  comparisons  between  the  DSST  and  Cowell  trajectories  over  both  the  fit  (DC  RMS)  and  predict  spans 
(from  the  R&D  GTDS  Ephemeris  Comparison  Program).  ORB1  output  was  produced  every  450  seconds  (a  45 
minute  time  interval  was  used  in  the  R&D  GTDS  Ephemeris  Comparison  Program). 

The  osculating  elements  required  for  the  Cowell  truth  EG  were  built  from  a  one  day  DSST  run  with  the 
following  inputs: 

Table  6.  DSST  Specifics  to  Build  Osculating  Elements,  MACC 


Epoch  Date  YYMMDD  HHMMSS.S 

820223  000000.0 

End  Date  YYMMDD  HHMMSS.S 

820224  000000.0 

Semimajor  Axis  (mean) 

7077.787  km 

Eccentricity  (mean) 

0.011542 

Inclination  (mean) 

98.250452  deg 

Longitude  of  Ascending  Node  (mean) 

158.15349  deg 

Argument  of  Perigee  (mean) 

89.4  deg 

Mean  Anomaly  (mean) 

312.90205  deg 

Input  Frame 

True  of  Date 

Integration  and  Output  Frame 

Earth  Equator  and 

Equinox  of  1950 
(Mean  of  1950) 

Integrator 

4th  Order  Runge-Kutta 

DSST  Step  Size 

43200.0  sec 

Gravity  Field 

8x8  GEM10B 

Lunar  /  Solar  Point  Masses 

Yes 

Drag  Model 

No 

Solar  Radiation  Pressure 

No 

The  osculating  elements  which  resulted  from  this  one  day  DSST  run  were  then  used  as  inputs  for  a  six  day 
Cowell  EG.  The  specifics  for  this  Cowell  EG  are  given  in  Table  7: 

Table  7.  Cowell  Truth  EG  Inputs,  MACC 


Epoch  Date  YYMMDD  HHMMSS.S 

820224  000000.0 

End  Date  YYMMDD  HHMMSS.S 

820302  000000.0 
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Table  7.  Cowell  Truth  EG  Inputs,  MACC  (Continued) 


Input,  Integration,  and  Output  Frame 

Earth  Equator  and 

Equinox  of  1950 
(Mean  of  1950) 

SST 

Integrator 

12th  Order  Summed  Cowell/ 
Adams  Predict-Partially  Correct 

Step  Size 

60.0  sec 

Gravity  Field 

21x21  GEM10B 

Lunar  /  Solar  Point  Masses 

Yes 

Drag  Model 

Jacchia-Roberts 

CD  =  2.0 

Solar  Radiation  Pressure 

No 

Average  Spacecraft  Cross-Sectional  Area 

1.0  m2 

Spacecraft  Mass 

100  kg 

i - — - — - - - 

The  same  initial  guess  of  the  state  vector  and  observation  noise  standard  deviation  were  used  for  all  of  the  D 
runs.  These  values  are  highlighted  in  the  following  table: 

Table  8.  DSST  DC  Specifics,  MACC 

i  - —  —  1 

Epoch  Date  YYMMDD  HHMMSS.S 

820224  000000.0 

Fit  Span  Begin  YYMMDD  HHMMSS.S 

820224  000000.0 

Fit  Span  End  YYMMDD  HHMMSS.S 

820227  000000.0 

Semimajor  Axis  (mean) 

7077.8  km 

Eccentricity  (mean) 

0.0011 

Inclination  (mean) 

98.2  deg 

Longitude  of  Ascending  Node  (mean) 

158.1  deg 

Argument  of  Perigee  (mean) 

89.4  deg 

Mean  Anomaly  (mean) 

176.0  deg 

DSST  Input  Frame 

True  of  Date 

DSST  Integration  and  Output  Frame 

Earth  Equator  and 
Equinox  of  1 950 
(Mean  of  1950) 

DSST  Step  Size 

43200.0  sec 

Geopotential  Model 

GEM10B 

Position  Standard  Deviation 

100  meters 

Velocity  Standard  Deviation 

10  cm/sec 
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In  the  DSST  DC  runs,  varying  configurations  corresponding  to  different  perturbation  modeling  were  used  to 
determine  the  optimal  balance  of  accuracy  and  computational  speed.  Results  of  the  DC’s  and  subsequent  EG  s  for  the 
MACC  are  given  in  Table  9;  again,  it  is  first  desirable  to  describe  some  of  the  notation  used  in  this  table  (and  other 
tables  in  this  section): 

•  GRAV  refers  to  the  geopotential  modeling  in  the  AOG. 

•  RESONPRD  is  a  GTDS  keyword  for  the  input  processor  to  override  the  default  border  for  the  modeling  of 
tesseral  resonance  in  the  AOG  (if  the  period  of  the  resonance  is  less  than  10  days  in  the  default  configuration, 
the  resonance  is  modelled  in  the  SPG  as  a  tesseral  linear  combination  high  frequency  term);  by  specifying  a 
value  with  this  keyword,  a  new  border  can  be  set  to  force  the  modeling  of  resonance  into  the  AOG  (i.e.,  for 
shallow  resonance). 

•  TESS  refers  to  the  tesseral  linear  combination  high  frequency  short  periodic  terms. 

•  SRP  is  an  abbreviation  for  solar  radiation  pressure. 

•  J-R  is  an  abbreviation  for  the  Jacchia-Roberts  atmospheric  density  model. 

•  H-P  is  an  abbreviation  for  the  Harris-Priester  atmospheric  density  model. 

•  L-S  is  an  abbreviation  for  lunar  /  solar  third  body  point  mass  effects. 

•  trunc  refers  to  truncated  modeling;  specifically,  the  maximum  power  of  the  eccentricity  was  set  to  2  and  the 
maximum  frequency  in  true  longitude  was  set  to  11  in  Fourier  coefficients  for  these  short  periodic 
expansions33. 

•  Total  Position  RMS  is  derived  over  the  three  day  predict  span  only  (from  the  R&D  GTDS  Ephemeris 
Comparison  Program). 

•  All  cases  which  attempt  to  solve  for  CD  include  first  order  drag  effects  in  the  AOG  partials  (if  CD  is  not  solved 
for,  then  all  drag  partials  in  the  AOG  are  shut  off).  All  cases  include  J2  partials  analytically  in  the  AOG.  No 
other  AOG  or  SPG  partials  were  included  in  any  of  the  runs. 

•  Most  cases  have  an  SSTESTFL  card  of  1  3;  this  setting  includes  all  desired  AOG  partials  in  the  A  matrix  (J2 
partials  are  computed  analytically;  all  other  partial  derivatives  are  computed  via  finite  differences).  In  addition, 
this  setting  includes  desired  AOG  partial  derivatives  for  non-conservative  perturbations  in  the  D  matrix  via 
finite  differences  (these  partials  were  also  included  in  the  A  matrix).  Cases  noted  by  **  have  an  SSTESTFL 
card  of  1  2;  this  setting  puts  only  J2  AOG  partials  into  the  A  matrix  (computed  analytically).  All  other  desired 
AOG  partials  for  non-conservative  perturbations  are  included  in  the  D  matrix  via  finite  differences  (with  this 
setting,  third  body  and  J22  AOG  partial  derivatives,  if  desired,  would  not  be  included  anywhere;  if  they  are 
desired,  the  SSTESTFL  setting  of  1  3  should  be  used). 


Table  9.  DSST  DC  and  Subsequent  EG  Results,  MACC 


GRAV 

(AOG) 

Zonal 

SP 

M- 

daily 

SP 

TESS 

SP 

h2 

SP 

J 2/ 

m-daily 

SP 

Drag 

SRP 

3rd 

Body 

Run 

Time 

DC  RMS 

Total 

Position 

RMS 

(Predict 

Span) 

21x21 

Yes 

Yes 

Yes 

Yes 

Yes 

J-R 

No 

L-S 

232.83 

23.50 

248.76 

J22 

AOG 

AOG 

sec 

meters 

meters 

SPG 

SPG 

(6  its) 

21x21 

Yes 

Yes 

Yes 

Yes 

Yes 

J-R 

No 

L-S 

311.6 

9.38 

260.40 

h2 

AOG 

AOG 

sec 

meters 

meters 

RESON 

SPG 

SPG 

(9  its) 

PRD  = 

1  day 
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Table  9.  DSST  DC  and  Subsequent  EG  Results,  MACC  (Continued) 


GRAV 

(AOG) 

Zonal 

SP 

M- 

daily 

SP 

TESS 

SP 

h1 

SP 

V 

m-daily 

SP 

Drag 

SRP 

3rd 

Body 

Run 

Time 

DC  RMS 

Total 

Position 

RMS 

(Predict 

Span) 

Yes 

Yes 

Yes 

Yes 

Yes 

J-R 

AOG 

SPG 

Solve 

CD 

No 

L-S 

AOG 

SPG 

674.48 

sec 

13.18 

meters 

(23  its) 

212.13 

meters 

21x21 

h2 

RESON 
PRD  — 

1  day 

Yes 

Yes 

Yes 

Yes 

Yes 

J-R 

AOG 

SPG 

Solve 

CD 

No 

L-S 

AOG 

SPG 

317.14 

sec 

9.20 

meters 

(9  its) 

236.57 

meters 

21x21 

h2 

RESON 
PRD  = 

2  days 

Yes 

Yes 

Yes 

Yes 

Yes 

J-R 

AOG 

SPG 

No 

L-S 

AOG 

SPG 

310.44 

sec 

8.28 

meters 

(9  its) 

259.63 

meters 

21x21 

T  2 

J2 

RESON 
PRD  = 

2  days 

Yes 

Yes 

Yes 

Yes 

•Yes 

J-R 

AOG 

SPG 

Solve 

cD 

No 

L-S 

AOG 

SPG 

393.10 

sec 

7.75 

meters 

(12  its) 

241.01 

meters 

i 

Yes 

Yes 

Yes 

Yes 

Yes 

J-R 

AOG 

SPG 

No 

L-S 

AOG 

SPG 

47.29 

sec 

123.48 

meters 

(7  its) 

342.54 

meters 

1 

Yes 

Yes 

Yes 

Yes 

Yes 

J-R 

AOG 

SPG 

Solve 

CD 

No 

L-S 

AOG 

SPG 

52.98 

sec 

117.47 

meters 

(8  its) 

386.02 

meters 

8x8 

T  2 

J2 

Yes 

Yes 

Yes 

Yes 

Yes 

J-R 

AOG 

SPG 

Solve 

CD 

No 

L-S 

AOG 

SPG 

72.72 

sec 

71.15 

meters 

(7  its) 

261.57 

meters 
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Table  9.  DSST  DC  and  Subsequent  EG  Results,  MACC  (Continued) 


GRAV 

(AOG) 

Zonal 

SP 

M- 

daily 

SP 

TESS 

SP 

h2 

SP 

V 

m-daily 

SP 

Drag 

SRP 

3rd 

Body 

Run 

Time 

DC  RMS 

Total 

Position 

RMS 

(Predict 

Span) 

8x8 

T  * 

J2 

Yes 

Yes 

No 

Yes 

Yes 

J-R 

AOG 

SPG 

Solve 

CD 

No 

L-S 

AOG 

SPG 

46.68 

sec 

92.70 

meters 

(7  its) 

279.60 

meters 

8x8 

h1 

Yes 

Yes 

No 

Yes 

No 

J-R 

AOG 

SPG 

Solve 

CD 

No 

L-S 

AOG 

SPG 

41.47 

sec 

92.96 

meters 

(7  its) 

277.11 

meters 

8x8 

h2 

Yes 

Yes 

No 

Yes 

No 

J-R 

AOG 

Solve 

CD 

No 

L-S 

AOG 

SPG 

40.26 

sec 

92.95 

meters 

(7  its) 

276.77 

meters 

8x8 

h2 

trunc 

8211 

trunc 

882 

No 

Yes 

No 

J-R 

AOG 

Solve 

CD 

No 

L-S 

AOG 

SPG 

39.33 

sec 

92.95 

meters 

(7  its) 

276.77 

meters 

8x8 

h2 

trunc 

8211 

trunc 

882 

No 

Yes 

No 

J-R 

AOG 

Solve 

CD 

No 

L-S 

AOG 

38.06 

sec 

92.92 

meters 

(7  its) 

276.70 

meters 

8x8 

h2 

trunc 

8211 

trunc 

882 

No 

No 

No 

J-R 

AOG 

No 

L-S 

AOG 

33.83 

sec 

97.50 

meters 

(7  its) 

259.80 

meters 

8x8 

h2 

trunc 

8211 

trunc 

882 

No 

No 

No 

J-R 

AOG 

Solve 

CD 

No 

L-S 

AOG 

37.08 

sec 

93.06 

meters 

(7  its) 

278.26 

meters 

8x8 

h2 
*  * 

trunc 

8211 

trunc 

882 

No 

No 

No 

J-R 

AOG 

Solve 

CD 

No 

■ 

34.28 

sec 

93.06 

meters 

(7  its) 

278.30 

meters 
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Table  9.  DSST  DC  and  Subsequent  EG  Results,  MACC  (Continued) 


GRAV 

(AOG) 

Zonal 

SP 

M- 

daily 

SP 

TESS 

SP 

h2 

SP 

V 

m-daily 

SP 

Drag 

SRP 

3rd 

Body 

Run 

Time 

DC  RMS 

Total 

Position 

RMS 

(Predict 

Span) 

21x21 

h2 
*  * 

RJESON 

PRD  = 

1  day 

(opt. 

deck) 

trunc 

8211 

trunc 

882 

No 

No 

No 

J-R 

AOG 

Solve 

CD 

No 

L-S 

AOG 

35.31 

sec 

76.51 

meters 

(7  its) 

196.85 

meters 

Spa 

ce  Comn 

land’s  SGP4  (solve  for  drag  term  ndot/2) 

28.67 

sec 

581.68 

meters 

(5  its) 

1412.3 

meters 

Space  Command’s  SGP4  (do  not  solve  for  drag  term  ndot/2) 

B*  =  0.0001  er  1 

36.03 

sec 

905.98 

meters 

(12  its) 

17876.0 

meters 

The  results  in  this  table  indicate  that  when  fitting  an  optimized  version  of  DSST  to  simulated  data  for  the  MACC, 
solving  for  Cp  appears  to  make  the  results  worse.  However,  solving  for  Cjj  in  the  “full-up”  DSST  cases  leads  to 
results  which  are  more  accurate  than  in  the  analogous  case  in  which  Cp  is  not  solved  for.  This  phenomena  can  be 
accredited  to  the  decreased  modeling  in  the  optimized  DSST  configuration.  Specifically,  error  characteristics 
resulting  from  perturbations  other  than  atmospheric  drag  may  bias  the  drag  effect.  For  example,  the  drag  model  in 
DSST  may  attempt  to  account  for  error  signatures  in  the  semimajor  axis  which  stem  from  un-modeled  tesseral 
resonance;  the  reduced  modeling  makes  it  more  difficult  for  the  OD  system  to  separate  which  source  perturbing 
effects  are  coming  from.  If  Cp  is  not  solved  for  in  these  cases,  its  value  defaults  to  the  same  value  used  in  the 
generation  of  the  truth  ephemeris;  hence,  better  results  can  be  expected  in  these  controlled  experiments  if  CD  is  not 
solved  for.  In  real  world  analysis,  the  actual  value  of  CD  may  not  be  known  very  accurately,  which  necessitates  a 
solution  for  Cp. 

Clearly,  the  modeling  of  the  tesseral  short  periodics  can  be  truncated  to  greatly  reduce  the  computer  run  time.  In 
fact,  for  the  8x8  case  in  which  the  tesseral  short  periodics  have  been  turned  off,  a  36%  savings  in  computer  run  time 
is  gained  at  the  expense  of  an  approximate  20  meter  fit  and  predict  difference  (as  compare  to  the  8x8  case  with  a  full 
tesseral  short  periodic  model).  Similarly,  the  zonal  and  m-daily  short  periodic  modeling  can  be  slightly  truncated  to 
benefit  the  computer  run  time  without  a  significant  loss  of  accuracy.  Third  body  and  drag  short  periodics  can  be 
removed  to  enhance  run  time  with  no  degradation  in  accuracy. 

Both  the  J22  and  the  J2  /  m-daily  model  serve  to  provide  accuracy  closure  with  Cowell,  but  at  an  expense  of 
computation  time.  For  reasons  of  simplicity  and  computational  time  savings,  these  models  can  be  neglected  in  an 
optimized  DSST  configuration. 


j 


i 
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The  period  of  resonance  for  the  14th  order  is  less  than  two  days  (approximately  1.84  days),  while  the  period  of 
resonance  for  the  15th  order  is  less  than  three  days  (approximately  2.18  days).  Therefore,  running  R&D  GTDS 
without  a  RESONPRD  card  to  change  the  default  resonance  border  from  ten  days  would  force  the  shallow  resonance 
terms  at  these  orders  into  the  tesseral  linear  combination  model  of  the  SPG.  A  RESONPRD  card  setting  the  resonance 
border  at  one  day  forces  the  shallow  resonance  terms  at  both  the  14th  and  15th  orders  into  the  AOG.  Using 
RESONPRD  with  a  value  of  two  days  only  puts  the  15th  order  shallow  resonance  into  the  AOG.  As  is  shown  in  the 
results  using  a  “full”  21x21  geopotential  field,  modeling  resonance  in  the  AOG  provides  for  modest  improvements  in 
fit  span  accuracy  over  modeling  resonance  in  the  SPG.  However,  improvements  in  predict  span  accuracy  are  not 
necessarily  gained.  With  the  optimized  version  of  DSST,  almost  100  meters  worth  of  improvement  in  predict 
accuracy  is  gained  by  including  the  resonant  terms  at  the  14th  and  15th  orders  in  the  AOG  (as  opposed  to  neglecting 
these  terms  with  an  optimized  version  of  DSST  using  an  8x8  base  gravity  field.  It  should  be  noted  some  of  this 
improvement  can  be  attributed  to  the  increased  AOG  zonal  modeling  which  results  from  the  method  used  to  include 
the  shallow  resonance  terms-refer  to  the  following  paragraph  for  details  of  this  method).  As  a  general  rule  of  thumb, 
it  is  probably  better  to  have  resonance  modeled  in  the  AOG  due  to  the  nature  of  the  theory.  In  this  manner,  resonance 
is  captured  in  the  mean  equations  of  motion  and  handled  numerically  (which  tends  to  provide  better  treatment  of  the 
non-linearities  and  large  magnitudes  associated  with  resonance). 


It  is  of  particular  importance  to  use  extreme  care  when  running  an  optimized  version  of  the  DSST.  For  example, 
if  the  tesseral  linear  combination  high  frequency  short  periodics  have  been  truncated  or  removed  for  an  orbit  in  which 
shallow  resonance  occurs,  resonance  terms  will  be  totally  neglected  (assuming  a  resonant  period  <10  days  and  no 
adjustment  with  the  RESONPRD  card).  In  these  cases,  the  user  must  experiment  with  the  AUTOFORC  and  RESNM 
cards  or  the  RESONPRD  card  and  hand  calculation  of  the  resonant  periods  to  ensure  the  appropriate  resonance 
modeling  occurs  (which  can  become  quite  tedious).  As  an  alternative,  the  user  can  specify  a  21x21  gravity  field  with 
truncated  or  no  tesseral  linear  combination  short  periodics  and  a  RESONPRD  of  one  day.  This  would  include  all 
shallow  resonant  terms  (and  obviously  deep  resonant  terms)  with  a  period  greater  than  one  day  into  the  AOG  at  the 
expense  of  additional  secular  and  long  period  zonal  calculations  in  the  AOG.  This  configuration,  which  is  used  in  the 
MACC  optimal  deck,  is  much  simpler  than  experimenting  with  the  AUTOFORC  and  RESNM  cards  or  hand 
calculation  of  resonant  period  (which  is  not  practical  for  operational  scenarios  with  multiple  satellites).  For  example, 
an  8x8  base  gravity  field  could  have  been  established  for  the  MACC  along  with  AUTOFORC  and  RESNM  cards  to 
include  all  the  resonant  terms  at  the  14th  and  15th  order  (which  requires  1  AUTOFORC  card  and  15  RESNM  cards). 
This  configuration  would  not  introduce  any  secular  and  long  period  zonal  calculations  in  the  AOG  beyond  the  8x8 
base  gravity  field.  However,  for  only  a  3%  increase  in  computer  run  time  (35.31  seconds  versus  34.28  seconds),  the 
21x21  configuration  with  a  RESONPRD  of  one  day  can  be  used  (which  does  introduce  secular  and  long  period  zonal 
calculations  in  the  AOG  beyond  the  8x8  base  gravity  field  up  to  the  21x21  limit)  and  remove  the  introduction  of  16 
R&D  GTDS  keyword  cards.  This  example  shows  the  generality  and  efficiency  with  which  the  recursions  handle  the 
calculation  of  zonal  effects  in  the  AOG  (it  should  be  noted  the  same  recursion  concept  is  used  for  all  portions  of  the 
geopotential  in  the  AOG  and  SPG). 

This  general  handling  of  resonant  terms  becomes  extremely  useful  when  considering  multiple  satellites  under  the 
broad  definition  of  the  MACC.  Instead  of  introducing  costly  “man-in-the-loop”  calculations,  this  general 
configuration  efficiently  captures  resonant  terms  (shallow  and  deep)  for  all  orbits  in  the  MACC  (up  to  the  21x21 
limit;  note  that  the  same  general  concept  can  be  used  with  the  50x50  version  of  R&D  GTDS). 

The  initial  drag  results  led  to  some  further  processing  with  the  Harris  Priester  drag  model.  In  contrast  to  the 
Jacchia-Roberts  drag  model,  the  Hanis-Priester  model  can  be  considered  “static”  in  nature.  The  “dynamic”  Jacchia- 
Roberts  model  implements  “observed”  or  predicted  values  for  drag-related  parameters  (solar  flux  and  geomagnetic 
data).  The  Harris-Priester  model  implements  established  “look-up”  tables  to  obtain  drag-related  parameters  (altitude 
versus  density).  In  this  manner,  the  Harris-Priester  model  presents  a  “smoother”  or  more  stable  picture  of  the 
atmosphere.  It  should  be  noted  various  Harris-Priester  density  tables  have  been  constructed  to  correspond  to  different 
levels  of  activity  in  the  atmosphere. 
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Results  for  a  full-up  Harris-Priester  DSST  fit  to  a  Harris-Priester  Cowell  truth  model  are  given  in  Table  10  (all 
test  procedures  for  items  other  than  atmospheric  drag  are  identical  to  the  other  test  procedures  described  in  this 
section): 

Table  10.  DSST  Harris-Priester  DC  and  Subsequent  EG  Results  I,  MACC 


GRAV 

(AOG) 

Zonal 

SP 

M- 

daily 

SP 

TESS 

SP 

h2 

SP 

V 

m -daily 

SP 

Drag 

SRP 

3rd 

Body 

Run 

Time 

DC  RMS 

Total 

Position 

RMS 

(Predict 

Span) 

BBSIB 

Yes 

Yes 

Yes 

Yes 

Yes 

H-P 

No 

L-S 

230.14 

16.20 

21.218 

■SB 

AOG 

AOG 

sec 

meters 

meters 

SPG 

SPG 

(6  its) 

■ 

F#=150 

Yes 

Yes 

Yes 

Yes 

Yes 

H-P 

No 

L-S 

283.36 

2.97 

4.7777 

h2 

AOG 

AOG 

sec 

meters 

meters 

RESON 
PRD  = 

1  day 

SPG 

F#=150 

SPG 

(8  its) 

Yes 

Yes 

Yes 

Yes 

Yes 

H-P 

No 

L-S 

358.83 

7.64 

280.20 

AOG 

AOG 

sec 

meters 

meters 
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SPG 

SPG 

(Hits) 

■ 

'• '  '' 

i 

i 

Solve 

cD 

■ 

F#=150 

21x21 

Yes 

Yes 

Yes 

Yes 

Yes 

H-P 

No 

L-S 

285.00 

2.96 

9.049 

h2 

AOG 

AOG 

sec 

meters 

meters 

RESON 
PRD  = 

1  day 

SPG 

Solve 

cD 

F#=150 

SPG 

(8  its) 

As  expected,  the  results  in  this  table  show  the  static  Harris-Priester  model  provides  for  a  tighter  DSST  fit  to  a 
Cowell  truth  model  (stressing  the  difficulties  which  arise  due  to  the  dynamic,  “noisy”  nature  of  the  atmosphere;  DCs 
with  the  Jacchia-Roberts  atmospheric  model  did  not  provide  the  same  tight  closure  as  with  the  “smooth”  Harris 
Priester  model.  Real  world  analysis  of  atmospheric  effects  proves  quite  challenging  due  to  the  noisy  and  uncertain 
nature  of  atmospheric  conditions,  especially  for  orbital  predictions  into  the  future).  Of  particular  importance  in  these 
results  is  the  accuracy  improvement  gained  by  modeling  resonance  in  the  AOG  in  cases  in  which  Cp  is  solved  for . 


Some  testing  was  also  done  to  try  and  simulate  real-world  uncertainties  with  atmospheric  drag.  In  these  tests,  a 
Jacchia-Roberts  drag  model  was  used  to  generate  a  Cowell  truth  EG.  Then,  a  Harris-Priester  DSST  configuration  was 
fit  over  a  three  day  arc  to  the  Jacchia-Roberts  Cowell  truth  trajectory.  Again,  a  three  day  predict  span  followed  the  fit 
span,  with  accuracy  metrics  derived  over  both  the  fit  (DC  RMS)  and  predict  spans  (R&D  GTDS  Ephemeris 
Comparison  Program).  The  results  of  this  testing  are  presented  in  Table  11 : 
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Table  11.  Simulation  of  Real  World  Drag  Uncertainties,  MACC 


GRAV 

(AOG) 

Zonal 

SP 

M- 

daily 

SP 

TESS 

SP 

J l 

SP 

V 

m -daily 

SP 

Drag 

SRP 

3rd 

Body 

Run 

Time 

DC  RMS 

Total 
Position 
RMS 
(Fit  -H 
Predict) 

21x21 

ti 
*  * 

RESON 

PRD  = 

1  day 

trunc 

8211 

trunc 

882 

No 

No 

No 

H-P 

AOG 

F#=275 

No 

L-S 

AOG 

40.37 

sec 

764.30 

meters 

(12  its) 

15754.0 

meters 

21x21 

h1 
*  * 

RESON 

PRD  — 

1  day 

(opt. 

deck) 

trunc 

8211 

trunc 

882 

No 

No 

No 

j 

H-P 

AOG 

Solve 

CD 

F#=275 

No 

L-S 

AOG 

35.04 

sec 

78.98 

meters 

(7  its) 

469.52 

meters 

21x21 

h1 
*  * 

RESON 

PRD  ~ 

1  day 

(opt. 

deck) 

trunc 

8211 

trunc 

882 

No 

No 

No 

H-P 

AOG 

Solve 

CD 

F#=150 

No 

L-S 

AOG 

34.77 

sec 

80.06 

meters 

(7  its) 

636.23 

meters 

Note  the  configurations  listed  in  Table  1 1  are  designated  as  “opt  deck”  even  though  they  use  the  Hams-Priester 
atmospheric  model  (the  optimal  configuration  uses  the  Jacchia-roberts  atmospheric  model).  This  designation  was 
given  because  all  geopotential,  partial  derivative,  solve-for,  and  force  model  options  for  the  AO G  and  SPG  match  the 
optimal  configuration  (even  though  the  particular  atmospheric  model  may  be  different). 

The  setup  for  the  MACC  which  provides  an  optimal  balance  between  accuracy  and  computational  speed  can  now 
be  given: 

•  21x21  geopotential  (AOG) 

•  RESONPRD  equal  to  one  day  (86400.0  seconds) 

•  J22  (AOG) 

•  Truncated  zonal  short  periodic  modeling  (the  maximum  power  of  the  eccentricity  was  set  to  2  and  the 
maximum  frequency  in  true  longitude  was  set  to  1 1  in  Fourier  coefficients  for  these  short  periodic  expansions) 
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•  Truncated  tesseral  m-daily  short  periodic  modeling  (maximum  power  of  the  eccentricity  was  set  to  2  in  Fourier 
Coefficients  for  these  short  periodic  expansions) 

•  Jacchia-Roberts  drag  (AOG) 

•  Solve  CD 

•  Lunar  /  solar  third  body  point  mass  effects  (AOG) 

•  SSTESTFL  120  0.0 

•  SSTAPGFL  100  1.0  0.0  0.0 

Sample  input  decks  for  the  MACC  are  given  in  the  appendices. 

4.0  High  Altitude  Circular  Case  (HACC) 

The  Explorer  27  (BE-C)  orbit  was  chosen  to  fulfill  the  high  altitude  circular  test  case.  In  the  ETAS  study,  the 
HACC  is  characterized  by  a  SMAH  range  of  1 000  -  2500  km  and  an  eccentricity  range  from  0.0  to  0.05.  For  this  test 
case,  multiple  runs  were  made  in  which  varying  configurations  of  DSST  were  fit  to  a  Cowell  truth  model  over  a  two 
day  arc.  A  three  day  DSST  EG  with  the  DC  solve-for  parameters  followed  the  fit  span.  Accuracy  metrics  were 
derived  from  comparisons  between  the  DSST  and  Cowell  trajectories  over  both  the  fit  (DC  RMS)  and  predict  spans 
(from  the  R&D  GTDS  Ephemeris  Comparison  Program).  ORB1  output  was  produced  every  60  seconds  for  the 
Cowell  truth  generation  (a  45  minute  time  interval  was  used  in  the  R&D  GTDS  Ephemeris  Comparison  Program). 

The  osculating  elements  required  for  the  Cowell  truth  EG  were  derived  from  a  five  day  Cowell  fit  to  an  SGP4 
.  ephemeris.  The  two  card  element  set  used  to  generate  the  SGP4  ephemeris  is  given  in  Table  12: 

Table  12.  Explorer  Two-Card  Element  Set,  HACC,  September  1994  SATCAT 


Epoch  Date  YYMMDD  HHMMSS.SSSS 

940826  074853.5101 

Mean  Motion 

13.37418275  rev/day 

Eccentricity 

0.0247108 

Inclination 

41.1904  deg 

Longitude  of  Ascending  Node 

349.3346  deg 

Argument  of  Perigee 

66.5461  deg 

Mean  Anomaly 

296.1065  deg 

B* 

0.000091  er'1 

ORB1  Output  Frequency 

Every  450.0  Seconds 

The  specifics  for  the  Cowell  DC  and  subsequent  EG  are  given  in  Table  13: 

Table  13.  Cowell  DC  and  Subsequent  EG  Inputs,  HACC 


Epoch  Date  YYMMDD  HHMMSS.SSSS 

940826  074853.5101 

Fit  Begin  Date  YYMMDD  HHMMSS.SSSS 

940826  074854.0000 

Fit  End  Date  YYMMDD  HHMMSS.SSSS 

940831  074854.0000 

EG  End  Date  YYMMDD  HHMMSS.SSSS 

940903  074854.0000 
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Table  13.  Cowell  DC  and  Subsequent  EG  Inputs,  HACC  (Continued) 


Input,  Integration,  and  Output  Frame 

Earth  Equator  and 

Equinox  of  1950 
(Mean  of  1950) 

Integrator 

12th  Order  Summed  Cowell/ 
Adams  Predict-Partially  Correct 

Step  Size 

60.0  sec 

Gravity  Field 

21x21  JGM2 

Lunar  /  Solar  Point  Masses 

Yes 

Drag  Model 

Jacchia-Roberts-Schatten 

(Solve  pj) 

Solar  Radiation  Pressure 

Yes 

(Solve  Cr) 

Average  Spacecraft  Cross-Sectional  Area 

1.0  m2 

Spacecraft  Mass 

100  kg 

Position  Standard  Deviation 

500  meters 

Velocity  Standard  Deviation 

50  cm/sec 

ORB1  Output  Frequency 

Every  60.0  Seconds 

The  same  initial  guess  of  the  state  vector  and  observation  noise  standard  deviation  were  used  for  all  of  the  DSST 
DC  runs.  These  values  are  highlighted  in  Table  14: 

Table  14.  DSST  DC  Specifics,  HACC 


Epoch  Date  YYMMDD  HHMMSS.SSSS 

940826  074853.5101 

Fit  Span  Begin  YYMMDD  HHMMSS.S 

940826  074854.0 

Fit  Span  End  YYMMDD  HHMMSS.S 

940828  074854.0 

Semimajor  Axis  (mean) 

7498.0  km 

Eccentricity  (mean) 

0.0247108 

Inclination  (mean) 

41.1904  deg 

Longitude  of  Ascending  Node  (mean) 

349.3346  deg 

Argument  of  Perigee  (mean) 

66.5461  deg 

Mean  Anomaly  (mean) 

296.1065  deg 

DSST  Input  Frame 

True  of  Date 

DSST  Integration  and  Output  Frame 

Earth  Equator  and 

Equinox  of  1950 
(Mean  of  1950) 

23 


Table  14.  DSST  DC  Specifics,  HACC  (Continued) 


DSST  Step  Size 

43200.0  sec 

Geopotential  Model 

JGM2 

Position  Standard  Deviation 

100  meters 

Velocity  Standard  Deviation 

10  cm/sec 

ORBl  Output  Frequency 

Eveiy  450.0  Seconds 

In  the  DSST  DC  runs,  varying  configurations  corresponding  to  different  perturbation  modeling  were  used  to 
determine  the  optimal  balance  of  accuracy  and  computational  speed.  Results  of  the  DC’s  and  subsequent  EG’s  for  the 
HACC  are  given  in  Table  1 5;  however,  it  is  first  desirable  to  highlight  a  few  notes  concerning  the  test  protocol  for  the 
HACC: 

•  Total  Position  RMS  is  derived  over  the  three  day  predict  span  only  (R&D  GTDS  Ephemeris  Comparison 
Program). 

•  R&D  GTDS  default  configuration  settings  are  used  unless  specifically  noted. 

•  RES  refers  to  the  inclusion  of  the  13th  and  14th  order  shallow  resonance  terms  with  the  AUTOFORC  and 
RESNM  keyword  cards. 

•  “trunc”  refers  to  truncated  modeling;  modeling  descriptions  corresponding  to  the  R&D  GTDS  Semianalytic 
Satellite  Theory  input  processor  are  specified  where  needed.  Additional  details  can  be  found  in  the  R&D 
GTDS  Semianalytic  Theory  Input  Processor  keyword  description  booklet33. 

•  All  cases  which  attempt  to  solve  for  CD  include  first  order  drag  effects  in  the  AOG  partials  (if  CD  is  not  solved 
for,  then  all  drag  partials  in  the  AOG  are  shut  off).  Similarly,  all  cases  which  attempt  to  solve  for  CR  include 
solar  radiation  pressure  effects  in  the  AOG  partials  (if  CR  is  not  solved  for,  then  all  solar  radiation  pressure 
partials  in  the  AOG  are  shut  off).  All  cases  include  J2  partials  analytically  in  the  AOG.  No  other  AOG  or  SPG 
partials  were  included  in  any  of  the  runs. 

•  Most  cases  have  an  SSTESTFL  card  of  1  3;  this  setting  includes  all  desired  AOG  partials  in  the  A  matrix  (J2 
partials  are  computed  analytically;  all  other  partial  derivatives  are  computed  via  finite  differences).  In  addition, 
this  setting  includes  desired  AOG  partial  derivatives  for  non-conservative  perturbations  in  the  D  matrix  via 
finite  differences  (these  partials  were  also  included  in  the  A  matrix).  Cases  noted  by  **  have  an  SSTESTFL 
card  of  1  2;  this  setting  puts  only  J2  AOG  partials  into  the  A  matrix  (computed  analytically).  All  other  desired 
AOG  partials  for  non-conservative  perturbations  are  included  in  the  D  matrix  via  finite  differences  (with  this 
setting,  third  body  and  AOG  partial  derivatives,  if  desired,  would  not  be  included  anywhere;  if  they  are 
desired,  the  SSTESTFL  setting  of  1  3  should  be  used). 
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Table  15.  DSST  DC  and  Subsequent  EG  Results,  HACC 


GRAV 

(AOG) 

Zonal 

SP 

M- 

daily 

SP 

TESS 

SP 

h2 

SP 

V 

m -daily 

SP 

Drag 

SRP 

3rd 

Body 

Run 

Time 

DC  RMS 

Total 

Position 

RMS 

(Predict 

Span) 

Yes 

Yes 

Yes 

Yes 

Yes 

J-R 

AOG 

L-S 

1165.3 

2.77 

74.141 

h2 

AOG 

SPG 

AOG 

sec 

meters 

meters 

RESON 

SPG 

Solve 

SPG 

(11  its) 

PRD  = 

1  day 

Solve 

cD 

CR 

21x21 

Yes 

Yes 

trunc 

Yes 

No 

J-R 

AOG 

L-S 

203.45 

13.93 

76.522 

h2 

6642 

AOG 

SPG 

AOG 

sec 

meters 

meters 

RESON 

-1010 

SPG 

Solve 

SPG 

(6  its) 

PRD  = 

1  day 

Solve 

CD 

Cr 

ESI 

Yes 

Yes 

trunc 

Yes 

No 

J-R 

AOG 

L-S 

156.76 

20.79 

71.134 

h2 

4442 

AOG 

SPG 

AOG 

sec 

meters 

meters 

RESON 

-8  8 

SPG 

Solve 

SPG 

(5  its) 

PRD  = 

1  day 

Solve 

cD 

cR 

21x21 

Yes 

Yes 

No 

Yes 

No 

J-R 

AOG 

L-S 

144.56 

39.88 

66.637 

h2 

AOG 

SPG 

AOG 

sec 

meters 

meters 

RESON 

SPG 

Solve 

SPG 

<6  its) 

PRD  = 

1  day 

Solve 

CD 

CR 

8x8 

Yes 

Yes 

No 

Yes 

No 

J-R 

AOG 

L-S 

125.23 

52.60 

149.71 

h2 

AOG 

SPG 

AOG 

sec 

meters 

meters 

SPG 

Solve 

SPG 

(6  its) 

Solve 

CD 

CR 

8x8 

Yes 

Yes 

No 

Yes 

No 

J-R 

No 

L-S 

103.81 

58.77 

332.22 

h2 

AOG 

AOG 

sec 

meters 

meters 

SPG 

SPG 

(4  its) 

Solve 

CD 
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Table  15.  DSST  DC  and  Subsequent  EG  Results,  HACC  (Continued) 


GRAV 

(AOG) 

Zonal 

SP 

M- 

daily 

SP 

TESS 

SP 

T  2 

J2 

SP 

V 

m -daily 

SP 

Drag 

SRP 

3rd 

Body 

Run 

Time 

DC  RMS 

Total 

Position 

RMS 

(Predict 

Span) 

8x8 

h2 

Yes 

Yes 

No 

Yes 

No 

J-R 

AOG 

SPG 

Solve 

CD 

AOG 

L-S 

AOG 

SPG 

102.21 

sec 

61.05 

meters 

(4  its) 

392.25 

meters 

8x8 

Yes 

Yes 

No 

Yes 

No 

J-R 

AOG 

SPG 

Solve 

CD 

AOG 

Solve 

CR 

L-S 

AOG 

SPG 

102.64 

sec 

53.19 

meters 

(4  its) 

155.03 

meters 

8x8 

T  ^ 

J2 

Yes 

Yes 

No 

Yes 

No 

No 

AOG 

Solve 

CR 

L-S 

AOG 

SPG 

98.03 

sec 

53.38 

meters 

(4  its) 

151.55 

meters 

8x8 

h2 

Yes 

Yes 

No 

Yes 

No 

No 

AOG 

Solve 

CR 

L-S 

AOG 

96.29 

sec 

i 

53.39 

meters 

(4  its) 

151.69 

meters 

8x8  j 

T  2 

J2  ! 

trunc 

8211 

trunc 

882 

9 

No 

Yes 

No 

No 

AOG 

Solve 

CR 

L-S 

AOG 

94.58 

sec 

53.39 

meters 

(4  its) 

151.65 

meters 

8x8 

T  2 

J2 

trunc 

8211 

trunc 

882 

No 

No 

No 

No 

AOG 

Solve 

Cr 

L-S 

AOG 

92.89 

sec 

53.59 

meters 

(4  its) 

151.41 

meters 

6x6 

ti 

trunc 

6211 

trunc 

662 

No 

No 

No 

No 

AOG 

Solve 

CR 

L-S 

AOG 

91.18 

sec 

80.34 

meters 

(4  its) 

285.94 

meters 

8x8 

h1 

RES 

trunc 

8211 

trunc 

882 

No 

No 

No 

No 

AOG 

Solve 

CR 

L-S 

AOG 

92.39 

sec 

48.53 

meters 

(4  its) 

111.23 

meters 

trunc 

trunc 

No 

No 

No 

No 

AOG 

L-S 

96.83 

48.71 

81.619 

T  2 

J2 

RESON 

PRD  ~ 

1  day 

82  11  ! 

882 

i 

Solve 

CR 

AOG 

sec 

meters 

(4  its) 

meters 
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Table  15-  DSST  DC  and  Subsequent  EG  Results,  HACC  (Continued) 


GRAV 

(AOG) 

Zonal 

SP 

M- 

daily 

SP 

TESS 

SP 

h1 

SP 

V 

m-daily 

SP 

Drag 

SRP 

3rd 

Body 

Run 

Time 

DC  RMS 

Total 

Position 

RMS 

(Predict 

Span) 

21x21 

T  ^ 

J2 

RESON 

PRD  = 

1  day 

He  He 

1/5 

day 

step 

size 

trunc 

8211 

trunc 

882 

No 

No 

No 

No 

AOG 

Solve 

Cr 

L-S 

AOG 

97.10 

sec 

48.71 

meters 

(4  its) 

81.693 

meters 

21x21 

h2 

RESON 

PRD  — 

1  day 

*  * 

(opt. 

deck) 

trunc 

8211 

trunc 

i 

882 

No 

No 

No 

No 

AOG 

Solve 

CR 

L-S 

AOG 

93.12 

sec 

48.71 

meters 

(4  its) 

81.619 

meters 

Space  Command’s  SGP4  (solve  for  drag  term  ndot/2) 

89.51 

sec 

415.75 

meters 

(4  its) 

1752.9 

meters 

Space  Command’s  SGP4  (do  not  solve  for  drag  term  ndot/2) 

B*  =  0.000091  er'1 

84.32 

sec 

422.17 

meters 

(4  its) 

694.49 

meters 

The  results  in  Table  15  indicate  the  importance  of  modeling  solar  radiation  pressure  in  the  AOG  and  solving  for 
CR  (-200  meters  over  the  three  day  predict).  In  addition,  the  impact  of  higher  degree  and  order  zonal  and  tesseral  m- 
daily  terms,  as  well  as  shallow  tesseral  resonance  terms,  can  be  clearly  seen;  specifically,  neglecting  the  contributions 
from  these  effects  above  the  8x8  field  (up  to  the  21x21  gravity  model  limit  used  in  this  testing)  more  than  doubles  the 
predict  error.  Furthermore,  at  altitudes  characteristic  of  this  orbit,  the  effects  of  atmospheric  drag  and  third  body  short 
periodics  are  negligible,  while  SRP  short  periodics  are  on  the  order  of  5  meters.  Finally,  the  tesseral  high  frequency 
short  periodic,  J22  short  periodic,  and  J2  /  m-daily  short  periodic  terms  serve  to  provide  accuracy  closure  with  Cowell, 
but  are  computationally  expensive.  For  these  orbits,  their  effects  can  be  neglected. 

Again,  it  is  of  particular  importance  to  use  extreme  care  when  running  an  optimized  version  of  the  DSST.  This 
particular  orbit  does  experience  shallow  resonance  effects.  In  the  cases  where  only  an  8x8  gravity  field  (AOG)  has 
been  included,  the  resonant  effects  were  neglected.  In  order  to  capture  these  effects,  one  DSST  case  implemented  the 
“RES”  modeling  (use  of  the  RESNM  and  AUTOFORC  cards  with  a  small  base  field),  while  another  case 
implemented  a  21x21  geopotential  field  with  a  RESONPRD  card  set  equal  to  one  day  and  geopotential  short  periodic 
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terms  either  removed  (as  with  the  tesseral  high  frequency  short  periodics)  or  truncated  (as  with  the  zonal  and  m-daily 
short  periodics).  Both  of  these  configurations  capture  the  shallow  resonance  contributions;  however,  as  was  stated  in 
the  MACC  section,  the  configuration  using  a  21x21  geopotential  model  and  a  RESONPRD  card  set  to  one  day  is  (1) 
simpler  than  using  the  AUTOFORC  and  RESNM  configuration  and  (2)  applicable  for  all  cases  under  the  broad  scope 
of  the  HACC.  For  these  reasons,  the  21x21  geopotential  model  with  RESONPRD  equal  to  one  day  (and  geopotential 
short  periodics  either  removed  or  truncated)  was  chosen  for  the  optimal  configuration.  This  configuration  was  only 
5%  slower  than  the  more  complicated  AUTOFORC  and  RESNM  combination,  but  30  meters  more  accurate  in  the 
predict  span  (a  benefit  gained  by  the  AOG  zonal  terms  above  the  8x8  base  configuration  up  to  the  21x21  limit).  In  all, 
roughly  40  meters  is  gained  in  the  three  day  predict  with  the  addition  of  the  shallow  resonance  terms. 

The  setup  for  the  HACC  which  provides  an  optimal  balance  between  accuracy  and  computational  speed  can  now 
be  given: 

•  21x21  geopotential  (AOG) 

•  RESONPRD  equal  to  one  day  (86400.0  seconds) 

•  J22  (AOG) 

•  Truncated  zonal  short  periodic  modeling  (the  maximum  power  of  the  eccentricity  was  set  to  2  and  the 
maximum  frequency  in  true  longitude  was  set  toll  in  Fourier  coefficients  for  these  short  periodic  expansions) 

•  Truncated  tesseral  m-daily  short  periodic  modeling  (maximum  power  of  the  eccentricity  was  set  to  2  in  Fourier 
Coefficients  for  these  short  periodic  expansions) 

•  Lunar  /  solar  third  body  point  mass  effects  (AOG) 

•  Solar  Radiation  Pressure  (AOG) 

•  Solve  for  CR 

•  SSTESTFL  1  2  0  0.0 

•  SSTAPGFL  1  0  0  0.0  0.0  1.0 

Sample  input  decks  for  the  HACC  are  given  in  the  appendices. 

5.0  Geosynchronous  Case  (GEO) 

The  BS  3A  orbit  was  chosen  to  fulfill  the  geosynchronous  test  case.  In  the  ETAS  study,  the  GEO  is  characterized 
by  a  SMAH  range  of  33,000  -  39,000  km  and  an  eccentricity  range  from  0.0  -  0.05.  For  this  test  case,  multiple  runs 
were  made  in  which  varying  configurations  of  DSST  were  fit  to  a  Cowell  truth  model  over  a  three  day  arc.  A  three 
day  DSST  EG  with  the  DC  solve-for  parameters  followed  the  fit  span.  Accuracy  metrics  were  derived  from 
comparisons  between  the  DSST  and  Cowell  trajectories  over  both  the  fit  (DC  RMS)  and  predict  spans  (from  the  R&D 
GTDS  Ephemeris  Comparison  Program).  ORB  1  output  was  produced  every  450  seconds. 

The  osculating  elements  required  for  the  Cowell  truth  EG  were  derived  from  a  ten  day  Cowell  fit  to  an  SGP4 
ephemeris.  The  two  card  element  set  used  to  generate  the  SGP4  ephemeris  is  given  in  Table  16: 

Table  16.  BS  3A  Two-Card  Element  Set,  GEO,  September  1994  SATCAT 


Epoch  Date  YYMMDD  HHMMSS.SSSS 

940825  170334.9932 

Mean  Motion 

1.00266891  rev/day 

Eccentricity 

0.0003102 

Inclination 

0.0467  deg 

Longitude  of  Ascending  Node 

243.1923  deg 

Argument  of  Perigee 

288.7716  deg 
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Table  16.  BS  3A  Two-Card  Element  Set,  GEO,  September  1994  SATCAT  (Continued) 


Mean  Anomaly 

167.5966  deg 

B* 

0.0  er'1 

The  specifics  for  the  Cowell  DC  and  subsequent  EG  are  given  in  Table  17: 

Table  17.  Cowell  DC  and  Subsequent  EG  Inputs,  GEO 


Epoch  Date  YYMMDD  HHMMSS.SSSS 

940825  170334.9932 

Fit  Begin  Date  YYMMDD  HHMMSS.SSSS 

940825  170335.0000 

Fit  End  Date  YYMMDD  HHMMSS.SSSS 

940904  170335.0000 

EG  End  Date  YYMMDD  HHMMSS.SSSS 

940905  170335.0000 

Input,  Integration,  and  Output  Frame 

Earth  Equator  and 

Equinox  of  1950 
(Mean  of  1950) 

Integrator 

12th  Order  Summed  Cowell/ 
Adams  Predict-Partially  Correct 

Step  Size 

60.0  sec 

Gravity  Field 

21x21  JGM2 

Lunar  /  Solar  Point  Masses 

Yes 

Drag  Model 

No 

Solar  Radiation  Pressure 

Yes 

(Solve  Cr) 

Average  Spacecraft  Cross-Sectional  Area 

1.0  m2 

Spacecraft  Mass 

100  kg 

Position  Standard  Deviation 

500  meters 

Velocity  Standard  Deviation 

50  cm/sec 

The  same  initial  guess  of  the  state  vector  and  observation  noise  standard  deviation  were  used  for  all  of  the  DSST 
DC  runs.  These  values  are  highlighted  in  Table  18: 

Table  18.  DSST  DC  Specifics,  GEO 


Epoch  Date  YYMMDD  HHMMSS.SSSS 

940825  170334.9932 

Fit  Span  Begin  YYMMDD  HHMMSS.S 

940825  170335.0000 

Fit  Span  End  YYMMDD  HHMMSS.S 

940828  170335.0000 

Semimajor  Axis  (mean) 

42167.16 

Eccentricity  (mean) 

0.000275 

Inclination  (mean) 

0.23  deg 
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Table  18.  DSST  DC  Specifics,  GEO  (Continued) 


Longitude  of  Ascending  Node  (mean) 

97.6  deg 

Argument  of  Perigee  (mean) 

75.3  deg 

Mean  Anomaly  (mean) 

166.09  deg 

DSST  Input,  Integration  and  Output  Frame 

Earth  Equator  and 

Equinox  of  1950 
(Mean  of  1950) 

DSST  Step  Size 

43200.0  sec 

Geopotential  Model 

JGM2 

Position  Standard  Deviation 

100  meters 

Velocity  Standard  Deviation 

10  cm/sec 

In  the  DSST  DC  runs,  vaiying  configurations  corresponding  to  different  perturbation  modeling  were  used  to 
determine  the  optimal  balance  of  accuracy  and  computational  speed.  Results  of  the  DC  s  and  subsequent  EG  s  for  the 
GEO  are  given  in  Table  1 9;  first,  a  few  notes  concerning  the  test  protocol  for  the  GEO: 

•  Total  Position  RMS  is  derived  over  the  three  day  predict  span  only  (from  the  R&D  GTDS  Ephemeris 
Comparison  Program). 

•  R&D  GTDS  default  configuration  settings  are  used  unless  specifically  noted. 

•  “NUM  SP”  refers  to  numerical  short  periodics  used  in  conjunction  with  weak  time  dependence  (WTD)  and 

third  body  perturbations.  '  - 

•  All  cases  which  attempt  to  solve  for  CR  include  solar  radiation  pressure  effects  in  the  AOG  partials  (if  CR  is 
not  solved  for,  then  all  solar  radiation  pressure  partials  in  the  AOG  are  shut  off).  All  cases  include  J2  partials 
analytically  in  the  AOG.  No  other  AOG  or  SPG  partials  were  included  in  any  of  the  runs. 

•  All  cases  have  an  SSTESTFL  card  of  1  2;  this  setting  puts  only  J2  AOG  partials  into  the  A  matrix  (computed 
analytically).  All  other  desired  AOG  partials  for  non-conservative  perturbations  are  included  in  the  D  matrix 
via  finite  differences.  For  these  runs,  only  solar  radiation  pressure  effects  (if  solving  for  C^  have  been 
included  in  the  D  matrix. 

•  Cases  with  a  in  the  DC  RMS  column  had  a  DC  convergence  criteria  of  1  .D-3  (rather  than  1  .D-4  which  was 
used  for  all  other  cases). 
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Table  19.  DSST  DC  and  Subsequent  EG  Results,  GEO 


GRAV 

(AOG) 

Zonal 

SP 

M- 

daily 

SP 

TESS 

SP 

h2 

SP 

V 

m-daily 

SP 

Drag 

SRP 

3rd 

Body 

Run 

Time 

DC  RMS 

Total 

Position 

RMS 

(Predict 

Span) 

■ 

Yes 

Yes 

No 

Yes 

No 

No 

AOG 

SPG 

Solve 

CR 

L-S 

AOG 

SPG 

137.54 

sec 

8.24 

meters 

(29  its) 

228.17 

meters 

8x8 

h2 

Yes 

Yes 

No 

Yes 

No 

No 

AOG 

SPG 

Solve 

Cr 

L-S 

AOG 

SPG 

78.49 

sec 

8.24 

meters 

(29  its) 

228.17 

meters 

| 

Yes 

Yes 

No 

Yes 

No 

No 

AOG 

SPG 

Solve 

CR 

L-S 

AOG 

SPG 

72.01 

sec 

8.94 

meters 

(30  its) 

209.87 

meters 

Yes 

No 

No 

Yes 

No 

No 

AOG 

SPG 

Solve 

CR 

L-S 

AOG 

SPG 

69.32 

sec 

8.94 

meters 

(29  its) 

209.80 

meters 

4x4 

v 

(opt. 

deck) 

Yes 

No 

No 

Yes 

No 

No 

1 

AOG 

SPG 

Solve 

CR 

L-S 

AOG 

NUMSP 

WTD 

30.19 

sec 

5.10 

meters 

(4  its) 

21.249 

meters 

ffl 

Yes 

No 

No 

Yes 

No 

No 

AOG 

Solve 

CR 

L-S 

AOG 

NUMSP 

WTD 

3.54 

meters 

(12  its) 

25.578 

meters 

| 

Yes 

No 

No 

Yes 

No 

No 

No 

L-S 

AOG 

NUMSP 

WTD 

27.80 

sec 

80.99 

meters 

(3  its) 

* 

391.75 

meters 

| 

Yes 

No 

No 

Yes 

No 

No 

AOG 

SPG 

Solve 

CR 

No 

35.86 

sec 

714.81 

meters 

(9  its) 

4046.5 

meters 
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Table  19.  DSST  DC  and  Subsequent  EG  Results,  GEO  (Continued) 


GRAV 

(AOG) 

Zonal 

SP 

M- 

daily 

SP 

TESS 

SP 

J22 

SP 

V 

m-daily 

SP 

Drag 

SRP 

3rd 

Body 

Run 

Time 

DC  RMS 

Total 

Position 

RMS 

(Predict 

Span) 

4x0 

h2 

Yes 

No 

No 

Yes 

No 

No 

AOG 

SPG 

Solve 

Cr 

L-S 

AOG 

NUM  SP 

WTD 

46.35 

sec 

24.74 

meters 

(16  its) 

7229.1 

meters 

Spa 

ce  Command’s  SGP4  (solve  for  drag  term  ndot/2) 

59.67 

sec 

159.50 

meters 

(30  its) 

* 

3934.5 

meters 

Space  Command’s  SGP4  (do  not  solve  for  drag  term  ndot/2) 

B*  =  0.0  er  1 

36.21 

sec 

479.85 

meters 

(12  its) 

* 

2060.2 

meters 

The  results  in  Table  19  indicate  a  small  degree  and  order  geopotential  model  (i.e.,  4x4)  can  be  used  without  a 
significant  loss  in  accuracy.  In  addition,  the  importance  of  third  body  effects  and  solar  radiation  pressure  (solving  for 
Cg)  can  clearly  be  seen.  Completely  neglecting  third  body  effects  leads  to  ~  4  km  worth  of  predict  error  over  the  three 
day  span',  neglecting  solar  radiation  pressure  in  both  the  AOG  and  SPG  leads  to  roughly  a  350-400  meter  error  over 
the  three  day  predict.  Neglecting  SRP  short  periodics  (but  still  modeling  SRP  in  the  AOG  and  solving  for  Cr),  does 
not  significantly  impact  accuracy  (~  5  m),  but  causes  extra  DC  iterations  for  convergence  (and,  hence,  longer  run 
times).  Furthermore,  weak  time  dependence,  which  accounts  for  the  motion  of  the  third  bodies  over  the  course  of  the 
satellite’s  orbit  in  the  equations  of  motion,  must  be  modeled.  A  weak  time  independent  theory  assumes  the  position  of 
the  third  body  does  not  move  over  the  course  of  an  averaging  interval  (the  orbital  period)  in  the  development  of  the 
equations  of  motion.  For  low  altitude  orbits,  this  assumption  works  well.  However,  for  geosynchronous  orbits  (with 
periods  of  roughly  24  hours),  these  results  indicate  using  a  weak  time  independent  theory  leads  to  extra  DC  iterations 
for  convergence  and  roughly  a  200  meter  predict  error  over  the  course  of  three  days. 

The  case  using  a  4x0  geopotential  field  illustrates  the  impact  of  resonance  on  the  geosynchronous  orbit.  For  the 
three  day  predict,  over  7  km  worth  of  error  is  introduced  by  neglecting  resonant  terms.  The  impact  of  using  a  subset 
of  resonant  terms  (rather  the  the  full  complement  for  the  1:1  resonance)  can  clearly  be  seen  in  the  SGP4  results.  The 
results  for  optimized  DSST  indicate  the  efficiency  with  which  a  lull  complement  of  resonant  terms  can  be  added  to  a 
propagation  theory. 

The  setup  for  the  GEO  which  provides  an  optimal  balance  between  accuracy  and  computational  speed  can  now 
be  given: 

•  4x4  geopotential  (AOG) 

•  J22  (AOG  and  SPG) 

•  Zonal  short  periodics  with  default  settings 

•  Lunar  /  solar  third  body  point  mass  effects  (AOG) 

•  Numerical  third  body  short  periodics  (WTD) 

•  Solar  radiation  pressure  (AOG  &  SPG) 


32 


•  Solve  for  CR 

•  SSTESTFL  120  0.0 

•  SSTAPGFL  100  0.0  0.01.0 

Sample  input  decks  for  the  GEO  are  given  in  the  appendices. 

6.0  Molniya  Case  (MOLY) 

Satellite  #9829  (SATCAT)  was  chosen  to  fulfill  the  Molniya  test  case14,20.  In  the  ETAS  study,  the  Molniya  / 
geostationary  transfer  orbit  is  characterized  by  a  SMAH  range  from  18,000  -  22,000  km  and  an  eccentricity  range 
from  0.05  to  1 .0  (this  analysis  will  strictly  focus  on  the  Molniya  orbit  which,  typically,  has  an  eccentricity  of  -  0.73). 

For  this  test  case,  multiple  runs  were  made  in  which  varying  configurations  of  DSST  were  fit  to  a  Cowell  truth 
model  over  a  three  day  arc.  A  three  day  DSST  EG  with  the  DC  solve-for  parameters  followed  the  fit  span.  Accuracy 
metrics  were  derived  from  comparisons  between  the  DSST  and  Cowell  trajectories  over  both  the  fit  (DC  RMS)  and 
predict  spans  (from  the  R&D  GTDS  Ephemeris  Comparison  Program).  ORB1  output  was  produced  every  450 
seconds. 

The  osculating  elements  required  for  the  Cowell  truth  EG  were  derived  from  a  six  month  DSST  fit  to  real  data14. 
The  six  month  DSST  fit  to  the  real  data  resulted  in  a  set  of  mean  Keplerian  elements  which  were  propagated  forward 
one  day  to  produce  the  osculating  elements  used  in  the  Cowell  truth  EG.  The  mean  Keplerian  elements,  as  well  as 
other  settings  used  in  the  one  day  DSST  propagation,  are  listed  in  Table  20: 

Table  20.  Mean  Element  Set  Used  To  Build  Osculating  Elements,  MOLY  (from  Reference  14) 


Epoch  Date  YYMMDD  HHMMSS.SSSS 

790803  234212.0000 

Semi-Major  Axis  (mean) 

26556.9582  km 

Eccentricity  (mean) 

0.6990986 

Inclination  (mean) 

63.173001  deg 

Longitude  of  Ascending  Node  (mean) 

190.619681  deg 

Argument  of  Perigee  (mean) 

281.59624  deg 

Mean  Anomaly  (mean) 

13.29315  deg 

Input,  Integration,  and  Output  Frame 

Earth  Equator  and 

Equinox  of  1950 
(Mean  of  1950) 

Integrator 

4th  Order  Runge-Kutta 

DSST  Step  Size 

43200.0  sec 

Gravity  Field 

8x8  GEM9 

Lunar  /  Solar  Point  Masses 

Yes 

Drag  Model 

Jacchia-Roberts 

Cd  =  2.0 

Solar  Radiation  Pressure 

Yes 

Cr-1.2 
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Table  20.  Mean  Element  Set  Used  To  Build  Osculating  Elements,  MOLY  (from  Reference  14) 


Global  Short  Periodic  Select 

Moderate  Accuracy  Molniya 
Configuration  (SPSHPER  =  5) 

Average  Spacecraft  Cross-Sectional  Area 

12.5  m2 

Spacecraft  Mass 

1250  kg 

The  specifics  for  the  Cowell  truth  EG  are  given  in  Table  21 : 

Table  21.  Cowell  Truth  EG  Inputs,  MOLY 


Epoch  Date  YYMMDD  HHMMSS.SSSS 

790804  234212.0000 

EG  End  Date  YYMMDD  HHMMSS.SSSS 

790818  234212.0000 

Input,  Integration,  and  Output  Frame 

Earth  Equator  and 

Equinox  of  1950 
(Mean  of  1950) 

Integrator 

Time  Regularized  Cowell 

Time  Regularization  Constant  for  Integrator 

1.5 

Step  Size 

200  Step  Per  Revolution 

Gravity  Field 

21x21  GEM  1 0B 

Lunar  /  Solar  Point  Masses 

Yes 

Drag  Model 

Jacchia-Roberts 

CD  =  2.0 

Solar  Radiation  Pressure 

Yes 

Cr=1-2 

Average  Spacecraft  Cross-Sectional  Area 

12.5  m2 

Spacecraft  Mass 

1250  kg 

The  same  initial  guess  of  the  state  vector  and  observation  noise  standard  deviation  were  used  for  all  of  the  DSST 
DC  runs.  These  values  are  highlighted  in  Table  22: 

Table  22.  DSST  DC  Specifics,  MOLY 


Epoch  Date  YYMMDD  HHMMSS.SSSS 

790804  234212.0000 

Fit  Span  Begin  YYMMDD  HHMMSS.S 

790804  234212.0000 

Fit  Span  End  YYMMDD  HHMMSS.S 

790807  234212.0000 

Semimajor  Axis  (mean) 

26572.176  km 

Eccentricity  (mean) 

0.699 

Inclination  (mean) 

63.2  deg 

Longitude  of  Ascending  Node  (mean) 

190.5  deg 
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Table  22.  DSST  DC  Specifics,  MOLY  (Continued) 


Argument  of  Perigee  (mean) 

281.6  deg 

Mean  Anomaly  (mean) 

15.429  deg 

DSST  Input,  Integration  and  Output  Frame 

Earth  Equator  and 

Equinox  of  1950 
(Mean  of  1950) 

Geopotential  Model 

GEM10B 

Position  Standard  Deviation 

100  meters 

Velocity  Standard  Deviation 

10  cm/sec 

In  the  DSST  DC  runs,  varying  configurations  corresponding  to  different  perturbation  modeling  were  used  to 
determine  the  optimal  balance  of  accuracy  and  computational  speed.  Results  of  the  DC’s  and  subsequent  EG’s  for  the 
MOLY  are  given  in  Table  23;  again,  it  is  first  desirable  to  highlight  a  few  notes  concerning  the  test  protocol  for  the 
MOLY: 

•  Total  Position  RMS  is  derived  over  the  three  day  predict  span  only  (from  R&D  GTDS  Ephemeris 
Comparison  Program). 

•  R&D  GTDS  default  configuration  settings  are  used  unless  specifically  noted. 

•  “trunc”  refers  to  truncated  modeling;  modeling  descriptions  corresponding  to  the  R&D  GTDS  Semianalytic 
Satellite  Theory  input  processor  are  specified  where  needed.  Additional  details  can  be  found  in  the  R&D 
GTDS  Semianalytic  Theory  Input  Processor  keyword  description  booklet33. 

•  “splunara”  and  “spsolara”  also  refer  to  truncated  modeling;  additional  details  can  be  found  in  the  R&D  GTDS 
Semianalytic  Theoiy  Input  Processor  keyword  description  booklet33. 

•  All  cases  which  attempt  to  solve  for  CD  include  first  order  drag  effects  in  the  AOG  partials  (if  CD  is  not  solved 
for,  then  all  drag  partials  in  the  AOG  are  shut  off).  Similarly,  all  cases  which  attempt  to  solve  for  CR  include 
solar  radiation  pressure  effects  in  the  AOG  partials  (if  CR  is  not  solved  for,  then  all  solar  radiation  pressure 
partials  in  the  AOG  are  shut  off).  All  cases  include  J2  partials  analytically  in  the  AOG.  No  other  AOG  or  SPG 
partials  were  included  in  any  of  the  runs. 

•  Most  cases  have  an  SSTESTFL  card  of  1  3;  this  setting  includes  all  desired  AOG  partials  in  the  A  matrix  (J2 
partials  are  computed  analytically;  all  other  partial  derivatives  are  computed  via  finite  differences).  In  addition, 
this  setting  includes  desired  AOG  partial  derivatives  for  non-conservative  perturbations  in  the  D  matrix  via 
finite  differences  (these  partials  were  also  included  in  the  A  matrix).  Cases  noted  by  **  have  an  SSTESTFL 
card  of  1  2;  this  setting  puts  only  J2  AOG  partials  into  the  A  matrix  (computed  analytically).  All  other  desired 
AOG  partials  for  non-conservative  perturbations  are  included  in  the  D  matrix  via  finite  differences  (with  this 
setting,  third  body  and  J22  AOG  partial  derivatives,  if  desired,  would  not  be  included  anywhere;  if  they  are 
desired,  the  SSTESTFL  setting  of  1  3  should  be  used). 
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Table  23.  DSST  DC  and  Subsequent  EG  Results,  MOLY 


GRAV 

(AOG) 

Zonal 

SP 

M- 

daily 

SP 

TESS 

SP 

J22 

SP 

V 

m-daily 

SP 

Drag 

SRP 

3rd 

Body 

Run 

Time 

DC  RMS 

Total 

Position 

RMS 

(Predict 

Span) 

8x8 

J22 

Yes 

Yes 

tranc 

88 

33  35 

-4141 

Yes 

No 

J-R 

AOG 

SPG 

Solve 

CD 

AOG 

SPG 

Solve 

CR 

L-S 

AOG 

SPG 

174.17 

sec 

20.62 

meters 

(5  its) 

149.19 

meters 

8x8 

T  2 

J2 

Yes 

Yes 

trunc 

44 

33  35 

-1515 

Yes 

No 

J-R 

AOG 

SPG 

Solve 

CD 

AOG 

SPG 

Solve 

CR 

L-S 

AOG 

SPG 

60.24 

sec 

33.15 

meters 

(5  its) 

146.27 

meters 

8x8 

J22 

Yes 

Yes 

No 

Yes 

No 

J-R 

AOG 

SPG 

Solve 

CD 

AOG 

SPG 

Solve 

CR 

L-S 

AOG 

SPG 

43.89 

sec 

117.30 

meters 

(5  its) 

220.10 

meters 

00  fS 
£  * 

Yes 

Yes 

No 

Yes 

! 

No 

J-R 

AOG  ^ 
SPG 

Solve 

CD 

AOG 

Solve 

CR 

L-S 

AOG 

SPG 

39.92 

sec 

117.65 

meters 

(4  its) 

220.08 

meters 

8x8 

J22 

Yes 

Yes 

No 

Yes 

No 

J-R 

AOG 

Solve 

CD 

AOG 

Solve 

CR 

L-S 

AOG 

SPG 

39.00 

sec 

117.65 

meters 

(4  its) 

220.08 

meters 

8x8 

T  2 

J2 

Yes 

Yes 

No 

Yes 

No 

J-R 

AOG 

Solve 

CD 

AOG 

Solve 

CR 

L-S 

AOG 

38.49 

sec 

208.28 

meters 

(4  its) 

357.40 

meters 

8x8 

J22 

Yes 

Yes 

No 

Yes 

No 

J-R 

AOG 

Solve 

CD 

AOG 

Solve 

CR 

No 

42.90 

sec 

1027.38 

meters 

(6  its) 

12670.0 

meters 

Table  23.  DSST  DC  and  Subsequent  EG  Results,  MOLY  (Continued) 


GRAV 

(AOG) 

Zonal 

SP 

M- 

daily 

SP 

TESS 

SP 

T  ^ 

J2 

SP 

J 2I 

m -daily 

SP 

Drag 

SRP 

3rd 

Body 

Run 

Time 

DC  RMS 

Total 

Position 

RMS 

(Predict 

Span) 

8x8 

h2 

Yes 

Yes 

No 

Yes 

No 

J-R 

AOG 

Solve 

CD 

AOG 

Solve 

CR 

L-S 

AOG 

splunara 

4482 

spsolara 

2242 

38.83 

sec 

116.66 

meters 

(4  its) 

219.04 

meters 

8x8 

ti 

Yes 

Yes 

No 

Yes 

No 

J-R 

AOG 

Solve 

CD 

No 

L-S 

AOG 

SPG 

38.62 

sec 

121.09 

meters 

(4  its) 

256.95 

meters 

8x8 

h2 

Yes 

Yes 

No 

Yes 

No 

No 

AOG 

Solve 

CR 

L-S 

AOG 

SPG 

34.77 

sec 

117.63 

meters 

(4  its) 

218.75 

meters 

8x8 

h2 

trunc 

6  510 

^^3 

No 

Yes 

No 

J-R 

AOG 

Solve 

CD 

AOG 

Solve 

CR 

L-S 

AOG 

SPG 

114.62 

meters 

(4  its) 

219.37 

meters 

8x8 

h2 

trunc 

43  8 

trunc 

442 

No 

Yes 

No 

J-R 

AOG 

Solve 

CD 

AOG 

Solve 

CR 

L-S 

AOG 

SPG 

36.59 

sec 

114.82 

meters 

(4  its) 

219.33 

meters 

8x8 

h2 
*  * 

(opt. 

deck) 

trunc 

438 

No 

Yes 

No 

J-R 

AOG 

Solve 

Cd 

AOG 

Solve 

CR 

L-S 

AOG 

SPG 

34.11 

sec 

114.82 

meters 

(4  its) 

219.33 

meters 

Space  Command’s  SGP4  (solve  for  drag  term  ndot/2) 

35.41 

sec 

2690.67 

meters 

(8  its) 

5519.0 

meters 

Space  Command’s  SGP4  (do  not  solve  for  drag  term  ndot/2) 

B*  =  0.0001  er 1 

28.29 

sec 

3082.99 

meters 

(4  its) 

4537.3 

meters 

The  results  in  Table  23  demonstrate  the  wide  variety  of  perturbative  effects  experienced  by  Molniya  orbits.  The 
low  perigee  heights  contribute  to  geopotential  and  atmospheric  effects,  while  the  high  apogee  heights  lead  to  third 
body  and  solar  radiation  pressure  effects.  Of  particular  importance  are  the  effects  of  tesseral  resonance  (Molniya 
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orbits  usually  maintain  a  repeating  groundtrack  and  complete  very  nearly  two  revolutions  per  day).  If  neglected, 
resonance  can  contribute  several  kilometers  worth  of  error  (as  is  evident  in  the  analytic  DCs,  which  only  model  a 
subset  of  the  resonant  terms)20.  Specifically,  the  following  notes  can  be  taken  from  Table  23  concerning  the  short 

periodic  models  of  DSST: 

•  drag  and  solar  radiation  pressure  short  periodics  do  not  significantly  affect  accuracy;  they  can  be  truncated  to 
reduce  computation  time 

•  neglecting  third  body  short  periodics  adds  ~  100-150  meters  worth  of  error  over  the  three  day  predict 

•  zonal  and  m-daily  short  periodics  can  be  truncated  to  decrease  computation  time  without  significantly 
impacting  accuracy 

•  tesseral  high  frequency  short  periodics  serve  to  provide  accuracy  closure  with  Cowell  (they  have  an  impact  on 
the  order  of  100  meters),  but  are  computationally  expensive 

These  results  also  indicate  neglecting  third  body  effects  (both  AOG  and  SPG)  leads  to  almost  13  km  worth  of 
error  over  the  three  day  predict. 

In  the  first  high  accuracy  DSST  run,  the  DC  RMS  error  is  on  the  order  of  20  meters.  This  error  can  be  attributed 
to  geopotential  effects  beyond  the  8x8  configuration  for  DSST  (up  to  the  21x21  limit  in  the  Cowell  truth  model) 
truncations  to  the  equations  of  motion  in  DSST’s  J22  model  (Ofe1]  in  the  AOG  and  0[e  ]  m  the  SPG),  neglecting  V 
m-daily  coupling,  an  incomplete  tesseral  high  frequency  short  periodic  model  at  the  8x8  limit  (SPTESSLC  8  8  33 
35-4141),  and  an  error  due  to  a  discrepancy  between  the  truth  and  DC  “solved-for”  values  of  CD  and  CR.  If  this  first 
high  accuracy  DSST  DC  is  repeated  with  the  exact  values  of  CD  and  CR  used  in  the  truth  (CD  =  2.0,  CR  =  1.2),  a  DC 
RMS  of  ~  13  meters  results.  This  implies  roughly  7  to  8  meters  of  the  original  ~  20  meter  DC  RMS  error  is  due  to 
having  poor  DC  solutions  for  CD  and  CR.  Further  analysis  implementing  an  8x8  DSST  fit  (SPTESSLC  =  8  8  33  35  - 
4141  J2  /  m-daily  off)  to  an  8x8  Cowell  truth  model  in  which  both  the  fit  and  truth  theories  used  a  small  value  of  J2 
(1.0  T J-6)  and  geopotential  effects  only  results  in  a  DC  RMS  of-  3  meters.  This  3  meter  error  can  be  attributed  to  the 
incomplete  tesseral  high  frequency  short  periqdic  model  of  DSST  for  the  8x8  configuration  (remember,  the  8x8 
DSST  configuration  isn’t  really  “full-up”;  rather,  an  SPTESSLC  of  8  8  33  35  -41  41  was  used).  If  this  same  test  is 
repeated  except  for  using  the  regular  value  of  J2,  an  error  of  ~  7  meters  results.  This  means  the  truncations  in  DSST 
J22  model  and  the  neglect  of  J2  /  m-daily  coupling  terms  results  in  -  4  meters  of  error.  This  leaves  an  approximate  5  to 
6  meter  error  due  to  geopotential  terms  beyond  the  8x8  configuration  up  to  the  21x21  limit  used  in  the  truth.  As  a 
final  note,  it  should  be  mentioned  this  specific  breakdown  of  errors  is  dependent  upon  the  geometry  of  this  particular 
oibit,  the  relationship  between  the  fit  span  and  perigee  locations,  etc.  Therefore,  these  error  quantifications  should  be 
used  as  a  “rule  of  thumb”  metrics. 

The  setup  for  the  MOLY  which  provides  an  optimal  balance  between  accuracy  and  computational  speed  can  now 
be  given: 

•  8x8  geopotential  (AOG) 

•  J22  (AOG  and  SPG) 

•  Truncated  central  body  zonal  harmonic  short  periodics  (SPZONALS  =  438) 

•  Truncated  m-daily  tesseral  harmonic  short  periodics  (SPMDAILY  =  442) 

•  Lunar  /  solar  third  body  point  mass  effects  (AOG  &  SPG) 

•  Atmospheric  drag  (AOG) 

•  Solve  for  CD 

•  Solar  radiation  pressure  (AOG) 

•  Solve  for  CR 

•  SSTESTFL  1  2  0  0.0 

•  SSTAPGFL  1  0  0  1.0  0.0  1.0 

Sample  input  decks  for  the  MOLY  are  given  in  the  appendices. 
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7.0  Near  Earth  Eccentric  Case  (NEEC) 


The  Vanguard  2  orbit  was  chosen  to  fulfill  the  near  earth  eccentric  test  case.  In  the  ETAS  study,  the  NEEC  is 
characterized  by  a  SMAH  range  from  0  -  2500  km  and  an  eccentricity  range  from  0.05  to  1.0.  For  this  test  case, 
multiple  runs  were  made  in  which  varying  configurations  of  DSST  were  fit  to  a  Cowell  truth  model  over  a  three  day 
arc.  A  three  day  DSST  EG  with  the  DC  solve-for  parameters  followed  the  fit  span.  Accuracy  metrics  were  derived 
from  comparisons  between  the  DSST  and  Cowell  trajectories  over  both  the  fit  (DC  RMS)  and  predict  spans  (from  the 
R&D  GTDS  Ephemeris  Comparison  Program).  ORB1  output  was  produced  every  450  seconds. 

The  osculating  elements  required  for  the  Cowell  truth  EG  were  derived  from  a  five  day  Cowell  fit  to  an  SGP4 
ephemeris.  The  two  card  element  set  used  to  generate  the  SGP4  ephemeris  is  given  in  Table  24: 

Table  24.  Vanguard  2  Two-Card  Element  Set,  NEEC,  September  1994  SATCAT 


Epoch  Date  YYMMDD  HHMMSS.SSSS 

940826  073513.6735 

Mean  Motion 

11.73921485  rev/day 

Eccentricity 

0.1522640 

Inclination 

32.8834  deg 

Longitude  of  Ascending  Node 

251.8592  deg 

Argument  of  Perigee 

10.8368  deg 

Mean  Anomaly 

352.1515  deg 

B* 

0.000154  er1 

The  specifics  for  the  Cowell  DC  and  subsequent  EG  are  given  in  Table  25: 

Table  25.  Cowell  DC  and  Subsequent  EG  Inputs,  NEEC 


Epoch  Date  YYMMDD  HHMMSS.SSSS 

940826  073513.6735 

Fit  Begin  Date  YYMMDD  HHMMSS.SSSS 

940826  073514.0000 

Fit  End  Date  YYMMDD  HHMMSS.SSSS 

940831  073514.0000 

EG  End  Date  YYMMDD  HHMMSS.SSSS 

940905  073514.0000 

Input,  Integration,  and  Output  Frame 

Earth  Equator  and 

Equinox  of  1950 
(Mean  of  1950) 

Integrator 

12th  Order  Summed  Cowell/ 
Adams  Predict-Partially  Correct 

Step  Size 

60.0  sec 

Gravity  Field 

21x21  JGM2 

Lunar  /  Solar  Point  Masses 

Yes 

Drag  Model 

Jacchia-Roberts-Schatten 

(Solve  px) 
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Table  25.  Cowell  DC  and  Subsequent  EG  Inputs,  NEEC  (Continued) 


Solar  Radiation  Pressure 

Yes 

(Solve  Cr) 

Average  Spacecraft  Cross-Sectional  Area 

1.0  m2 

Spacecraft  Mass 

100  kg 

Position  Standard  Deviation 

500  meters 

Velocity  Standard  Deviation 

50  cm/sec 

The  same  initial  guess  of  the  state  vector  and  observation  noise  standard  deviation  were  used  for  all  of  the  DSST 
DC  runs.  These  values  are  highlighted  in  Table  26: 

Table  26.  DSST  DC  Specifics,  NEEC 


Epoch  Date  YYMMDD  HHMMSS.SSSS 

940826  073513.6735 

Fit  Span  Begin  YYMMDD  HHMMSS.S 

940826  073514.0 

Fit  Span  End  YYMMDD  HHMMSS.S 

940829  073514.0 

Semimajor  Axis  (mean) 

8177.913  km 

Eccentricity  (mean) 

0.1522640 

Inclination  (mean) 

32.8834  deg 

Longitude  of  Ascending  Node  (mean) 

251.8592  deg 

Argument  of  Perigee  (mean) 

10.8368  deg 

Mean  Anomaly  (mean) 

352.1515  deg 

DSST  Input  Frame 

True  of  Date 

DSST  Integration  and  Output  Frame 

Earth  Equator  and 

Equinox  of  1950 
(Mean  of  1950) 

Geopotential  Model 

JGM2 

Position  Standard  Deviation 

100  meters 

Velocity  Standard  Deviation 

10  cm/sec 

In  the  DSST  DC  runs,  varying  configurations  corresponding  to  different  perturbation  modeling  were  used  to 
determine  the  optimal  balance  of  accuracy  and  computational  speed.  Results  of  the  DC’s  and  subsequent  EG’s  for  the 
NEEC  are  given  in  Table  27;  however,  it  is  first  desirable  to  highlight  a  few  notes  concerning  the  test  protocol  for  the 
NEEC: 


•  Total  Position  RMS  is  derived  over  the  three  day  predict  span  only  (from  R&D  GTDS  Ephemeris 
Comparison  Program). 

•  R&D  GTDS  default  configuration  settings  are  used  unless  specifically  noted. 

•  “trunc”  refers  to  truncated  modeling;  modeling  descriptions  corresponding  to  the  R&D  GTDS  Semianalytic 
Satellite  Theory  input  processor  are  specified  where  needed.  Additional  details  can  be  found  in  the  R&D 
GTDS  Semianalytic  Theory  Input  Processor  keyword  description  booklet33. 
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•  RES  refers  to  the  inclusion  of  the  11th  and  12th  order  shallow  resonance  terms  with  the  AUTOFORC  and 
RESNM  keyword  cards. 

•  All  cases  which  attempt  to  solve  for  CD  include  first  order  drag  effects  in  the  AOG  partials  (if  CD  is  not  solved 
for,  then  all  drag  partials  in  the  AOG  are  shut  off).  Similarly,  all  cases  which  attempt  to  solve  for  Cr  include 
solar  radiation  pressure  effects  in  the  AOG  partials  (if  Cr  is  not  solved  for,  then  all  solar  radiation  pressure 
partials  in  the  AOG  are  shut  off).  All  cases  include  J2  partials  analytically  in  the  AOG.  No  other  AOG  or  SPG 
partials  were  included  in  any  of  the  runs. 

•  Most  cases  have  an  SSTESTFL  card  of  1  3;  this  setting  includes  all  desired  AOG  partials  in  the  A  matrix  (J2 
partials  are  computed  analytically;  all  other  partial  derivatives  are  computed  via  finite  differences).  In  addition, 
this  setting  includes  desired  AOG  partial  derivatives  for  non-conservative  perturbations  in  the  D  matrix  via 
finite  differences  (these  partials  were  also  included  in  the  A  matrix).  Cases  noted  by  **  have  an  SSTESTFL 
card  of  1  2;  this  setting  puts  only  J2  AOG  partials  into  the  A  matrix  (computed  analytically).  All  other  desired 
AOG  partials  for  non-conservative  perturbations  are  included  in  the  D  matrix  via  finite  differences  (with  this 
setting,  third  body  and  J22  AOG  partial  derivatives,  if  desired,  would  not  be  included  anywhere;  if  they  are 
desired,  the  SSTESTFL  setting  of  1  3  should  be  used). 


Table  27.  DSST  DC  and  Subsequent  EG  Results,  NEEC 


GRAV 

(AOG) 

Zonal 

SP 

M- 

daily 

SP 

TESS 

SP 

ti 

SP 

V 

m-dai!y 

SP 

Drag 

SRP 

3rd 

Body 

Run 

Time 

DC  RMS 

Total 

Position 

RMS 

(Predict 

Span) 

21x21 

Yes 

Yes 

trunc 

Yes 

Yes 

J-R 

AOG 

L-S 

541.29 

41.76 

223.16 

J22 

2121 

AOG 

SPG 

AOG 

sec 

meters 

meters 

RESON 

PRD  = 

1  day 

21  15 

-37  37 

SPG 

Solve 

CD 

Solve 

Cr 

SPG 

(6  its) 

21x21 

Yes 

Yes 

trunc 

Yes 

Yes 

J-R 

AOG 

L-S 

202.46 

42.76 

234.47 

h1 

6  6  21 

AOG 

SPG 

AOG 

sec 

meters 

meters 

RESON 

PRD  = 

1  day 

15 

-27  27 

SPG 

Solve 

CD 

Solve 

CR 

SPG 

! 

(9  its) 

21x21 

Yes 

Yes 

trunc 

Yes 

No 

J-R 

AOG 

L-S 

132.27 

42.83 

232.69 

h1 

6621 

AOG 

SPG 

AOG 

sec 

meters 

meters 

RESON 

PRD  = 

1  day 

15 

-27  27 

SPG 

Solve 

CD 

Solve 

Cr 

SPG 

(8  its) 

21x21 

Yes 

Yes 

No 

Yes 

No 

J-R 

AOG 

L-S 

70.74 

55.08 

239.42 

ti 

AOG 

SPG 

AOG 

sec 

meters 

meters 

RESON 

SPG 

Solve 

SPG 

(9  its) 

PRD  = 

1  day 

Solve 

CD 

CR 
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Table  27.  DSST  DC  and  Subsequent  EG  Results,  NEEC  (Continued) 


GRAV 

(AOG) 

Zonal 

SP 

M- 

daily 

SP 

TESS 

SP 

J22 

SP 

V 

m -daily 

SP 

Drag 

SRP 

3rd 

Body 

Run 

Time 

DC  RMS 

Total 

Position 

RMS 

(Predict 

Span) 

12x12 

Yes 

Yes 

No 

Yes 

No 

J-R 

AOG 

L-S 

57.42 

55.37 

456.24 

h2 

AOG 

SPG 

AOG 

sec 

meters 

meters 

RESON 

SPG 

Solve 

SPG 

(9  its) 

PRD  = 

1  day 

Solve 

CD 

CR 

8x8 

Yes 

Yes 

No 

Yes 

No 

J-R 

AOG 

L-S 

48.50 

58.63 

389.63 

J22 

AOG 

SPG 

AOG 

sec 

meters 

meters 

SPG 

Solve 

SPG 

(8  its) 

Solve 

CD 

CR 

8x8 

Yes 

Yes 

No 

Yes 

No 

J-R 

AOG 

L-S 

47.46 

61.08 

207.29 

h2 

AOG 

SPG 

AOG 

sec 

meters 

RES 

SPG 

Solve 

SPG 

(6  its) 

Solve 

CD 

CR 

Trunc. 

Trunc. 

No 

Yes 

No 

J-R 

AOG 

L-S 

47.67 

61.06 

229.95 

h2 

87  17 

886 

AOG 

SPG 

AOG 

sec 

meters 

meters 

RESON 

SPG 

Solve 

SPG 

(6  its) 

PRD  = 

1  day 

Solve 

CD 

CR 

21x21 

Trunc. 

No 

Yes 

No 

J-R 

AOG 

L-S 

45.92 

61.06 

229.95 

h2 

886 

AOG 

SPG 

AOG 

sec 

meters 

meters 

RESON 

PRD  = 

1  day 

1 

Solve 

CD 

Solve 

CR 

SPG 

(6  its) 

Trunc. 

No 

Yes 

No 

J-R 

AOG 

L-S 

45.54 

61.06 

229.94 

j22 

886 

AOG 

Solve 

AOG 

sec 

meters 

meters 

RESON 

■ 

Solve 

CR 

SPG 

(6  its) 

PRD  — 

1  day 

CD 

21x21 

Trunc. 

No 

Yes 

No 

J-R 

AOG 

L-S 

44.10 

60.94 

228.52 

T  ^ 

J2 

886 

AOG 

Solve 

AOG 

sec 

meters 

meters 

RESON 

PRD  “ 

1  day 

1 

Solve 

CD 

CR 

(6  its) 
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Table  27.  DSST  DC  and  Subsequent  EG  Results,  NEEC  (Continued) 


GRAV 

(AOG) 

Zonal 

SP 

M- 

daily 

SP 

TESS 

SP 

h2 

SP 

V 

m-daily 

SP 

Drag 

SRP 

3rd 

Body 

Run 

Time 

DC  RMS 

Total 

Position 

RMS 

(Predict 

Span) 

21x21 

h1 

RESON 

PRD  = 

1  day 

Trunc. 

664 

No 

Yes 

No 

J-R 

AOG 

Solve 

CD 

AOG 

Solve 

CR 

L-S 

AOG 

42.63 

sec 

69.54 

meters 

(6  its) 

217.87 

meters 

21x21 

h2 

RESON 

PRD  ~ 

1  day 

| 

No 

Yes 

No 

J-R 

AOG 

Solve 

cD 

AOG 

Solve 

CR 

L-S 

AOG 

46.30 

sec 

78.95 

meters 

(8  its) 

238.46 

meters 

21x21 

h2 

it  it 

RESON 

PRD  = 

1  day 

(opt. 

deck) 

Tninc. 

6513 

Trunc. 

664 

No 

Yes 

No 

J-R 

AOG 

Solve 

CD 

AOG 

Solve 

CR 

L-S 

AOG 

37.89 

sec 

69.42 

meters 

(6  its) 

218.71 

meters 

21x21 

T  2 

J2 

it  it 

RESON 

PRD  = 

1  day 

^^9 

Trunc. 

664 

No 

Yes 

No 

No 

AOG 

Solve 

CR 

L-S 

AOG 

i 

i 

! 

35.97 

sec 

78.48 

meters 

(6  its) 

312.51 

meters 

21x21 

J22 

it  it 

RESON 

PRD  = 

1  day 

a 

Trunc. 

664 

No 

Yes 

No 

J-R 

AOG 

Solve 

CD 

No 

L-S 

AOG 

36.75 

sec 

67.82 

meters 

(6  its) 

267.96 

meters 

21x21 

h2 

it  it 

RESON 

PRD  = 

1  day 

Trunc. 

6513 

Trunc. 

664 

No 

Yes 

No 

J-R 

AOG 

Solve 

CD 

AOG 

Solve 

CR 

No 

i 

38.06 

sec 

72.15 

meters 

(6  its) 

266.27 

meters 
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Table  27.  DSST  DC  and  Subsequent  EG  Results,  NEEC  (Continued) 


GRAV 

(AOG) 

Zonal 

SP 

M- 

daily 

SP 

TESS 

SP 

h2 

SP 

V 

m-daily 

SP 

Drag 

SRP 

3rd 

Body 

Run 

Time 

DC  RMS 

Total 

Position 

RMS 

(Predict 

Span) 

Space  Command’s  SGP4  (solve  for  drag  term  ndot/2) 

28.27 

sec 

909.86 

meters 

(5  its) 

1080.7 

meters 

Space  Command’s  SGP4  (do  not  solve  for  drag  term  ndot/2) 

B*  =  0.000154  er'1 

26.69  | 

sec 

910.34 

meters 

(4  its) 

3284.2 

meters 

The  results  in  Table  27  indicate  the  importance  of  modeling  shallow  resonance  terms;  adding  the  11th  and  12th 
order  shallow  resonance  terms  reduces  the  predict  error  by  roughly  200  meters.  The  short  periodic  models  for  drag, 
SRP,  and  third  body  have  negligible  accuracy  contributions,  but  each  add  a  few  seconds  of  run  time;  therefore,  these 
short  periodic  models  can  be  neglected  to  reduce  computation  time.  Like  many  other  cases,  the  tesseral  high 
frequency  short  periodic  and  J2  /  m-daily  short  periodic  terms  serve  to  provide  closure  with  Cowell,  but  are 
computationally  expensive.  For  these  orbits,  their  effects  can  be  neglected. 

These  results  also  highlight  the  contribution  of  several  of  the  AOG  terms: 

•  not  modeling  drag  in  the  AOG  (as  well  as  not  solving  for  Cp)  adds  roughly  100  meters  worth  of  predict  error 

•  not  modeling  solar  radiation  pressure  in  the  AOG  (as  well  as  not  solving  for  Cr)  adds  roughly  50  meters  worth 
of  predict  error 

•  not  modeling  third  body  effects  in  the  AOG  adds  roughly  50  meters  worth  of  predict  error 

In  the  first  high  accuracy  DSST  run,  the  DC  RMS  error  is  on  the  order  of  40  meters.  This  error  can  mainly  be 
attributed  to  truncations  to  the  equations  of  motion  in  DSST’s  J22  model  (Ofe1]  in  the  AOG  and  O[e0]  in  the  SPG),  an 
incomplete  tesseral  high  frequency  short  periodic  model,  and  discrepancies  between  the  truth  and  DC  “solved-fof’ 
values  of  CD  and  CR.  If  this  same  DSST  case  is  fit  to  another  Cowell  truth  model  containing  21x21  geopotential 
effects  only,  the  DC  RMS  is  reduced  to  -  21  meters  (in  other  words,  -  20  meters  can  be  attributed  to  the  poor  DC 
solution  for  CD  and  Cr).  If  this  test  is  taken  one  step  further  such  that  both  the  21x21  Cowell  truth  and  DSST  use  a 
small  value  of  J2  (1.0  D’6),  the  error  is  reduced  to  ~  1.7  meters  (i.e.,  errors  due  to  truncations  in  DSST’s  J22  model  are 
on  the  order  of  -  20  meters).  This  means  ~l-2  meters  of  error  can  be  attributed  to  not  having  a  complete  21x21 
tesseral  high  frequency  short  periodic  model  (remember,  an  SPTESSLC  of  21  21  21  15  -37  37  was  used).  Finally,  it 
should  be  mentioned  this  specific  breakdown  of  errors  is  dependent  upon  the  geometry  of  this  particular  orbit,  the 
relationship  between  the  fit  span  and  perigee  locations,  etc.  Therefore,  these  error  quantifications  should  be  used  as  a 
“rule  of  thumb”  metrics. 


Again,  it  is  of  particular  importance  to  use  extreme  care  when  running  an  optimized  version  of  the  DSST.  Note 
how  the  optimized  deck  for  this  case  implements  a  21x21  geopotential  field  with  a  RESONPRD  card  set  to  one  day 
and  truncated  short  periodics  (rather  than  the  RESNM  and  AUTOFORC  combination).  Again,  some  accuracy  and 
run-time  is  sacrificed  for  a  deck  which  is  (1)  simpler  to  implement  than  the  RESNM  and  AUTOFORC  combination 
and  (2)  general  enough  to  capture  resonant  teims  which  might  be  experienced  by  other  satellites  under  the  broad 
definition  of  the  NEEC. 

The  setup  for  the  NEEC  which  provides  an  optimal  balance  between  accuracy  and  computational  speed  can  now 
be  given: 
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•  21x21  geopotential  (AOG) 

•  RESONPRD  equal  to  one  day  (86400.0  seconds) 

•  J22  (AOG  and  SPG) 

•  Truncated  zonal  short  periodics  (SPZONALS  =  6  5  13)  . 

•  Truncated  m-daily  short  periodics  (SPMDAILY  =  664) 

•  Lunar  /  solar  third  body  point  mass  effects  (AOG) 

•  Atmospheric  drag  (AOG) 

•  Solve  for  CD 

•  Solar  radiation  pressure  (AOG) 

•  Solve  for  CR 

•  SSTESTFL  1  20  0.0 

•  SSTAPGFL  1  0  0  1.0  0.0  1.0 

Sample  input  decks  for  the  NEEC  are  given  in  the  appendices. 

8.0  Conclusions 

This  paper  describes  R&D  GTDS  input  decks  which  provide  a  balance  between  speed  and  accuracy  when  using 
the  DSST  orbit  propagation  theory.  The  orbit  classes  studied  in  this  effort  are  similar  to  those  investigated  in  the 
Ephemeris  Theory  Accuracy  Study  (ETAS)  performed  by  Kaman  Sciences5,  and  include  low,  medium,  and  high 
altitude  circular  orbits,  as  well  as  Molniya,  geosynchronous,  and  near  earth  eccentric  orbits.  In  addition  to  the 
optimized  decks  (refer  to  Appendix  A),  standardized  test  cases  have  been  supplied  for  each  of  the  orbit  types 
analyzed  herein  (refer  to  Appendix  B).  All  results  were  obtained  from  batch  differential  correction  fits  and 
subsequent  ephemeris  comparisons  to  Cowell  generated  truth  data.  The  input  decks  specified  in  this  study  provide  the 
R&D  GTDS  user  with  insight  concerning  the  accuracy  and  computational  expense  of  the  various  force  model  options 
(AOG  &  SPG)  available  in  DSST  for  each  orbital  class.  The  optimal  decks  also  provide  the  basis  for  an  automatic 
force  model  option  which  can  be  added  to  current  DSST  implementations. 

9.0  Recommendations 

Clearly,  this  paper  does  not  address  all  of  the  orbit  types  currently  in  use.  Therefore,  if  optimal  DSST  input  decks  are 
desired  for  other  orbit  types,  further  research  is  required.  In  addition,  software  modifications  should  be  made  to  R&D 
GTDS  to  allow  the  user  to  automatically  implement  these  optimal  input  decks  (if  so  desired).  These  setups  could  be 
specified  on  the  SPSHPER  keyword  card.  Finally,  an  automatic  resonance  selection  capability  (based  on  the  input 
semimajor  axis)  should  be  considered  for  R&D  GTDS. 
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Appendix  A 


Optimized  DSST  Input  Decks 
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LACC  OPTIMIZED  INPUT  DECK 


CONTROL 

DC 

SATNAME 

EPOCH 

820224.0 

0.0 

ELEMENT 1 

1 

6 

1 

6628.457 

0.0089 

64.84 

ELEMENT2 

224.51 

271.89 

115.16 

OBSINPUT 

15 

820224000000.0 

820224013000.0 

ORB TYPE 

5 

1 

1 

43200.0 

1.0 

DMOPT 

OBSDEV 

21 

22 

23 

100. 

100. 

100. 

OBSDEV 

24 

25 

26 

o 

H 

10. 

10. 

END 

DCOPT 

PRINTOUT 

1 

4 

CONVERG 

30 

1 

l.D-4 

END 

OGOPT 

NCBODY 

1 

DRAG 

1 

1.0 

ATMOS DEN 

1 

SCPARAM 

l.D-6 

100. DO 

SPGRVFRC 

1 

1 

1 

3.0 

1.0 

3.0 

SPTESSLC 

6 

6 

4 

2.0 

-10.0 

10.0 

POTFIELD 

1 

6 

MAXDEGEQ 

1 

8. 

MAXORDEQ 

1 

8. 

STATEPAR 

3 

STATETAB 

1 

2 

3 

4.0 

5.0 

6.0 

DRAGPAR 

1 

DRAGPAR2 

1 

1 

SSTESTFL 

1 

2 

0 

o 

o 

SSTAPGFL 

1 

0 

0 

1.0 

6.0 

0.0 

END 

FIN 

CONTROL 

EPHEM 

OUTPUT 

SATNAME 

OUTPUT 

1 

2 

1 

820224. 

030000.0 

86400.0 

ORBTYPE 

5 

1 

1 

43200.0 

1.0 

OGOPT 

NCBODY 

1 

DRAGPAR 

0 

DRAG 

1 

1.0 

ATMOSDEN 

1 

SCPARAM 

l.D-6 

100. DO 

SPGRVFRC 

1 

1 

1 

3.0 

1.0 

3.0 

SPTESSLC 

6 

6 

4 

2.0 

o 

o 

H 

1 

10.0 

POTFIELD 

1 

6 

MAXDEGEQ 

1 

8. 

MAXORDEQ 

1 

8. 

OUTOPT 

21 

820224000000.0 

820224030000.0 

450.0 

END 

FIN 

1234567 


1234567 
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MACC  OPTIMIZED  INPUT  DECK 


CONTROL 

DC 

LNDSAT-4 

EPOCH 

820224.0 

0.0 

ELEMENT 1 

3 

6 

1 

7077.8 

0.0011 

98.2 

ELEMENT2 

158.1 

89.4 

176.0 

OBS INPUT 

15 

820224000000.0 

820227000000.0 

ORB  TYPE 

5 

1 

1 

43200.0 

1.0 

DMOPT 

OBSDEV 

21 

22 

23 

100. 

100. 

100. 

OBSDEV 

24 

25 

26 

10. 

10. 

10. 

END 

DCOPT 

PRINTOUT 

1 

4 

CONVERG 

30 

1 

l.D-4 

END 

OGOPT 

DRAG 

1 

1 

ATMOS DEN 

1 

SPDRAG 

0 

SCPARAM 

l.D-6 

100. DO 

SPGRVFRC 

1 

1 

3 

3.0 

3.0 

3.0 

SPZONALS 

8 

2 

11 

SPMDAILY 

8 

8 

2 

POTFIELD 

1 

6 

MAXDEGEQ 

1 

21. 

MAXORDEQ 

1 

21. 

RESONPRD 

86400.0 

STATEPAR 

3 

STATETAB 

1 

2 

3 

4.0 

5.0 

6.0 

DRAGPAR 

1 

DRAGPAR2 

1 

1 

SSTESTFL 

1 

2 

SSTAPGFL 

1 

0 

0 

1.0 

0.0 

o 

o 

END 

FIN 

CONTROL 

EPHEM 

OUTPUT 

LNDSAT-4 

OUTPUT 

1 

2 

1 

820302. 

0.0 

86400.0 

ORBTYPE . 

5 

1 

1 

43200.0 

1.0 

OGOPT 

DRAG 

1 

1 

DRAGPAR 

0 

v 

ATMOSDEN 

1 

SPDRAG 

0 

SCPARAM 

l.D-6 

100. DO 

SPGRVFRC 

1 

1 

3 

3.0 

3.0 

3.0 

SPZONALS 

8 

2 

11 

SPMDAILY 

8 

8 

2 

POTFIELD 

1 

6 

MAXDEGEQ 

1 

21. 

MAXORDEQ 

1 

21. 

RESONPRD 

86400.0 

OUTOPT 

21 

820224000000.0 

820302000000.0 

450.0 

END 

FIN 

8207201 


8207201 
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HACC  OPTIMIZED  INPUT  DECK 


CONTROL 

DC 

EXPLORER 

EPOCH 

940826.0 

074853.5101 

ELEMENT 1 

3 

6 

1 

7498.0 

0.0247108 

41.1904 

ELEMENT2 

349.3346 

66.5461 

296.1065 

OBS INPUT 

15 

940826074854.0 

940828074854.0 

ORB TYPE 

5 

1 

1 

43200.0 

1.0 

DMOPT 

OBSDEV 

21 

22 

23 

100. 

100. 

100. 

OBSDEV 

24 

25 

26 

10. 

10. 

10. 

END 

DCOPT 

PRINTOUT 

1 

4 

CONVERG 

30 

1 

l.D-4 

END 

OGOPT 

SOLRAD 

1 

1.0 

SPSRP 

0 

SCPARAM 

l.D-6 

100. DO 

SPGRVFRC 

1 

1 

3 

3.0 

3.0 

3.0 

SPZONALS 

8 

2 

11 

SPMDAILY 

8 

8 

2 

POTFIELD 

1 

2 

MAXDEGEQ 

1 

21. 

MAXORDEQ 

1 

21. 

RESONPRD 

86400.0 

STATEPAR 

3 

STATETAB 

1 

2 

3 

4.0 

5.0 

6.0 

SOLRDPAR 

1 

SSTESTFL 

1 

2 

0 

o 

o 

SSTAPGFL 

T 

0 

0 

0.0 

0.0 

1.0 

END 

.  ‘ 

FIN 

CONTROL 

EPHEM 

OUTPUT 

EXPLORER 

OUTPUT 

1 

2 

1 

940831. 

074854.0 

86400.0 

ORB TYPE 

5 

1 

1 

43200.0 

1.0 

OGOPT 

SOLRAD 

1 

1.0 

SPSRP 

0 

SCPARAM 

l.D-6 

100. DO 

SOLRDPAR 

0 

SPGRVFRC 

1 

1 

3 

3.0 

3.0 

3.0 

SPZONALS 

8 

2 

11 

SPMDAILY 

8 

8 

2 

POTFIELD 

1 

2 

MAXDEGEQ 

1 

21. 

MAXORDEQ 

1 

21. 

RESONPRD 

86400.0 

OUTOPT 

21 

940826074854.0 

940831074854.0 

450.0 

END 

FIN 

65032A 


65032A 


4 
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GEOSYNCHRONOUS  OPTIMIZED  INPUT  DECK 


CONTROL 

DC 

BS3A 

EPOCH 

940825.0 

170334.9932 

ELEMENT 1 

1 

6 

1 

42167.16 

0.000275 

0.23 

ELEMENT2 

97.6 

75.3 

166.09 

OBS INPUT 

15 

940825170335.0 

940828170335.0 

ORB TYPE 

5 

1 

1 

43200.0 

1.0 

DMOPT 

OBSDEV 

21 

22 

23 

100. 

100. 

100. 

OBSDEV 

24 

25 

26 

10. 

10. 

10. 

END 

DCOPT 

PRINTOUT 

1 

4 

CONVERG 

30 

1 

l.D-4 

END 

OGOPT 

SOLRAD 

1 

1.0 

SCPARAM 

l.D-6 

100. DO 

SPGRVFRC 

1 

3 

3 

2.0 

1.0 

3.0 

SPNUMGRV 

7 

1 

8 

2.0 

2.0 

3600.0 

POTFIELD 

1 

2 

MAXDEGEQ 

1 

4. 

MAXORDEQ 

1 

4. 

STATEPAR 

3 

STATETAB 

1 

2 

3 

4.0 

5.0 

6.0 

SOLRDPAR 

1 

SSTESTFL 

1 

2 

0 

0.0 

SSTAPGFL 

1 

0 

0 

0.0 

0.0 

1.0 

END 

FIN 

CONTROL 

EPHEM 

OUTPUT 

BS3A 

OUTPUT 

1 

2 

1 

940905. 

170335.0 

86400. 

ORB TYPE 

5 

1 

1 

43200.0 

1.0 

OGOPT 

SOLRAD 

1 

1.0 

SOLRDPAR 

0 

SCPARAM 

l.D-6 

100. DO 

SPGRVFRC 

1 

3 

3 

2.0 

1.0 

3.0 

SPNUMGRV 

7 

1 

8 

2.0 

2.0 

3600.0 

POTFIELD 

1 

2 

MAXDEGEQ 

1 

4  . 

MAXORDEQ 

1 

4. 

OUTOPT 

21 

940825170335.0 

940905170335.0 

450.0 

END 

FIN 


20771 


20771 
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MOLNIYA  OPTIMIZED  INPUT  DECK 


CONTROL 

DATAMGT 

NSSC 

OGOPT 

POTFIELD 

1 

6 

END 

FIN 

CONTROL 

DC 

NSSC 

EPOCH 

790804.0 

234212.0 

ELEMENT1 

1 

6 

1 

26572.176 

0.699 

63  .2 

ELEMENT2 

190.5 

281.6 

15.429 

OBS INPUT 

15 

790804234212.0 

790807234212.0 

ORB TYPE 

5 

1 

1 

43200.0 

1.0 

DMOPT 

OBSDEV 

21 

22 

23 

100. 

100. 

100. 

OBSDEV 

24 

25 

26 

10. 

10. 

10. 

END 

OGOPT 

DRAG 

1 

1.0 

ATMOSDEN 

1 

SPDRAG 

0 

SCPARAM 

12.5D-6 

1250. DO 

SOLRAD 

1 

1.0 

SPSRP 

0 

SPGRVFRC 

1 

1 

3 

1. 

1. 

3. 

SPZONALS 

4 

3 

8 

SPMDAILY 

4 

4 

2 

MAXDEGEQ 

1 

8. 

MAXORDEQ 

1 

8. 

STATEPAR 

3 

STATETAB 

1 

2 

3 

4.0 

5.0 

6.0 

DRAGPAR 

1 

DRAGPAR2 

1 

:  1 

SOLRDPAR 

1 

SSTESTFL 

1 

2 

0 

0.0 

SSTAPGFL 

1 

0 

0 

1.0 

o 

o 

o 

H 

END 

DCOPT 

PRINTOUT 

1 

4 

CONVERG 

30 

1 

l.D-3 

END 

FIN 

CONTROL 

EPHEM 

OUTPUT 

NSSC 

OUTPUT 

1 

2 

1 

790818.0 

234212.0 

172800 

ORB TYPE 

5 

1 

1 

43200.0 

1.0 

OGOPT 

DRAGPAR 

0 

DRAG 

1 

H 

O 

ATMOSDEN 

1 

SPDRAG 

0 

SCPARAM 

12.5D-6 

1250. DO 

SOLRAD 

1 

1.0 

SPSRP 

0 

SOLRDPAR 

0 

SPGRVFRC 

1 

1 

3 

1. 

1. 

3. 

SPZONALS 

4 

3 

8 

SPMDAILY 

4 

4 

2 

MAXDEGEQ 

1 

8. 

MAXORDEQ 

1 

8. 

OUTOPT 

21 

790804234212.0 

790818234212.0 

450. 

END 

FIN 

9829 


9829 


9829 
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NEEC  OPTIMIZED  INPUT  DECK 


CONTROL 

DC 

VANGARD2 

EPOCH 

940826.0 

073513.6735 

ELEMENT 1 

3 

6 

1 

8177.913 

0.1522640 

32.8834 

ELEMENT 2 

251.8592 

10.8368 

352.1515 

OBS INPUT 

15 

940826073514.0 

940829073514.0 

ORBTYPE 

5 

1 

1 

43200.0 

1.0 

DMOPT 

OBSDEV 

21 

22 

23 

100. 

100. 

100. 

OBSDEV 

24 

25 

26 

10. 

10. 

10. 

END 

DCOPT 

PRINTOUT 

1 

4 

CONVERG 

30 

1 

l.D-4 

END 

OGOPT 

DRAG 

1 

1.0 

ATMOSDEN 

1 

SPDRAG 

0 

SCPARAM 

l.D-6 

100. DO 

SOLRAD 

1 

1.0 

SPSRP 

0 

SPGRVFRC 

1 

1 

3 

3.0 

1.0 

3.0 

SPZONALS 

6 

5 

13 

SPMDAILY 

6 

6 

4 

POTFIELD 

1 

2 

MAXDEGEQ 

1 

21. 

MAXORDEQ 

1 

21. 

RESONPRD 

86400.0 

STATEPAR 

3 

STATETAB 

1 

2 

3 

4.0 

5.0 

6.0 

DRAGPAR 

1 

,  v  .. 

DRAGPAR2 

1 

1 

SOLRDPAR 

1 

SSTESTFL 

1 

2 

0 

0.0 

SSTAPGFL 

1. 

0 

0 

1.0 

0.0 

1.0 

END 

FIN 

CONTROL 

EPHEM 

OUTPUT 

VANGARD2 

OUTPUT 

1 

2 

1 

940905. 

073514.0 

86400.0 

ORBTYPE 

5 

1 

1 

43200.0 

1.0 

OGOPT 

DRAGPAR 

0 

DRAG 

1 

1.0 

ATMOSDEN 

1 

SPDRAG 

0 

SCPARAM 

l.D-6 

100. DO 

SOLRAD 

1 

1.0 

SPSRP 

0 

SOLRDPAR 

0 

SPGRVFRC 

1 

1 

3 

3.0 

1.0 

3.0 

SPZONALS 

6 

5 

13 

SPMDAILY 

6 

6 

4 

POTFIELD 

1 

2 

MAXDEGEQ 

1 

21. 

MAXORDEQ 

1 

21. 

RESONPRD 

86400.0 

OUTOPT 

21 

940826073514.0 

940905073514.0 

450.0 

59001A 


5900 1 A 
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Appendix  B 


Optimal  Deck  Test  Case  Summaries 


This  appendix  details  information  required  to  reproduce  results  for  each  of  the  optimal  DSST  input  deck  test 
cases  described  in  this  document.  Specifically,  the  following  information  is  given: 


•  The  input  decks  required  to  generate  the  Cowell  truth  orbit.  For  cases  starting  with  a  two  card  elements  set, 
SGP4  ephemeris  generation  decks  arc  given,  along  with  the  Cowell  DC  and  subsequent  Cowell  ephemeris 
generation  decks  that  were  used  to  fit  the  SGP4  ephemeris  and  produce  the  truth  trajectory.  For  cases  starting 
with  DSST  mean  Keplerian  elements,  the  DSST  ephemeris  generation  decks  used  to  convert  the  mean  elements 
to  osculating  elements  are  given,  along  with  the  Cowell  ephemeris  generation  decks  which  used  the  osculating 
elements  as  inputs  to  produce  a  truth  trajectory.  In  all  cases,  output  summaries  are  included  to  verify  the 
reproduction  of  truth  data. 


•  The  optimal  DSST  DC  input  decks.  Vital  statistics  from  the  DC  output  report  are  included  to  verify 
reproduction  of  the  test  cases.  GTDS  input  decks  and  output  summaries  are  also  given  for  the  comparison  of  the 
predicted  “optimal”  ephemerides  to  the  truth  data. 
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LACC  Test  Case 
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TRUTH  DATA  GENERATION 


PROCEDURE : 

1)  Start  with  mean  Keplerian  elements 

2)  Propagate  forward  one  day  using  DSST  to  obtain  osculating  elements 

3)  Propagate  osculating  elements  forward  using  Cowell  theory  to  generate  truth  data 


CONTROL 

EPHEM 

SATNAMB 

EPOCH 

820223.0 

o 

o 

ELEMENT1 

3 

6 

1 

6635.0814 

0.010201164 

64.9567 

ELEMENT2 

228.6393 

271.2229 

88.164558 

OUTPUT 

1 

2 

1 

820224.0 

o 

o 

43200.0 

ORB TYPE 

5 

1 

1 

43200.0 

1.0 

OGOPT 

SPSHPER 

3 

DRAG 

1 

1.0 

ATMOSDEN 

1 

SCPARAM 

l.D-6 

100. DO 

POTFIELD 

1 

6 

MAXDEGEQ 

1 

8. 

MAXORDEQ 

1 

8. 

END 

FIN 

CONTROL 

EPHEM 

OUTPUT 

SATNAME 

OUTPUT 

1 

2 

1 

820302.0 

0.0 

86400.0 

ORE TYPE 

2 

1 

1 

60.0 

OGOPT 

DRAG 

1 

1.0 

ATMOSDEN 

1 

SCPARAM 

l.D-6 

100. DO 

POTFIELD 

1 

6 

MAXDEGEQ 

1 

21.0 

MAXORDEQ 

1 

21.0 

OUTOPT 

1 

820224000000.0 

820302000000.0 

450.0 

TRUTH  DATA  GENERATION 


COWELL  EPHEMERIS  OUTPUT  SUMMARY: 
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GTDS  FINAL  REPORT 


SATELLITE  NAME  SATNAME 

SATELLITE  NUMBER  1234567 


RUN  REFERENCE  DATE 

FEB 

23,  1982 

0  HRS 

0  MINS 

0.00000  SECONDS 

RUN  EPOCH  DATE 

FEB 

24,  1982 

0  HRS 

0  MINS 

0.00000  SECONDS 

RUN  FINAL  TIME 

FEB 

27,  1982 

8  HRS 

30  MINS 

0.00000  SECONDS 

TOTAL  TIME  OF  FLIGHT 

3  DAYS 

8  HRS 

30  MINS 

0.00000  SECONDS 

CAUSE  OF  TERMINATION  IMPACT  ON  CENTRAL  BODY 


CENTRAL  BODY  IS  EARTH  (INERTIAL 
X  0.5593 19 986957752 9E+04 

VX  -0 . 1067207195712157E+01 

SMA  0 . 5566864380626048E+04 

LAN  0 . 2119745 94 2 73 5 612E+ 03 

EA  0 . 1961928616234382E+03 

PR  0 . 4680456054501052E+04 

APH  0 . 751347 0675104418E+02 

RA  0 . 2 856256 3544 54737E+02 

AZ  0 . 1543948811906195E+03 


***END  CONDITIONS*** 

SYSTEM) 

Y  0 . 304477235863 8 768E+0 4 

VY  0 .29843 83598 9 04645E+01 

ECC  0 . 1592293 73 2 19491 0E+00 

AP  0.3382687050241062E+03 

P  0 . 114821 82424 9102 0E+01 

APR  0.6453272706751 04 4E+04 

C3  -0 . 3580116316352344E+02 

DEC  0 . 71467714 99993 070E+ 01 

RMAG  0 . 6413107503189311E+04 


MEAN  OF  1950.0  --  EARTH  EQUATOR 

Z  0 . 79 84 86333 9 60 7 9 82E+0 3 
VZ  “0 . 65240 63462834691E+01 

INC  0 . 64 60 801957824 047E+02 
MA  0 .198737053 68 08 2 72E+03 
SLR  0 . 542572213 844 06 27E+04 
PH  -0 . 1 6 97 6 8 194 54 9 8 94 8E+ 04 
TA  0 . 1938156243439430E+03 

VPA  0 . 925753355 7949 9 78E+02 
VMAG  0 .7253 197965 8 00011E+01 


*** INITIAL  CONDITIONS*** 

CENTRAL  BODY  IS  EARTH  (INERTIAL  SYSTEM)  MEAN  OF  1950.0  --  EARTH  EQUATOR 


X  -0.3260654138207152E+04 

VX  0.4586987250728 03 3E+01 

SMA  0 . 662 845721232248 8E+04 

LAN  0 . 2 245 064 3 24 6050 5 7E+ 03 

EA  0 . 11562294263 111 8 8E+ 03 

PR  0 . 65693 129072 5 0206E+04 

APH  0 . 3094635173947700E+03 

RA  0 . 2372324679669332E+03 

AZ  0 . 28008 9735645 8445E+02 


Y  -0 . 5065846974834642E+04 
VY  0 . 4241328473629827E+00 

ECC  0 . 8922785978361727E-02 
AP  0 . 271 892244504 8535E+03 
P  0 . 1491856897264069E+01 

APR  0 . 66 676015 173 94770E+04 
C3  -0 . 3006736162217285E+02 

DEC  0 . 2512 3 93 544617 055E+02 
RMAG  0 . 6654033 9 7956044 8E+ 04 


Z  0 . 2825154370949725E+04 

VZ  0 . 6200975003 2012 67E+01 

INC  0 . 64 83 7894 25 043 53 3E+02 

MA  0 . 115161980 083 9504E+03 

SLR  0 . 6 627 92 94 6 03 4 64 8 9E+ 04 

PH  0 . 191174 9 0725 02 06 2E+03 

TA  0 . 1160830242166191E+03 

VPA  0 . 89539029 04 810165E+02 

VMAG  0 . 7724793311207796E+01 


***  SECTIONING  SUMMARY  *** 

NUMBER  OF  SECTIONS  SCHEDULED  ■  1 

NUMBER  OF  SECTIONS  COMPLETED  =  1 

SECTION  CENTRAL  BODY  TIME  OF  CROSSING  TIME  INTO  FLIGHT  CAUSE  OF  CROSSING 

H  M  S  H  M  S 

1  EARTH  FEB  27  8  30  0.000  80  30  0.000  IMPACT  ON  CENTRAL  BODY 


***  FILE  GENERATION  SUMMARY  *** 

AN  ORB-1  FILE  HAS  BEEN  GENERATED 

START  TIME  OF  THE  FILE  **  820224000000.000 

END  TIME  OF  THE  FILE  »  820227083000.000 


***  NUMBER  OF  GTDS  ERRORS  ENCOUNTERED**  0** 
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LACC  OPTIMIZED  INPUT  DECK 


CONTROL 

DC 

SATNAME 

EPOCH 

820224.0 

0.0 

ELEMENT1 

1 

6 

1 

6628.457 

0.0089 

64.64 

ELEMENT2 

224.51 

271.89 

115.16 

OBS INPUT 

15 

820224000000.0 

820224013000.0 

ORBTYPE 

5 

1 

1 

43200.0 

1.0 

DMOPT 

OBSDEV 

21 

22 

23 

100. 

100. 

100. 

OBSDEV 

24 

25 

26 

10. 

10. 

10. 

END 

DCOPT 

PRINTOUT 

1 

4 

CONVERG 

30 

1 

l.D-4 

END 

OGOPT 

NCBODY 

1 

DRAG 

1 

1.0 

ATMOSDEN 

1 

SCPARAM 

l.D-6 

100. DO 

SPGRVFRC 

1 

1 

1 

3.0 

1.0 

3.0 

SPTESSLC 

6 

6 

4 

2.0 

o 

o 

H 

1 

10.0 

POTFIELD 

1 

6 

MAXDEGEQ 

1 

8. 

MAXORDEQ 

1 

8. 

STATEPAR 

3 

STATETAB 

1 

2 

3 

4.0 

5.0 

6.0 

DRAG PAR 

1 

DRAGPAR2 

1 

1 

SSTESTFL 

1 

2 

0 

0.0 

SSTAPGFL 

1 

0 

0 

1.0 

6.0 

0.0 

END 

FIN 

CONTROL 

EPHEM 

OUTPUT 

SATNAME 

OUTPUT 

1 

2 

1 

820224. 

030000.0 

86400.0 

ORBTYPE 

5 

1 

1 

43200.0 

1.0 

OGOPT 

NCBODY 

1 

DRAGPAR 

0 

DRAG 

1 

1.0 

ATMOSDEN 

1 

SCPARAM 

l.D-6 

100. DO 

SPGRVFRC 

i 

1 

1 

3.0 

1.0.  - 

3.0 

SPTESSLC 

6 

6 

4 

2.0 

o 

o 

10.0 

POTFIELD 

1 

6 

MAXDEGEQ 

1 

8. 

MAXORDEQ 

1 

8. 

OUTOPT 

21 

820224000000.0 

820224030000.0 

450.0 

END 

FIN 


1234567 


1234567 
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SUMMARY  OF  LACC  OPTIMIZED  DSST  FIT  TO  COWELL  TRUTH: 


SATELLITE  SATNAME  1234567 

EPOCH  820224  0.000 

COORD.  SYSTEM  MEAN  OF  1950.0 

CENTRAL  BODY  EARTH 


ITERATION  NUMBER  6 

*********************************************** 

*  POSITION  ERROR  RMS  (M)  17.404205  * 

*  CURRENT  WEIGHTED  RMS  0.11590288  *  START  TIME  820224  0.000 

*  PREDICTED  WEIGHTED  RMS  0.11590311  * 

*  PREVIOUS  WEIGHTED  RMS  0.11590339  *  END  TIME  820224  12230.000 

*  SMALLEST  WEIGHTED  RMS  0.11590339  * 

*  RELATIVE  CHANGE  IN  RMS  4 .40389818E-06  *  NO.  OBS .  AVAILABLE  72 

*  PENALTY  0.00000000  * 

*  DC  HAS  CONVERGED  *  NO.  OBS.  INCLUDED  72 

*  * 

***********************************************  NO.  OBS.  ACCEPTED  72  100  PCT. 
OBSERVATION  SUMMARY  BY  TYPE 


TYPE 

TOTAL  NO. 

NO.  ACCEPTED 
WEIGHTED  RMS 
MEAN  RESDUAL 
STANDARD  DEV 


X  Y  Z  XDOT  YDOT  ZDOT 

12  12  12  12  12  12 

12  (100%)  12  (100V)  12  (1001)  12  (100%)  12  (100%)  12  (100%) 

7 . 0606E-02  0.1319  8.8872E-02  0.1270  0.1428  0.1152 

0.5897  -0.1745  1.190  -0.4787  0.7263  0.2646 

7.036  13.19  8.807  1.176  1.229  1.122 


0 INTERVAL 
TOTAL  NO. 

NO.  ACCEPTED 
(100%) 

0 INTERVAL 
TOTAL  NO. 

NO.  ACCEPTED 
(100%) 


00-20  20-40 

0  6 
0(0%)  6  (100%) 


OBSERVATION  SUMMARY  BY  20  DEGREE  TRUE  ANOMALY  INTERVAL 
40-60'  60-80  80-100  100-120  120-140  140-160  160-180 

6  0  6  6  0  6  6 

6  (100%)  0(0%)  6  (100%)  6  (100%)  0(0%)  6  (100%)  6 


180-200  200-220  220-240  240-260  260-280  280-300  300-320  320-340  340-360 

066066066 
0(0%)  6  (100%)  6  (100%)  0(0%)  6  (100%)  6  (100%)  0(0%)  6  (100%)  6 


THERE  ARE  7  SOLVE  PARAMETERS 


NUMBER  OF  DYNAMIC  SOLVE  PARAMETERS  =  7 

A  jl  K  P  Q  LAM  CSUBD 

NUMBER  OF  STATION  LOCATIONS  BEING  SOLVED  FOR  *  0 

NUMBER  OF  BIASES  BEING  SOLVED  FOR  *  0 


***  FILE  GENERATION  SUMMARY  *** 
NO  FILES  HAVE  BEEN  GENERATED 

*  NUMBER  OF  GTDS  ERRORS  ENCOUNTERED  TO  DATE  »  0  * 


60 


COMPARE  DECK  FOR  LACC; 


CONTROL  COMPARE  SATNAME  1234567 

COMPOPT 

CMPEPHEM  1102102  820224013000.0  820224030000.0  7.5 

END 

FIN 


LACC  COMPARE  SUMMARY: 


EPHEMERIS  COMPARISON  SUMMARY  REPORT 
Number  of  Points  Compared  13 


MINIMUM  POSITION  DIFFERENCE  MAXIMUM  POSITION  DIFFERENCE 


YYMMDD 

HHMMSS.SSS 

0cm) 

YYMMDD 

HHMMSS.SSS 

(km) 

RADIAL 

620224 

24500.000 

8 . 498042E-04 

820224 

21500.000 

2 . 135948E-02 

CROSS  TRACK 

820224 

21500.000 

1 . 065784E-02 

820224 

15230.000 

4.588219E-02 

ALONG  TRACK 

820224 

13000.000 

5 . 055184E-03 

820224 

30000.000 

1.013102E-01 

TOTAL 

820224 

13000.000 

1 . 904066E-02 

820224 

30000.000 

1 . 052537E-01 

MINIMUM  VELOCITY  DIFFERENCE  MAXIMUM  VELOCITY  DIFFERENCE 


YYMMDD 

HHMMSS.SSS 

(km/sec) 

YYMMDD 

HHMMSS.SSS 

(km/sec) 

RADIAL 

820224 

22230.000 

1.024987E-05 

820224 

30000.000 

1.508931E-04 

CROSS  TRACK 

820224 

24500.000 

3.131833E-06 

820224 

13730.000 

5 . 871673E-05 

ALONG  TRACK 

820224 

24500.000 

1.064040E-07 

820224 

13730.000 

2.326795E-05 

TOTAL 

820224 

20730.000 

4.330394E-05 

820224 

30000.000 

1 . 585244E- 04 

POSITION  RMS 

(km) 

VELOCITY  RMS 

(km/sec) 

RADIAL 

1 . 2699E-02 

7.2032E-0S 

CROSS  TRACK 

3 . 0955E-02 

3 . 6669E-05 

ALONG  TRACK 

5.6565E-02 

1 . 3224E-05 

TOTAL 

6.5720E-02 

GTDS  COMPARE  PROGRAM 

8.1903E-05 

PAGE 

NORMAL  COMPLETION  OF  JOB 
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MACC  Test  Case 
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TRUTH  DATA  GENERATION 

PROCEDURE : 

1)  Start  with  mean  Keplerian  elements 

2)  Propagate  forward  one  day  using  DSST  to  obtain  osculating  elements 

3)  propagate  osculating  elements  forward  using  Cowell  theory  to  generate  truth  data 


CONTROL  EPHEM 
EPOCH  820223*0 

ELEMENT1  361  7077.787 

ELEMENT2  1S8. 15349 

OUTPUT  121  820224.0 

ORBTYPE  511  43200.0 

OGOPT 

MAXDEGEQ  1  8 

MAXORDEQ  1  8 

DRAG  1  2 

SCPARAM  l.D-6 

POTFIELD  1  6 

END 

FIN 

CONTROL  EPHEM 
OUTPUT  1  2  1  820302.0 

ORBTYPE  2  1  1  60.0 

OGOPT 

DRAG  1  1 

ATMOSDEN  1 

MAXDEGEQ  1  21.0 

MAXORDEQ  1  21.0 

OUTOPT  1  820224000000.0 

SCPARAM  l.D-6 

POTFIELD  1  6 

END 

FIN 


LNDSAT-4  8207201 

0.0 

0.011542  98.250452 

89.4  312.90205 

0.0  43200.0 

1.0 


100. DO 


OUTPUT  LNDSAT-4  8207201 

0.0  86400.0 


820302000000.0  450.0 

100. DO 
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TRUTH  DATA  GENERATION 


COWELL  EPHEMERIS  OUTPUT  SUMMARY: 
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GTDS  FINAL  REPORT 


SATELLITE  NAME 

LNDSAT-4 

SATELLITE  NUMBER 

8207201 

RUN  REFERENCE  DATE 

FEB 

23, 

1982 

RUN  EPOCH  DATE 

FEB 

24, 

1982 

RUN  FINAL  TIME 

MAR 

2, 

1982 

TOTAL  TIME  OF  FLIGHT 

6  DAYS 

CAUSE  OF  TERMINATION 

SPECIFIED  TIME 

0 

HRS 

0  MINS 

0.00000 

SECONDS 

0 

HRS 

0  MINS 

0.00000 

SECONDS 

0 

HRS 

0  MINS 

0.00000 

SECONDS 

0 

HRS 

0  MINS 

0.00000 

SECONDS 

FLIGHT  REACHED 


CENTRAL  BODY  IS  EARTH  (INERTIAL  SYSTEM) 


X 

vx 


-0 . 63 24 74 09 6 8 94 2 S53E+ 04 
0 . 27099237242 89485E+ 01 


Y 

VY 


**END  CONDITIONS*** 

0 . 2 072 81 8 6 57 9943 22E+ 04 
0.32350 8674042329 6E+00 


MEAN  OF  1950.0  --  EARTH  EQUATOR 

Z  0 . 22 682174436793 85E+04 

VZ  0 . 7 04737 74 94 5 03 17 1E+01 


SMA  0 . 7064998645530947E+04 
IAN  0 . 1647012051719553E+03 
EA  0 . 3095957 802717 8 89E+03 
PR  0.7001256250617576E+04 
APH  0 .  790603 0404443181E+03 


ECC  0.1181967691217463E-01 
AP  0 . 6995308761596635E+02 
P  0.1648 610790422 61 8E+01 
APR  0. 716874104 04443 18E+ 04 
C3  -0 .2812988822 8802 13E+ 02 


INC  0.9829175379969750E+02 
MA  0 , 3 10117 6 17 19 3215 0E+0 3 
SLR  0 . 7084008 8374 79219E+04 
PH  0 . 6231182506175765E+03 
TA  0 . 3090719555792481E+03 


RA  0 . 1618543 854 217003E+03 
AZ  0 . 3512363878729003E+03 


DEC  0.18818 64 70769269 6E+02 
RMAG  0. 703 16 23986318 16 8E+ 04 


VPA  0 . 90 5 2185 8 94543 64 6E+ 02 
VMAG  0 . 7557372162575918E+01 


*** INITIAL  CONDITIONS*** 
CENTRAL  BODY  IS  EARTH  (INERTIAL  SYSTEM) 

X  0 . 2434985999885991E+04  Y  -0 . 1956779793707955E+04 

VX 

-0 . 6438412655705597E+01 

VY 

0 . 2019 610450 187177E+01 

SMA 

0 .7072044035079919E+04 

ECC 

0 . 12 6 53 14 83 59 77 568 E- 01 

LAN 

0 . 15875 64854 2 92450E+03 

AP 

0 . 8334407619753735E+02 

EA 

0 . 16 18 3 85110459 561E+03 

P 

0 . 1644 09 12 3 8 75 5 010E+ 01 

PR 

0.6982560412697187E+04 

APR 

0.7161527657462651E+04 

APH 

0.7833896574626 5 15E+03 

%  C3 

-0 . 2818141671791043E+02 

RA 

0 . 321214283 13 843 52E+0 3 

DEC 

-0 . 64121S1715921503E+02 

AZ 

0.1993584452782684E+03 

RMAG 

0.7157069741883023E+04 

MEAN  OF  1950.0  --  EARTH  EQUATOR 

Z  -0 . 6439371344275231E+04 
VZ  -0 . 3080855896509433E+01 


INC  0 . 983 1861805268 92 0E+02 
MA  0 . 16161253 99350 340E+03 
SLR  0 . 7 07 09 117 8 55 3014 0E+ 04 
PH  0 . 6044224126971867E+03 
TA  0 . 1620631405221004E+03 

VPA  0 . 8977401196957886E+02 
VMAG  0 . 7417793 2669 16529E+ 01 


***  SECTIONING  SUMMARY  *** 


SECTION  CENTRAL  BODY 


NUMBER  OF  SECTIONS  SCHEDULED 
NUMBER  OF  SECTIONS  COMPLETED 
TIME  OF  CROSSING 

H  M  S 

MAR  2  0  0  0.000 


1 

1 

TIME  INTO  FLIGHT 
H  M  S 
144  0  0.000 


CAUSE  OF  CROSSING 

SPECIFIED  TIME  OF  FLIGHT  REACHED 


***  FILE  GENERATION  SUMMARY  *** 


AN  ORB-1  FILE  HAS 
START  TIME  OF  THE 
END  TIME  OF  THE 


BEEN  GENERATED 

FILE  «  820224000000.000 

FILE  -  820302000000.000 


***  NUMBER  OF  GTDS  ERRORS  ENCOUNTERED  =  0** 
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MACC  OPTIMIZED  INPUT  DECK 


CONTROL 

DC 

LNDSAT-4 

EPOCH 

820224.0 

0.0 

ELEMENT1 

3 

6 

1 

7077.8 

0.0011 

98.2 

ELEMENT2 

158.1 

89.4 

176.0 

OBS INPUT 

15 

820224000000.0 

820227000000.0 

ORB TYPE 

5 

1 

1 

43200.0 

1.0 

DMOPT 

OBSDEV 

21 

22 

23 

100. 

100. 

100. 

OBSDEV 

24 

25 

26 

10. 

10. 

10. 

END 

DCOPT 

PRINTOUT 

1 

4 

CONVERG 

30 

1 

l.D-4 

END 

OGOPT 

DRAG 

1 

1 

ATMOSDEN 

1 

SPDRAG 

0 

SCPARAM 

l.D-6 

100. DO 

SPGRVFRC 

1 

1 

3 

3.0 

3.0 

3.0 

SPZONALS 

8 

2 

11 

SPMDAILY 

8 

8 

2 

POTFIELD 

1 

6 

MAXDEGEQ 

1 

21. 

MAXORDEQ 

1 

21. 

RESONPRD 

86400.0 

STATEPAR 

3 

STATETAB 

1 

2 

3 

4.0 

5.0 

6.0 

DRAG PAR 

1 

DRAG PAR 2 

1 

1 

SSTESTFL 

1 

2 

SSTAPGFL 

1 

0 

0 

1.0 

0.0 

0.0 

END 

FIN 

CONTROL 

EPHEM 

OUTPUT 

LNDSAT-4 

OUTPUT 

1 

2 

1 

820302. 

0.0 

86400.0 

ORB TYPE 

5 

1 

1 

43200.0 

1.0 

OGOPT 

DRAG 

1 

1 

DRAG PAR 

0 

ATMOSDEN 

1 

SPDRAG 

0 

^  : 

SCPARAM 

l.D-6 

100. DO 

SPGRVFRC 

1 

1 

3 

3.0 

3.0 

3.0 

SPZONALS 

8 

2 

11 

SPMDAILY 

8 

8 

2 

POTFIELD 

1 

6 

MAXDEGEQ 

1 

21. 

MAXORDEQ 

1 

21. 

RESONPRD 

86400.0 

OUTOPT 

21 

820224000000.0 

820302000000.0 

450.0 

END 

FIN 

8207201 


8207201 
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SUMMARY  OF  MACC  OPTIMIZED  DSST  FIT  TO  COWELL  TRUTH: 


ITERATION  NUMBER  7 


LNDSAT-4  8207201 


*********************************************** 

*  POSITION  ERROR  RMS  (M)  76.510799  * 

*  CURRENT  WEIGHTED  RMS  0.45418732  * 

*  PREDICTED  WEIGHTED  RMS  0.45580603  * 


START  TIME  820224  0.000 


EPOCH 

820224  0.000 

* 

PREVIOUS  WEIGHTED  RMS 

0.45582550 

* 

END 

TIME 

820226 

235230 

* 

SMALLEST  WEIGHTED  RMS 

0.45582550 

* 

COORD.  SYSTEM 

MEAN  OF  1950.0 

* 

RELATIVE  CHANGE  IN  RMS 

3 . 59387939E-03 

* 

NO. 

OBS. 

AVAILABLE 

3456 

* 

PENALTY 

0.00000000 

* 

CENTRAL  BODY 

EARTH 

* 

DC  HAS  CONVERGED 

* 

* 

NO. 

OBS. 

INCLUDED 

3456 

*********************************************** 

NO. 

OBS. 

ACCEPTED 

3423  99 

OBSERVATION  SUMMARY  BY  TYPE 


TYPE 

X 

Y 

Z 

XD0T 

YDOT 

ZDOT 

TOTAL  NO. 

576 

576 

576 

576 

576 

576 

NO.  ACCEPTED 

570  ( 

98%) 

576  (100%) 

575  ( 

99%) 

560  (  97%) 

575  (  99%) 

567  (  98%) 

WEIGHTED  RMS 

0.4544 

0.3822 

0.4861 

0.4854 

0.3506 

0.5407 

MEAN  RESDUAL 

0.3981 

0.2098 

0.4400 

9 . 6646E-02  - 

■3.4598E-02  - 

•3 . 1819E-02 

. STANDARD  DEV 

45.44 

38.22 

48.61 

4.853 

3.506 

5.407 

OBSERVATION  SUMMARY  BY  20  DEGREE 

TRUE  ANOMALY 

0 INTERVAL 

00-20 

20-40 

40-60 

60-80 

80-100 

100-120  : 

TOTAL  NO. 

156 

168 

198 

234 

204 

168 

NO.  ACCEPTED 

156  (100%) 

167  (  99%) 

197  ( 

99%) 

234  (100%) 

204  (100%) 

166  (  98%) 

98%} 

0 INTERVAL 

180-200 

200-220 

220-240 

240-260 

260-280 

280-300  : 

TOTAL  NO. 

240 

204 

174 

168 

174 

216 

NO.  ACCEPTED 

237  ( 

98%) 

202  (  99%) 

172  ( 

98%) 

165  (  98%) 

171  (  98%) 

213  (  98%) 

99%) 

THERE  ARE 

7  SOLVE  PARAMETERS 

NUMBER  OF  DYNAMIC  SOLVE 

PARAMETERS  * 

7 

A 

H 

K 

P 

Q 

LAM 

CSUBD 

120-140 

156 


300-320 

228 


140-160 

192 


320-340 

192 


160-180 

228 


340-360 

156 


NUMBER  OF  STATION  LOCATIONS  BEING  SOLVED  FOR 
NUMBER  OF  BIASES  BEING  SOLVED  FOR-  0 


***  FILE  GENERATION  SUMMARY  *** 
NO  FILES  HAVE  BEEN  GENERATED 

*  NUMBER  OF  GTDS  ERRORS  ENCOUNTERED  TO  DATE  -  0  * 
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COMPARE  DECK  FOR  MACC: 


CONTROL 

COMPOPT 

COMPARE 

LNDSAT-4 

8207201 

CMPEPHEM 

END 

FIN 

1102102  820227000000.0 

620302000000.0 

45. 

MACC  COMPARE  SUMMARY: 


EPHEMERIS  COMPARISON  SUMMARY  REPORT 
Number  of  Points  Compared  97 


MINIMUM  POSITION  DIFFERENCE  MAXIMUM  POSITION  DIFFERENCE 


YYMMDD 

HHMMSS.SSS 

(km) 

YYMMDD 

HHMMSS.SSS 

(km) 

RADIAL 

820228 

94500.000 

1 .  374394E-03 

820227 

184500.000 

9.115588E-02 

CROSS  TRACK 

820301 

141500.000 

5 . 420682E-04 

820301 

21500.000 

6 . 921743E-02 

ALONG  TRACK 

620227 

163000.000 

2.999669E-03 

820301 

193000.000 

3 . 680619E-01 

TOTAL 

820227 

163000.000 

1.286429E-02 

820301 

193000.000 

3.764689E-01 

MINIMUM  VELOCITY  DIFFERENCE 

MAXIMUM  VELOCITY 

DIFFERENCE 

YYMMDD 

HHMMSS.SSS 

(km/sec) 

YYMMDD 

HHMMSS.SSS 

(km/sec) 

RADIAL 

820228 

81500.000 

5.361974E-06 

820301 

124500.000 

3.909677E-04 

CROSS  TRACK 

820228 

21500.000 

5 . 331724E-07 

820227 

43000.000 

1 . 131149E-04 

ALONG  TRACK 

820228 

163000.000 

3 . 636120E-08 

820228 

73000.000 

9 . 822645E-05 

TOTAL 

820227 

133000.000 

2.279860E-05 

820301 

124500.000 

3 . 942775E-04 

POSITION  RMS 

(km) 

VELOCITY  RMS 

(km/sec) 

RADIAL 

4 . 0244E-02 

2 . 0643E-04 

CROSS  TRACK 

3 . 2424E-02 

3 .3844E-05 

ALONG  TRACK 

1 . 8994E-01 

4 . 4235E-05 

TOTAL 

1 . 9685E-01 

GTDS  COMPARE  PROGRAM 

2 . 1381E-04 

PAGE 

NORMAL  COMPLETION  OF  JOB 
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HACC  Test  Case 
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TRUTH  DATA  GENERATION 


PROCEDURE : 

1)  Start  with  2 -card  element  set 

2)  Propagate  2 -card  element  set  forward  using  SGP4  to  obtain  SGP4 -based  ephemeris 

3)  Fit  Cowell  theory  to  SGP4 -based  ephemeris 

4)  Propagate  Cowell  DC  solve- for  vector  forward  to  produce  truth  data 


CONTROL 

EPHEM 

EXPLORER 

65032A 

EPOCH 

940826.0 

074853.5101 

ELEMENT1 

6 

18 

1 

13.37418275 

0.0247108 

41.1904 

ELEMENT2 

349.3346 

66.5461 

296.1065 

ELEMENT3 

-0.00000029 

0.000000 

0.000091 

OUTPUT 

1 

2 

1 

940831.0 

074854.0 

86400.0 

ORB  TYPE 

14 

1 

8 

1.0 

OGOPT 

POTFIELD 

1 

7 

OUTOPT 

1 

940826074854.0 

940831074854.0 

450.0 

END 

FIN 


SGP4  EPHEMERIS  OUTPUT  SUMMARY: 


PAGE  12 


GTDS 

FINAL  REPORT 

SATELLITE  NAME 

EXPLORER 

SATELLITE  NUMBER 

65032 

RUN  REFERENCE  DATE 

AUG 

26, 

1994 

0 

HRS 

0 

MINS 

0.00000 

SECONDS 

RUN  EPOCH  DATE 

AUG 

26, 

1994 

7 

HRS 

48 

MINS 

53.51010 

SECONDS 

RUN  FINAL  TIME 

AUG 

31, 

1994 

7 

HRS 

48 

MINS 

54.00000 

SECONDS 

TOTAL  TIME  OF  FLIGHT 

5  DAYS 

0 

HRS 

0 

MINS 

0.48990 

SECONDS 

CAUSE  OF  TERMINATION 

SPECIFIED  TIME 

OF  FLIGHT 

REACHED 

***END  CONDITIONS*** 

CENTRAL  BODY  IS  EARTH  (INERTIAL  SYSTEM) 

X  0 . 48875493 06874350E+04  Y  -0 . 5498692794044931E+04 

MEAN  OF  1950.0 
Z 

--  EARTH  EQUATOR 

-0 . 176891328008 8411E+04 

VX 

0.4786071196Q55583E+01 

VY 

0.3047635044022132E+01 

vz 

0 . 4475570612493452E+01 

SMA 

0. 750190 8 902328932E+04 

'ecc 

0 . 24 9906230 13 83 419E-01 

INC 

0 . 4120299091118378E+02 

LAN 

0 . 32 797 707 05 618 6 8 3E+ 03 

AP 

0.90989096337 14 63 8E+02 

MA 

0 . 2511718748703 972E+03 

EA 

0 . 2498278467174672E+03 

P 

0 . 17 9 624 6 619 16 8 8 53E+ 01 

SLR 

0 . 7497223725870149E+04 

PR 

0 . 7 31443 1525 0667 02E+04 

APR 

0 . 7689386279591161E+04 

PH 

0 . 93 62965250 667 021E+03 

APH 

0 . 1311251279591161E+04 

C3 

-0 . 2656662492104218E+02 

TA 

0 . 24 84 8 943 77 6 14 734E+ 03 

RA 

0 .3 1224853 013 61 856E+ 03 

DEC 

-0 . 13 3 5 14 0772 8 1109 0E+02 

VPA 

0 . 913 442 013 7 012 03 2E+ 02 

AZ 

0 . 5064941371860270E+02 

RMAG 

0 . 756655898257545 0E+04 

VMAG 

0 . 7226706663664657E+01 

*** INITIAL  CONDITIONS*** 

CENTRAL  BODY  IS  EARTH  (INERTIAL  SYSTEM)  "NORAD"  TRUE  OF  REF.  --  EARTH  EQUATOR 


X 

0 . 72944688993 68215E+04 

Y 

-0 . 1366234653003144E+04 

Z 

0.6456847303414492E+01 

VX 

0.8537442543527309E+00 

VY 

0.547866115356 85 3 2E+01 

VZ 

0. 4850032676429656 E+01 

SMA 

0 . 749910 616 9 847256E+04 

ECC 

0. 2471 080000000 14 8E -01 

INC 

0 .4119040000000000E+02 

LAN 

0 . 3493346000000000E+03 

AP 

0 . 66 54 609 9999 9994 QE+0 2 

MA 

0 .29610650 00000 007E+ 03 

EA 

0 . 2 94 82146897 75155E+03 

P 

0 . 17 95240 08965 1203E+01 

SLR 

0 . 7494527038366274E+04 

PR 

0 . 7313797257105383E+04 

APR 

0 .76 84415 08 25 8912 8E+04 

PH 

0 .9356622571053831E+03 

APH 

0 . 13062 80082589128 E+ 04 

C3 

-0 .2657655399004164E+02 

TA 

0 . 2935295945886971E+03 

RA 

0 .34 9391562104 6252E+03 

DEC 

0 .4984967616062 8 8 IE- 01 

VPA 

0 . 912 852 07 94 12137 7E+0 2 

AZ 

0 . 4 8809624779 68625E+02 

RMAG 

0 . 74213 14 93 3 473 3 03E+ 04 

VMAG 

0 . 7366642671522530E+01 

SECTION  CENTRAL  BODY 
1  EARTH 


***  SECTIONING  SUMMARY'*** 


NUMBER  OF  SECTIONS  SCHEDULED  - 
NUMBER  OF  SECTIONS  COMPLETED  - 
TIME  OF  CROSSING 

H  M  S 

AUG  31  7  48  54.000 


1 

1 

TIME  INTO  FLIGHT 
H  M  S 
120  0  0.490 


CAUSE  OF  CROSSING 

SPECIFIED  TIME  OF  FLIGHT  REACHED 


***  FILE  GENERATION  SUMMARY  *** 


AN  ORB-1  FILE  HAS 
START  TIME  OF  THE 
END  TIME  OF  THE 


BEEN  GENERATED 

FILE  -  940826074854.000 

FILE  -  940831074654.000 


***  NUMBER  OF  GTDS  ERRORS  ENCOUNTERED-  0** 
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TRUTH  DATA  GENERATION 


COWELL  DC  &  EPHEM  DECK  (FIT  TO  SGP4-BASED  EPHEMERIS)  TO  GENERATE  TRUTH 


CONTROL 

DC 

EPOCH 

940826 . 0 

074853.5101 

ELEMENT 1 

1 

2 

1 

7502.98 

0.02566 

ELEMENT2 

349.0 

65.1 

OBS INPUT 

15 

940826074854.0 

940831074854. 

ORB TYPE 

2 

1 

1 

60.0 

DMOPT 

OBSDEV 

21 

22 

23 

500.0 

500.0 

OBSDEV 

24 

25 

26 

50.0 

50.0 

END 

* 

DCOPT 

PRINTOUT 

1 

4 

CONVERG 

30 

1 

l.D-3 

END 

OGOPT 

DRAG 

1 

1.0 

ATMOSDEN 

1 

SCPARAM 

l.D-6 

100. DO 

SOLRAD 

1 

1.0 

POTFIELD 

1 

2 

MAXDEGEQ 

1 

21.0 

MAXORDEQ 

1 

21.0 

STATEPAR 

1 

STATETAB 

1 

2 

3 

4.0 

5.0 

DRAGPAR 

1 

SOLRDPAR 

1 

END 

FIN 

CONTROL 

EPHEM 

OUTPUT 

OUTPUT 

1 

2 

1 

940903.0 

074854.0 

ORBTYPE 

2 

1 

1 

60.0 

OGOPT 

DRAG 

1 

1.0 

DRAGPAR 

0 

ATMOSDEN 

1 

SCPARAM 

l.D-6 

100. DO 

SOLRAD 

1 

1.0 

'  , 

SOLRDPAR 

0 

POTFIELD 

1 

2 

MAXDEGEQ 

1 

21.0 

MAXORDEQ 

1 

21.0 

OUTOPT 

22 

2 

1 

940826074854.0 

940903074854 

END 

FIN 

41.2 

297.22 


500. 0 
50.0 


EXPLORER 

86400.0 


DATA: 


65032A 


65032A 
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TRUTH  DATA  GENERATION 


COWELL  EPHEMERIS  OUTPUT  SUMMARY: 


ITERATION  NUMBER  5 


*********************************************** 


• 

POSITION  ERROR  RMS  (M) 

428.11032 

* 

SATELLITE 

EXPLORER  65032 

* 

CURRENT  WEIGHTED  RMS 

0.45208753 

* 

START  TIME  940826 

75624.000 

* 

PREDICTED  WEIGHTED  RMS 

0.45480147 

* 

EPOCH 

940826  74853.510 

* 

PREVIOUS  WEIGHTED  RMS 

0.45491167 

* 

END 

TIME 

940831 

74854.000 

* 

SMALLEST  WEIGHTED  RMS 

0.45491167 

* 

COORD.  SYSTEM 

MEAN  OF  1950.0 

* 

RELATIVE  CHANGE  IN  RMS 

6.20810804E-03 

* 

NO. 

OBS. 

AVAILABLE 

5760 

* 

PENALTY 

0.00000000 

* 

CENTRAL  BODY 

EARTH 

* 

* 

DC  HAS  CONVERGED 

* 

* 

NO. 

OBS. 

INCLUDED 

5760 

NO.  OBS.  ACCEPTED  5732  99  PCT 


OBSERVATION  SUMMARY  BY  TYPE 


TYPE 

X 

Y 

Z 

XDOT 

YDOT 

ZDOT 

TOTAL  NO. 

960 

960 

960 

960 

960 

960 

NO.  ACCEPTED 

949  (  98%) 

948  (  98%) 

955  (  99%) 

960  (100%) 

960  (100%) 

960  (100%) 

WEIGHTED  RMS 

0.5014 

0.4876 

0.4891 

0.4272 

0.3988 

0.3968 

MEAN  RESDUAL 

-69.59 

-129.6 

144.0 

-2.610 

0.4184 

2.070 

STANDARD  DEV 

240.9 

206.5 

197.7 

21.20 

19.94 

19.73 

OBSERVATION  SUMMARY 

BY  20  DEGREE 

TRUE  ANOMALY 

’  INTERVAL 

0 INTERVAL 

00-20 

20-40 

40-60 

60-80 

80-100 

100-120 

120-140 

140-160 

160-180 

TOTAL  NO. 

306 

306 

318 

336 

318 

324 

318 

318 

342 

NO.  ACCEPTED 

306  (100%) 

306  (100%) 

318  (100%) 

336  (100%) 

318  (100%) 

321  (  99%) 

317  (  99%) 

315  (  99%) 

342 

(100%) 

0 INTERVAL 

180-200 

200-220 

220-240 

240-260 

260-280 

280-300 

300-320 

320-340 

340-360 

TOTAL  NO. 

336 

330 

324 

318 

306 

312 

336 

306 

306 

NO.  ACCEPTED 

332  (  98%) 

327  (  99%) 

320  (  98%) 

317  (  99%) 

304  (  99%) 

308  (  98%) 

334  (  99%) 

305  (  99%) 

306 

(100%) 

THERE  ARE  B  SOLVE  PARAMETERS 

NUMBER  OP  DYNAMIC  SOLVE  PARAMETERS  -  8 

X  Y  2  VX  VY  VZ  RHOl  SOLRAD 

NUMBER  OP  STATION  LOCATIONS  BEING  SOLVED  FOR  -  0 

NUMBER  OF  BIASES  BEING  SOLVED  FOR  *  0 


***  FILE  GENERATION  SUMMARY  *** 
NO  FILES  HAVE  BEEN  GENERATED 

*  NUMBER  OF  GTDS  ERRORS  ENCOUNTERED  TO  DATE  -  5  * 
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HACC  OPTIMIZED  INPUT. DECK 


CONTROL  DC 
EPOCH 

ELEMENT1  361 
ELEMENT2 
OBSINPUT  15 
ORBTYPE  511 
DMOPT 

OBSDEV  21  22  23 
OBSDEV  24  25  26 
END 
DCOPT 

PRINTOUT  1  4 

CONVERG  30  1 

END 
OGOPT 

SOLRAD  1 
SPSRP  0 
SCPARAM 

SPGRVFRC  113 
SPZONALS  8  2  11 

SPMDAILY  882 
POTFIELD  1  2 

MAXDEGEQ  1 
MAXORDEQ  1 
RESONPRD 
STATEPAR  3 
STATETAB  123 
SOLRDPAR  1 
SSTESTFL  120 
SSTAPGFL  100 
END 
FIN 

CONTROL  EPHEM 
OUTPUT  121 
ORBTYPE  511 
OGOPT 

SOLRAD  1 
SPSRP  0 
SCPARAM 
.  SOLRDPAR  0 
SPGRVFRC  11  3 

SPZONALS  82  11 
SPMDAILY  882 
POTFIELD  1  2 

MAXDEGEQ  1 
MAXORDEQ  1 
RESONPRD 
OUTOPT  21 
END 
FIN 


940826.0 

7498.0 

349.3346 

940826074854.0 

43200.0 

100. 

10. 


l.D-4 


1.0 

l.D-6 

3.0 


21. 

21. 

86400.0 

4.0 

0.0 

0.0 


940831. 

43200.0 

1.0 

l.D-6 

3.0 


21. 

21. 

86400.0 

940826074854.0 


074853.5101 

0.0247108 

66.5461 

940828074854.0 

1.0 

100. 

10. 


100. DO 
3.0 


5.0 


0.0 


OUTPUT 

074854.0 

1.0 


100. DO 
3.0 


940831074854.0 


EXPLORER 

41.1904 

296.1065 


100. 

10. 


3.0 


1.0 


EXPLORER 

86400.0 


450.0 


65032A 


65032A 
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SUMMARY  OF  HACC  OPTIMIZED  DSST  FIT  TO  COWELL  TRUTH: 


ITERATION  NUMBER  4 


******************************************* 

•  *** 

* 

POSITION  ERROR  RMS  (M) 

48.705827  - 

* 

SATELLITE 

EXPLORER  65032 

* 

CURRENT  WEIGHTED  RMS 

0.28884820 

★ 

START  TIME  940826 

74954.000 

* 

PREDICTED  WEIGHTED  RMS 

0.30075942 

* 

EPOCH 

940826  74853.510 

* 

PREVIOUS  WEIGHTED  RMS 

0.30075980 

* 

END 

TIME 

940828 

74854.000 

* 

SMALLEST  WEIGHTED  RMS 

0.30075980 

* 

COORD.  SYSTEM 

MEAN  OF  1950.0 

* 

RELATIVE  CHANGE  IN  RMS 

3 . 96050355E-02 

* 

NO. 

OBS. 

AVAILABLE 

**** 

* 

PENALTY 

0.00000000 

* 

CENTRAL  BODY 

EARTH 

* 

DC  HAS  CONVERGED 

* 

* 

NO. 

OBS. 

INCLUDED 

**** 

******************************************* 

**** 

NO. 

OBS. 

ACCEPTED 

****  99  PCT 

OBSERVATION  SUMMARY  BY  TYPE 


TYPE  X 

TOTAL  NO.  2880 

NO.  ACCEPTED  2851  ( 
WEIGHTED  RMS  0.2873 
MEAN  RESDUAL  0.6894 
STANDARD  DEV  28.72 


Y 

2880 

98%)  2878  (  99%) 
0.3243 
0.3899 
32.42 


Z  XDOT  YDOT  ZDOT 

2880  2880  2B80  2880 

2880  (100%)  2815  (  97%)  2861  (  99%)  2880  (100%) 
0.2218  0.2951  0.3432  0.2431 

0.4194  8 . 5625E-02  -4.9562E-02  -6.2918E-02 

22.18  2.950  3.432  2.430 


OINTERVAL  00-20 

TOTAL  NO.  864 

NO.  ACCEPTED  864  (100%) 
(100%) 

OINTERVAL  180-200 

TOTAL  NO.  1026 

NO.  ACCEPTED  1026  (100%) 
99%) 


20-40 

954 

954  (100%) 

200-220 

948 

938  (  98%) 


OBSERVATION  SUMMARY  BY  20  DEGREE  TRUE  ANOMALY  INTERVAL 


40-60 

60-80 

80-100 

100-120 

120-140 

140-160 

160-180 

948 

1014 

972 

972 

972 

972 

1026 

942  ( 

99%) 

1005  ( 

99%) 

954  (  98%) 

963  (  99%) 

972  (100%) 

972  (100%) 

1026 

220-240 

240-260 

260-280 

280-300 

300-320 

320-340 

340-360 

948 

924 

888 

942 

1002 

936 

972 

938  ( 

98%) 

918  ( 

99%) 

683  (  99%) 

931  (  98%) 

987  (  98%) 

922  (  98%) 

970  ( 

THERE  ARE  7  SOLVE  PARAMETERS 


NUMBER  OF  DYNAMIC  SOLVE  PARAMETERS  -  7 

A  H  K  p  Q  LAM  SOLRAD 

NUMBER  OF  STATION  LOCATIONS  BEING  SOLVED  FOR  »  0 

NUMBER  OF  BIASES  BEING  SOLVED  FOR  -  0 


***  FILE  GENERATION  SUMMARY  *** 
NO  FILES  HAVE  BEEN  GENERATED 

*  NUMBER  OF  GTDS  ERRORS  ENCOUNTERED  TO  DATE  «  0  * 
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COMPARE  DECK  FOR  HACC: 


CONTROL  COMPARE 
COMPOPT 

CMPEPHEM  1102102  940828074854.0 

END 

FIN 


HACC  COMPARE  SUMMARY: 


EXPLORER  65032A 


940831074854.0  45. 


EPHEMERIS  COMPARISON  SUMMARY  REPORT 
Number  of  Points  Compared 


MINIMUM  POSITION  DIFFERENCE 


MAXIMUM  POSITION  DIFFERENCE 


YYMMDD 

HHMMSS.SSS 

Ocm) 

YYMMDD 

HHMMSS.SSS 

(km) 

RADIAL 

940829 

181854.000 

1.677306E-03 

940831 

23354.000 

7.173356E-02 

CROSS  TRACK 

940828 

220354.000 

1.377805E-04 

940831 

40354.000 

5.726527E-02 

ALONG  TRACK 

940828 

211854.000 

3.657224E-05 

940830 

134854.000 

2 . 025899E-01 

TOTAL 

940628 

220354.000 

4 . 484148E- 03 

940830 

134854.000 

2 . 03  3114E-01 

MINIMUM  VELOCITY 

DIFFERENCE 

MAXIMUM  VELOCITY 

DIFFERENCE 

YYMMDD 

HHMMSS.SSS 

{km/sec) 

YYMMDD 

HHMMSS.SSS 

(km/sec) 

RADIAL 

940828 

190354.000 

2 . 955573E-07 

940830 

121854.000 

1 . 702762E-04 

CROSS  TRACK 

940831 

23354.000 

2 . 772951E-08 

940831 

74854.000 

4 . 844473E-05 

ALONG  TRACK 

940828 

91854.000 

6 . 024580E-08 

940830 

181854.000 

6.031790E-05 

TOTAL 

940828 

211854.000 

1.714331E-05 

940830 

121854.000 

1 . 734637E-04 

RADIAL 
CROSS  TRACK 
ALONG  TRACK 
TOTAL 


POSITION  RMS 


2 . 8012E-02 
2 . 6375E-02 
7 . 19S1E-02 
8 . 1619E-02 

GTDS  COMPARE  PROGRAM 


VELOCITY  RMS 


7.0302E-05 
2 . 5504E-05 
2 . 8725E-05 
8 . 0112E-05 


NORMAL  COMPLETION  OF  JOB 
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GEO  Test  Case 
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TRUTH  DATA  GENERATION 


PROCEDURE: 

1)  Start  with  2-card  element  set 

2)  Propagate  2-card  element  set  forward  using  SGP4  to  obtain  SGP4-based  ephemeris 

3)  Fit  Cowell  theory  to  SGP4 -based  ephemeris 

4)  Propagate  Cowell  DC  solve-for  vector  forward  to  produce  truth  data 


CONTROL 

EPHEM 

EPOCH 

940825.0 

170334.9932 

ELEMENT1 

8 

18 

1 

1.00266891 

0.0003102 

ELEMENT2 

243.1923 

288.7716 

ELEMENT3 

0.00000000 

0.000000 

OUTPUT 

1 

2 

1 

940905.0 

170335.0 

ORBTYPE 

14 

1 

8 

1.0 

OGOPT 

POTFIELD 

1 

7 

OUTOPT 

1 

940825170335.0 

940905170335.0 

END 

FIN 


BS3A  20771 

0.0467 

167.5966 

0.00000 

86400.0 


450.0 


SGP4  EPHEMERIS  OUTPUT  SUMMARY: 
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GTDS  FINAL  REPORT 


SATELLITE  NAME 

BS3A 

SATELLITE  NUMBER 

20771 

RUN  REFERENCE  DATE 

AUG  25,  1994 

0  HRS 

0  MINS 

0.00000  SECONDS 

RUN  EPOCH  DATE 

AUG  25,  1994 

17  HRS 

3  MINS 

34.99320  SECONDS 

RUN  FINAL  TIME 

SEPT  5,  1994 

17  HRS 

3  MINS 

35.00000  SECONDS 

TOTAL  TIME  OF  FLIGHT 

11  DAYS 

0  HRS 

0  MINS 

0.00680  SECONDS 

CAUSE  OF  TERMINATION 

SPECIFIED  TIME 

OF  FLIGHT 

REACHED 

***END  CONDITIONS*** 

CENTRAL  BODY  IS  EARTH  (INERTIAL  SYSTEM) 


X 

0 ,4147498 5636 868 65E+05 

Y 

-0 . 767951853783 934 1E+ 04 

vx 

0 . 5 59624 64 03 20 947 5E+ 00 

VY 

0 . 3022273313587103E+01 

SMA 

0 . 4216 8764422 8 1726E+05 

ECC 

0. 2 74 613 16 86 64 405 7E -03 

LAN 

0 . 18743 913 81048706E+03 

AP 

0 . 3449220111341053E+03 

EA 

0 . 1777200 177 8165 12E+03 

P 

0 . 2393 8 3712 866472 6E+02 

PR 

0 . 421571S432480045E+05 

APR 

0 .4218034 4 5208 34 08E+05 

APH 

0 .358022 09520 834 08E+05 

C3 

-0.4726256572321093E+01 

RA 

0 . 3500817939879975E+03 

DEC 

0 . 62 9 19 9 1122 2 03 05 IE- 02 

AZ 

0 . 9002 013 0274 142 5 9E+02 

RMAG 

0 . 4218033535351635E+05 

CENTRAL  BODY  IS  EARTH  (INERTIAL 

X  0 . 3952629 83 87 00647E+05 

*** INITIAL  CONDITIONS*** 
SYSTEM) 

Y  -0 . 14 7 24 7 2 717 5 9 32 63E+0 5 

VX 

0 . 1073 17 24 9 65 6543 2E+01 

VY 

0. 28801818 07058570E+01 

SMA 

0 . 4216716110762471E+05 

ECC 

0 . 3 1019999 99597714E- 03 

LAN 

0 . 2431922999999999E+03 

AP 

0.2 8877159 99999954E+03 

EA 

0 .1676004163972050E+03 

P 

0 . 2 3937 00 6 044 266 92E+ 02 

PR 

0 . 4215408 085425082E+05 

APR 

0 . 42180241360 99 860E+05 

APH 

0.3580210 63 609985 9E+05 

C3 

-0 . 4726436278015460E+01 

RA 

0 . 33 95 6813431687 20E+03 

DEC 

0. 46411153393 00120E- 01 

AZ 

0. 9000518602416244 E+ 02 

RMAG 

0 .4217993 624 8902 80E+05 

MEAN  OF  1950.0  --  EARTH  EQUATOR 

Z  -0 .176776 690243 93 55E+03 
VZ  -0.3373151171135438E-02 


INC  0.2109068722704865E-01 
MA  0 . 1777193 918337379E+03 
SLR  0 . 4216876124276986E+05 
PH  0 .35779049324 80045E+05 
TA  0.1777206436437087E+03 


VPA  0. 8999937405206303 E+ 02 
VMAG  0 . 307365 04 53186 14 5E+01 


"NORAD"  TRUE  OF  REF.  --  EARTH  EQUATOR 
Z  0 . 34166 90188895451E+02 

VZ  -0 . 2780373401265742E-03 

INC  0 . 467000000000 00O0E-01 

MA  0 . 1675966000000055E+03 

SLR  0. 42 16715705 013 01 1E+ 05 

PH  0 . 35775945 8 542 508 1E+0 5 

TA  0 . 1676042322164656E+03 

VPA  0 . 8 9 996 1836026224 7E+02 

VMAG  0 . 3073621077231286E+01 


SECTION  CENTRAL  BODY 
1  EARTH 


***  SECTIONING  SUMMARY 

*** 

NUMBER 

OF  SECTIONS 

SCHEDULED  « 

1 

NUMBER 

OF  SECTIONS 

COMPLETED  - 

1 

TIME 

OF  CROSSING 

TIME 

INTO 

FLIGHT 

H  M 

S 

H 

M 

S 

SEPT 

5 

17  3 

35.000 

264 

0 

0.007 

CAUSE  OF  CROSSING 

SPECIFIED  TIME  OF  FLIGHT  REACHED 


***  FILE  GENERATION  SUMMARY  *** 

AN  ORB-1  FILE  HAS  BEEN  GENERATED 

START  TIME  OF  THE  FILE  «  940825170335.000 

END  TIME  OF  THE  FILE  -  940905170335.000 


***  NUMBER  OF  GTDS  ERRORS  ENCOUNTERED-  0** 
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TRUTH  DATA  GENERATION 


COWELL  DC  &  EPHEM  DECK  (FIT  TO  SGP4 -BASED  EPHEMERIS)  TO  GENERATE  TRUTH 


CONTROL 

DC 

EPOCH 

940825.0 

170334.9932 

ELEMENT1 

1 

2 

1 

42167.16 

0.000275 

ELEMENT2 

97.6 

75.3 

OBS INPUT 

15 

940825170335.0 

940904170335.0 

ORETYPE 

2 

1 

1 

60.0 

DMO  PT 

OBSDEV 

21 

22 

23 

500.0 

500.0 

OBSDEV 

24 

25 

26 

50.0 

50.0 

END 

DCOPT 

PRINTOUT 

1 

4 

CONVERG 

30 

1 

l.D-3 

END 

OGOPT 

SOLRAD 

1 

1.0 

SCPARAM 

l.D-6 

100. DO 

POTFIELD 

1 

2 

MAXDEGEQ 

1 

21.0 

MAXORDEQ 

1 

21.0 

STATEPAR 

1 

STATETAB 

1 

2 

3 

4.0 

5.0 

SOLRDPAR 

1 

END 

FIN 

CONTROL 

EPHEM 

OUTPUT 

OUTPUT 

1 

2 

1 

940905.0 

170335.0 

ORB TYPE 

2 

1 

1 

60.0 

OGOPT 

SOLRAD 

1 

1.0 

SCPARAM 

l.D-6 

100. DO 

SOLRDPAR 

0 

POTFIELD 

1 

2 

MAXDEGEQ 

1 

21.0 

MAXORDEQ 

1 

21.0 

OUTOPT 

21 

940825170335.0 

940905170335.0 

END 

FIN 

0.23 

166.09 


500.0 

50.0 


BS3A 

86400.0 


DATA: 


20771 


20771 
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TROTH  DATA  GENERATION 


COWELL  EPHEMERIS  OUTPUT  SUMMARY: 


ITERATION  NUMBER  3 


*********************************************** 


* 

POSITION  ERROR  RMS  (M) 

2453.3437 

* 

SATELLITE 

BS3A  20771 

* 

CURRENT  WEIGHTED  RMS 

1 .9976520 

* 

START  TIME  940825 

171105.000 

* 

PREDICTED  WEIGHTED  RMS 

2.5423798 

* 

EPOCH 

940825  170334.993 

* 

PREVIOUS  WEIGHTED  RMS 

2.5424042 

* 

END 

TIME 

940904 

170335.000 

* 

SMALLEST  WEIGHTED  RMS 

2.5424042 

* 

COORD.  SYSTEM 

MEAN  OF  1950.0 

* 

RELATIVE  CHANGE  IN  RMS 

0.21426264 

* 

NO. 

OBS. 

AVAILABLE 

**** 

* 

PENALTY 

0.00000000 

* 

CENTRAL  BODY 

EARTH 

* 

* 

DC  HAS  CONVERGED 

* 

* 

NO. 

OBS. 

INCLUDED 

**** 

***********************************************  NO.  OBS .  ACCEPTED  ****  97  PCT 


OBSERVATION  SUMMARY  BY  TYPE 


TYPE 

X 

Y 

Z 

XDOT 

YDOT 

ZDOT 

TOTAL  NO. 

1920 

1920 

1920 

1920 

1920 

1920 

NO.  ACCEPTED 

1876  (  97%) 

1920  (100%) 

1698  ( 

88%) 

1920  (100%) 

1920  (100%) 

1920  (100%) 

WEIGHTED  RMS 

2.180 

2.191 

3.943 

0.1615 

0.1515 

0.3833 

MEAN  RESDUAL 

-113.6 

-2.593 

-102.3 

0.1289 

-0.1550 

-0.1233 

STANDARD  DEV 

1084. 

1095. 

1969. 

8.075 

7.574 

19.17 

OBSERVATION  SUMMARY  BY  20  DEGREE 

TRUE  ANOMALY 

’  INTERVAL 

0 INTERVAL 

00-20 

20-40 

40-60 

60-80 

80-100 

100-120 

120-140 

140-160 

160-180 

TOTAL  NO. 

690 

684 

684 

678 

672 

660 

630 

594 

564 

NO.  ACCEPTED 

690  (100%) 

682  (  99%) 

670  ( 

97%) 

649  (  95%) 

633  (  94%) 

621  (  94%) 

599  {  95%) 

594  (100%) 

564 

(100%) 

0  INTERVAL 

180-200 

200-220 

220-240 

240-260 

260-280 

280-300 

300-320 

320-340 

340-360 

TOTAL  NO. 

582 

588 

612 

612 

630 

636 

648 

672 

684 

NO.  ACCEPTED 

582  (100%) 

584  (  99%) 

602  ( 

98%) 

596  {  97%) 

602  (  95%) 

607  (  95%) 

629  (  97%) 

666  (  99%) 

684 

{100%) 

THERE  ARE 

7  SOLVE  PARAMETERS 

NUMBER  OF  DYNAMIC  SOLVE 

PARAMETERS  - 

7 

X 

Y 

Z 

VX 

VY 

VZ 

SOLRAD 

NUMBER  OP  STATION  LOCATIONS  BEING  SOLVED  FOR  -  0 


NUMBER  OP  BIASES  BEING  SOLVED  FOR  -  0 


***  FILE  GENERATION  SUMMARY  *** 
NO  FILES  HAVE  BEEN  GENERATED 

*  NUMBER  OF  GTDS  ERRORS  ENCOUNTERED  TO  DATE  =  0  * 
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GEOSYNCHRONOUS  OPTIMIZED  INPUT  DECK 


CONTROL 

DC 

EPOCH 

940825.0 

ELEMENT1 

1 

6 

1 

42167.16 

ELEMENT2 

97.6 

OBS INPUT 

15 

940825170335.0 

ORB TYPE 

5 

1 

1 

43200.0 

DMOPT 

OBSDBV 

21 

22 

23 

100. 

OBSDEV 

24 

25 

26 

10. 

END 

DCOPT 

PRINTOUT 

1 

4 

CONVERG 

30 

1 

l.D-4 

BS3A  20771 

170334.9932 

0.000275  0.23 

75.3  166.09 

940828170335.0 

1.0 

100.  100- 

10.  10. 


OGOPT 


SOLRAD 

1 

1.0 

SCPARAM 

l.D-6 

100. DO 

SPGRVFRC 

1 

3 

3 

2.0 

1.0 

3.0 

SPNUMGRV 

7 

1 

6 

2.0 

2.0 

3600.0 

POTFIELD 

1 

2 

MAXDEGEQ 

1 

4. 

MAXORDEQ 

1 

4. 

STATEPAR 

3 

STATETAB 

1 

2 

3 

4.0 

5.0 

6.0 

SOLRDPAR 

1 

SSTESTFL 

1 

2 

0 

o 

o 

SSTAPGFL 

1 

0 

0 

0.0 

0.0 

1.0 

END 

FIN 

CONTROL 

EPHEM 

OUTPUT 

BS3A 

OUTPUT 

1 

2 

1 

940905. 

170335.0 

86400.0 

ORB TYPE 

5 

1 

1 

43200.0 

1.0 

OGOPT 

SOLRAD 

1 

1.0 

SOLRDPAR 

0 

SCPARAM 

l.D-6 

100. DO 

SPGRVFRC 

1 

3 

3 

2.0 

1.0 

3.0 

SPNUMGRV 

7 

1 

8 

2.0 

2.0 

3600.0 

POTFIELD 

1 

2 

MAXDEGEQ 

1 

4. 

MAXORDEQ 

1 

4. 

OUTOPT. 

21 

940825170335.0 

940905170335.0 

450.0 

END 

FIN 


20771 
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SUMMARY  OF  GEOSYNCHRONOUS  OPTIMIZED  DSST  FIT  TO  COWELL  TRUTH: 


ITERATION  NUMBER  4 


SATELLITE 

EPOCH 

COORD.  SYSTEM 
CENTRAL  BODY 


BS3A  20771 

940825  170334.993 
MEAN  OF  1950.0 
EARTH 


********************************************* 

** 

* 

POSITION  ERROR  RMS  (M) 

5.1035665 

* 

* 

CURRENT  WEIGHTED  RMS 

2.20243188E-02 

* 

START  TIME  940825 

171105.000 

* 

PREDICTED  WEIGHTED  RMS 

2 . 20766417E-02 

* 

* 

PREVIOUS  WEIGHTED  RMS 

2 . 20767451E-02 

* 

END 

TIME 

940828 

170335.000 

* 

SMALLEST  WEIGHTED  RMS 

2 . 2076745 IE- 02 

* 

* 

RELATIVE  CHANGE  IN  RMS 

2.37472948E-03 

* 

NO. 

OBS. 

AVAILABLE 

3456 

* 

PENALTY 

0.00000000 

* 

* 

DC  HAS  CONVERGED 

* 

* 

NO. 

OBS. 

INCLUDED 

3456 

*********************************************** 

NO. 

OBS. 

ACCEPTED 

3454  99  PCT 

OBSERVATION  SUMMARY  BY  TYPE 


TYPE  X 

Y 

TOTAL  NO.  576 

576 

NO.  ACCEPTED  576  (100%) 

574  {  99%) 

WEIGHTED  RMS  3.6111E-02 

3 . 6011E-02 

MEAN  RESDUAL-4 . 2038E-02 

4 . 1731E-02 

STANDARD  DEV  3.611 

3.601 

0 INTERVAL 

00-20 

20-40 

TOTAL  NO. 

210 

198 

NO.  ACCEPTED 
98%) 

210  (100%) 

198  (100%) 

0 INTERVAL 

180-200 

200-220 

TOTAL  NO. 

174 

186 

NO.  ACCEPTED 
(100%) 

174  (100%) 

186  (100%) 

Z  XDOT  YDOT 

576  576  576 

576  (100%)  576  (100%)  576  (100%) 

8 . 9331E-04  1.2533E-02  1.2431E-02 

-4 . 7060E-02  2 . 9880E-03  1.7902E-03 

7.5930E-02  0.1253  0.1243 

OBSERVATION  SUMMARY  BY  20  DEGREE 
40-60  60-80  80-100 

186  180  180 

186  (100%)  180  (100%)  180  (100%) 

220-240  240-260  260-280 

204  210  216 

204  (100%)  210  (100%)  216  (100%) 


ZDOT 

576 

576  (100%) 

6 . 6651E-05 
6 . 5442E-05 
6 . 6329E-04 

TRUE  ANOMALY  INTERVAL 
100-120  120-140  140-160 

180  180  180 

180  (100%)  180  (100%)  180  (100%) 

280-300  300-320  320-340 

204  198  204 

204  (100%)  198  (100%)  204  (100%) 


160-180 
168 
166  ( 

340-360 

198 

198 


THERE  ARE  7  SOLVE  PARAMETERS 


NUMBER  OF  DYNAMIC  SOLVE  PARAMETERS  =  7 

A  H  K  P  Q  LAM  SOLRAD 

NUMBER  OF  STATION  LOCATIONS  BEING  SOLVED  FOR  =■=  0 

:  NUMBER  OF  BIASES  BEING  SOLVED  FOR  =  0  ... 


***  FILE  GENERATION  SUMMARY  *** 
NO  FILES  HAVE  BEEN  GENERATED 

*  NUMBER  OF  GTDS  ERRORS  ENCOUNTERED  TO  DATE  *  0  * 
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COMPARE  DECK  FOR  GEOSYNCHRONOUS  CASE: 


CONTROL  COMPARE 
COMPOPT 

CMPEPHEM  1102102  940828170335.0  940831170335.0 

END 

FIN 


20771 


GEOSYNCHRONOUS  COMPARE  SUMMARY: 


EPHEMERIS  COMPARISON  SUMMARY  REPORT 
Number  of  Points  Compared  577 


MINIMUM  POSITION  DIFFERENCE 

MAXIMUM  POSITION  DIFFERENCE 

YYMMDD 

HHMMSS.SSS 

(km) 

YYMMDD 

HHMMSS.SSS 

(km) 

RADIAL 

CROSS  TRACK 
ALONG  TRACK 

TOTAL 

940829 

940829 

940828 

940828 

111835.000 

51835.000 

170335.000 

170335.000 

4 . 641092E-03 
7.078641E-08 

1 . 198785E-03 

5 . 880411E-03 

940831 

940829 

940831 

940831 

5605.000 

222605.000 

170335.000 

170335.000 

8 . 159907E-03 
3 . 100871E-04 
3.501914E-02 
3.581756E-02 

MINIMUM  VELOCITY  DIFFERENCE 

MAXIMUM  VELOCITY  DIFFERENCE 

YYMMDD 

HHMMSS.SSS 

(km/sec) 

YYMMDD 

HHMMSS.SSS 

(km/sec) 

RADIAL 

940828 

170335.000 

1.860573E-06 

940831 

170335.000 

4.400656E-06 

CROSS  TRACK 

940831 

134835.000 

1 . 779926E-11 

940830 

181105.000 

1 . 410646E-08 

ALONG  TRACK 

940831 

115605.000 

2 . 699794E-07 

940831 

10335.000 

5 . 553332E-07 

TOTAL 

940828 

170335.000 

1. 901702E-06 

940831 

170335.000 

4.426126E-06 

POSITION  RMS 

VELOCITY  R t 

(km) 

(km/sec) 

RADIAL 

6 . 2625E-03 

3 . 0904E-06 

CROSS  TRACK 

1 . 6640E-04 

9.0059E-09 

ALONG  TRACK 

2 . 0304E-02 

4.1029E-07 

TOTAL 

2.1249E-02 

GTDS  COMPARE  PROGRAM 

3 . 1175E-06 

NORMAL  COMPLETION  OF  JOB 


16 
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MOLY  Test  Case 
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TRUTH  DATA  GENERATION 


PROCEDURE: 

1)  Start  with  mean  Keplerian  elements 

2)  Propagate  forward  one  day  using  DSST  to  obtain  osculating  elements 

3)  Propagate  osculating  elements  forward  using  Cowell  theory  to  generate  truth  data 


CONTROL 

EPHEM 

NSSC 

EPOCH 

790803.0 

234212.0 

ELEMENT1 

1 

6 

1 

26556.9582 

0.6990986 

63.173001 

ELEMENT2 

190.619681 

281.59624 

13.29315 

OUTPUT 

1 

2 

1 

790804.0 

234212.0 

43200.0 

ORB TYPE 

5 

1 

1 

43200.0 

1.0 

OGOPT 

SPSHPER 

5 

DRAG 

1 

1.0 

ATMOS DEN 

1 

SCPARAM 

12.5D-6 

1250. DO 

SOLRAD 

1 

1.0 

POTFIELD 

1 

5 

MAXDEGEQ 

1 

8. 

MAXORDEQ 

1 

8. 

END 

FIN 

CONTROL 

EPHEM 

OUTPUT 

NSSC 

OUTPUT 

1 

2 

1 

790818.0 

234212.0 

172800.0 

ORBTYPE 

1 

1 

1 

200.0 

1.5 

OGOPT 

DRAG 

1 

1.0 

ATMOSDEN 

1 

SCPARAM 

12.5D-6 

1250. DO 

SOLRAD 

1 

1.0 

POTFIELD 

1 

6 

MAXDEGEQ 

1 

21. 

MAXORDEQ 

1 

21. 

OUTOPT 

1 

790804234212.0 

790818234212.0 

450. 

END 

FIN 

9829 


9829 
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TRUTH  DATA  GENERATION 


COWELL  EPHEMERIS  OUTPUT  SUMMARY: 


PAGE  26 


GTDS  FINAL  REPORT 


SATELLITE  NAME 
SATELLITE  NUMBER 
RUN  REFERENCE  DATE 
RUN  EPOCH  DATE 
RUN  FINAL  TIME 
TOTAL  TIME  OF  FLIGHT 
CAUSE  OF  TERMINATION 


NSSC 

9829 

AUG 

3,  1979 

0  HRS 

0  MINS 

0.00000 

SECONDS 

AUG 

4,  1979 

23  HRS 

42  MINS 

12.00000 

SECONDS 

AUG 

18,  1979 

23  HRS 

42  MINS 

12.00000 

SECONDS 

14  DAYS 

0  HRS 

0  MINS 

0.00000 

SECONDS 

SPECIFIED  TIME  OF  FLIGHT  REACHED 


CENTRAL  BODY  IS  EARTH  (INERTIAL  SYSTEM) 

X  -0 . 1371825581792548E+05  Y 

VX  0.8606074345775577E+00  VY 


*  *  * END  CONDITIONS  *  * * 

-0 . Ill 9 942 65 2 170 OS 0E+ 05 
-0.1717010361617677E+01 


MEAN  OF  1950.0  --  EARTH  EQUATOR 

Z  0 . 1774 66569 7 1454 02E+05 

VZ  0 .  3619315325579901E+01 


SMA  0 . 26556804798674 68E+05 
LAN  0 . 1888106505833530E+03 
EA  0 . 8 540744767 50 662 1E+02 
PR  0.8001 94 8767031005E+04 
APH  0. 387335228303183 6E+05 


ECC  0 . 6986855599650172E+00 
AP  0 . 2815323 01483 0572E+0 3 
P  0.119638S453041661E+02 
APR  0 . 451116608303 183 6E+05 
C3  -0.7504676165332438E+01 


INC  0 . 6319595630208168E+02 
MA  0 . 4550424397783793E+02 
SLR  0 . 1359279482213544E+05 
PH  0 . 1623810767031005E+04 
TA  0 . 13 094046115 44 9 84E+03 


RA  0.2192278053194159E+03 
AZ  0 .  3 96726163 9170771E+02 


DEC  0 .  450 6044124614 617E+02 
RMAG  0.2507112942736071E+05 


VPA  0 . 4577045721988872E+02 
VMAG  0.4097342207374211E+01 


CENTRAL  BODY  IS  EARTH  (INERTIAL 
X  -0 . 12429 83594562453E+05 
VX  -0 . 248 8 03453042124 7E+01 

SMA  0 . 2 657 2 17 5 9 14 71 07 0E+ 05 
LAN  0 . 1904 8233 87631476E+03 
EA  0.4248 63 46205003 5 1E+02 
PR  0.7993192062777778E+04 
APH  0 .3877302176664363E+05 

RA  0.1936893756438893E+03 
AZ  0 . 269994384960544BE+02 


*** INITIAL  CONDITIONS*** 
SYSTEM) 

Y  -0 .3 027 62566273 8 115E+ 04 
VY  -0 .3 27707117 66 95446E+01 

ECC  0 . 699 1894 044193 557E+ 00 
AP  0 . 2816075013505303E+03 
P  0 . 119742 7 3 103 119 3 0E+ 02 

APR  0.4515115976664363E+05 
C3  -0 . 75 00334 960889100E+01 

DEC  0 . 6313668242350110E+01 
RMAG  0 . 12 8713 2 143 3 4 54 07E+0 5 


1950.0 

--  EARTH  EQUATOR 

Z 

0 . 1415477535 8457 10E+04 

VZ 

0 . 54777759852 67707E+01 

INC 

0.63 17 748516631303 E+ 02 

MA 

0 . 1542 8 83 506600 653E+02 

SLR 

0 . 135 8 194726056 090E+05 

PH 

0 . 1615 054 062777 7 7 8E+04 

TA 

0 .8547104510736438E+02 

VPA 

0 . 5 6553 7 2 54 5 8632 63E+ 02 

VMAG 

0 . 685 095183 65313 40E+01 

***  SECTIONING  SUMMARY 

*  +  * 

NUMBER 

OF  SECTIONS 

SCHEDULED  » 

1 

NUMBER 

OF  SECTIONS 

COMPLETED  - 

1 

SECTION 

CENTRAL  BODY 

TIME  OF  CROSSING 

TIME 

INTO 

FLIGHT 

H  M 

S 

H 

M 

S 

1 

EARTH 

AUG  18 

23  42 

12.000 

336 

0 

0.000 

CAUSE  OF  CROSSING 

SPECIFIED  TIME  OF  FLIGHT  REACHED 


***  FILE  GENERATION  SUMMARY  *** 


AN  ORB-1  FILE  HAS  BEEN  GENERATED 

START  TIME  OF  THE  FILE  «  790804234212.000 

END  TIME  OF  THE  FILE  *  790818234212.000 


***  NUMBER  OF  GTDS  ERRORS  ENCOUNTERED  *  1** 
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molniya  optimized  input  deck 


CONTROL 

OGOPT 

POTFIELD 

END 

FIN 

CONTROL 

EPOCH 

DATAMGT 

1  6 

DC 

790804.0 

ELEMENT1  1 
ELEMENT2 

OBS INPUT  15 

6 

1 

26572.176 

190.5 

790804234212.0 

ORB TYPE 

DMOPT 

5 

1 

1 

43200.0 

OBSDEV 

21 

22 

23 

100. 

OBSDEV 

END 

24 

25 

26 

10. 

NSSC 

234212.0 

0.699  63.2 

281.6  15.429 

790807234212.0 

1.0 


100.  100. 

10.  10. 


OGOPT 

DRAG  1 

ATMOSDEN  1 

SPDRAG  0 
SCPARAM 
SOLRAD  1 
SPSRP  0 
SPGRVFRC  113 
SPZONALS  430 
SPMDAILY  442 
MAXDEGEQ  1 
MAXORDEQ  1 
STATEPAR  3 
STATETAB  123 
DRAG  PAR  1 
DRAGPAR2  1  1 

SOLRDPAR  1 
SSTESTFL  120 
SSTAPGFL  100 
END 
DCOPT 

PRINTOUT  1  4 

CONVERG  30  1 

END 
FIN 

CONTROL  EPHEM 
OUTPUT  121 
ORB TYPE.  5  11 

OGOPT 

DRAGPAR  0 
DRAG  1 

ATMOSDEN  1 

SPDRAG  0 
SCPARAM 
SOLRAD  1 
SPSRP  0 
SOLRDPAR  0 
SPGRVFRC  113 
SPZONALS  438 
SPMDAILY  442 
MAXDEGEQ  1 
MAXORDEQ  1 
OUTOPT  21 
END 
FIN 


1.0 


12.5D-6 

1.0 


0.0 

1.0 


l.D-3 


790818.0 

43200.0 


12.5D-6 

1.0 


8. 

8. 

790804234212.0 


1250. DO 


5.0 


OUTPUT 

234212.0 

1.0. 


1250. DO 


790818234212.0 


1.0 


NSSC 

172800.0 


450. 


9829 


9829 


9829 
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SUMMARY  OF  MOLNIYA  OPTIMIZED  DSST  FIT  TO  COWELL  TRUTH; 


ITERATION  NUMBER  4 


SATELLITE 

EPOCH 

COORD.  SYSTEM 
CENTRAL  BODY 


NSSC  9829 

790804  234212.000 
MEAN  OF  1950.0 
EARTH 


*********************************************** 

*  POSITION  ERROR  RMS  (M)  114.82174  * 

*  CURRENT  WEIGHTED  RMS  0.45742789  * 

*  PREDICTED  WEIGHTED  RMS  0.51345696  * 

*  PREVIOUS  WEIGHTED  RMS  0.51349399  * 

*  SMALLEST  WEIGHTED  RMS  0.51349399  * 

*  RELATIVE  CHANGE  IN  RMS  0.10918550  * 

*  PENALTY  0.00000000  * 

*  DC  HAS  CONVERGED  * 

*  * 
*********************************************** 

OBSERVATION  SUMMARY  BY  TYPE 


START  TIME  790804  234212.000 
END  TIME  790807  233442.000 
NO.  OBS.  AVAILABLE  3456 
NO.  OBS.  INCLUDED  3456 
NO.  OBS.  ACCEPTED  3382  97  PCT 


TYPE 

X 

Y 

Z 

XDOT 

YDOT 

ZDOT 

TOTAL  NO. 

576 

576 

576 

576 

576 

576 

NO.  ACCEPTED 

502  (  87%) 

576  (100%) 

576  (100%) 

576  (100%) 

576  (100%) 

576  (100%) 

WEIGHTED  RMS 

0.9218 

0.4468 

0.5014 

0.1340 

9 . 2191E-02 

0.1024 

MEAN  RESDUAL 

-18.03 

-5.572 

19.20 

-0.2010 

■1 . 0291E-02  - 

•3 . 9840E-02 

STANDARD  DEV 

90.40 

44.34 

46.32 

1.325 

0.9219 

1.023 

OBSERVATION  SUMMARY  BY  20  DEGREE 

TRUE  ANOMALY 

0 INTERVAL 

00-20 

20-40 

40-60 

60-80 

80-100 

100-120  : 

TOTAL  NO. 

24 

24 

30 

48 

84 

120 

NO.  ACCEPTED 

24  (100%) 

24  (100%) 

30  (100%) 

48  (100%) 

84  (100%) 

120  (100%) 

99%) 

0 INTERVAL 

180-200 

200-220 

220-240 

240-260 

260-280 

280-300 

TOTAL  NO. 

714 

450 

234 

120 

72 

48 

NO.  ACCEPTED 

669  (  93%) 

425  (  94%) 

234  (100%) 

120  (100%) 

72  (100%) 

48  (100%) 

(100%) 

THERE  ARE 

8  SOLVE  PARAMETERS 

NUMBER  OF  DYNAMIC  SOLVE 

PARAMETERS  * 

8 

A 

H 

K 

P 

Q 

LAM 

CSUBD 

120-140 

240 


140-160 

450 


160-180 

700 


300-320  320-340  340-360 

36  24  30 

36  (100%)  24  (100%)  30 


NUMBER  OF  STATION  LOCATIONS  BEING  SOLVED  FOR 

\ 

NUMBER  OF  BIASES  BEING'  SOLVED  FOR  •  0 


***  FILE  GENERATION  SUMMARY  *** 


NO  FILES  HAVE  BEEN  GENERATED 

*  NUMBER  OF  GTDS  ERRORS  ENCOUNTERED  TO  DATE  -  0  * 
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COMPARE  DECK  FOR  MOLNIYA  CASE: 


CONTROL  COMPARE 
COMPOPT 

CMPEPHEM  1102102  790807234212.0  790810234212.0 

END 

FIN 


9829 


MOLNIYA  COMPARE  SUMMARY: 


EPHEMERIS  COMPARISON  SUMMARY  REPORT 
Number  of  Points  Compared  577 


MINIMUM  POSITION  DIFFERENCE 


MAXIMUM  POSITION  DIFFERENCE 


YYMMDD 

HHMMSS.SSS 

(km) 

YYMMDD 

HHMMSS.SSS 

(km) 

RADIAL 

CROSS  TRACK 
ALONG  TRACK 
TOTAL 

790810 

790808 

790810 

790808 

104942.000 

82712.000 

201212.000 

121942.000 

4 . 686057E-04 
5.364643E-06 

1 . 358405E-03 

3 . 791407E-02 

790810 

790810 

790810 

790810 

220442.000 

45712.000 

224212.000 

224212.000 

3 . 501438E-01 
5.219956E-02 
7.213969E-01 
7.334274E-01 

MINIMUM  VELOCITY  DIFFERENCE 

MAXIMUM  VELOCITY 

DIFFERENCE 

YYMMDD 

HHMMSS.SSS 

(km/sec) 

YYMMDD 

HHMMSS.SSS 

(km/sec) 

RADIAL 

CROSS  TRACK 
ALONG  TRACK 
TOTAL 

790808 

790808 

790810 

790808 

200442.000 

44942.000 

81212.000 

131942.000 

1 . 490001E-09 

1 . 377708E-09 
1.646474E-08 

1 . 893925E-06 

790810 

790810 

790808 

790810 

224942.000 

224212.000 

224942.000 

224942.000 

4 . 816413E-04 
1.892690E-05 
5 . 052250E-05 
4 . 820318E-04 

POSITION  RMS 

VELOCITY  R> 

(km) 

(km/sec) 

RADIAL 

1 . 2727E-01 

7 . 3599E-05 

CROSS  TRACK 

2 . 2510E-02 

‘  3 . 4218E-06 

ALONG  TRACK 

1.7720E-01 

1.3041E-05 

TOTAL 

2.1933E-01 

GTDS  COMPARE  PROGRAM 

7 . 4824E-05 

PAGE 


0 


NORMAL  COMPLETION  OF  JOB 


16 
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NEEC  Test  Case 
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TRUTH  DATA  GENERATION 


PROCEDURE: 

X)  Start  with  2-card  element  set 

2)  Propagate  2-card  element  set  forward  using  SGP4  to  obtain  SGP4-based  ephemeris 

3)  Pit  Cowell  theory  to  SGP4 -based  ephemeris 

4)  Propagate  Cowell  DC  solve- for  vector  forward  to  produce  truth  data 


CONTROL 

EPHEM 

VANGARD2 

59001A 

EPOCH 

940826.0 

073513.6735 

ELEMENT1 

8 

18 

1 

11.73921485 

0.1522640 

32.8834 

ELEMENT2 

251.8592 

10.8368 

352.1515 

ELEMENT3 

-0.00000000 

0.000000 

0.000154 

OUTPUT 

1 

2 

1 

940905.0 

073514.0 

86400.0 

OEBTYPE 

14 

1 

8 

1.0 

OGOPT 

POTFIELD 

1 

7 

OUTOPT 

1 

940826073514.0 

940905073514.0 

450.0 

END 

PIN 


SGP4  EPHEMERIS  OUTPUT  SUMMARY: 

GTDS  FINAL  REPORT 

PAGE  13 


SATELLITE  NAME 

SATELLITE  NUMBER 

RUN  REFERENCE  DATE 

VANGARD2 

59001 

AUG  26,  1994 

0 

HRS 

0 

MINS 

0.00000 

SECONDS 

RUN  EPOCH  DATE 

AUG  26,  1994 

7 

HRS 

35 

MINS 

13.67350 

SECONDS 

RUN  FINAL  TIME 

SEPT  5,  1994 

7 

HRS 

35 

MINS 

14.00000 

SECONDS 

TOTAL  TIME  OF  FLIGHT 

10  DAYS 

0 

HRS 

0 

MINS 

0.32650 

SECONDS 

CAUSE  OF  TERMINATION 

SPECIFIED  TIME 

OF  FLIGHT 

REACHED 

*« 

*END  CONDITIONS*** 

CENTRAL  BODY  IS  EARTH  {INERTIAL  SYSTEM) 

MEAN  OF  1950.0 

--  EARTH  EQUATOR 

X 

0 . 42269055303 87705E+04 

Y 

0. 766 84 016 6744583 8E+ 04 

Z 

-0 . 25512729 943 07648E+04 

vx 

-0 .4 861278344 2 13523E+01 

VY 

0 . 2411752650711627E+01 

VZ 

-0. 30400895710 74028 E+01 

SMA 

0 . 8181674036773987E+04 

•  ECC 

0 . 1520921345305507E+00 

INC 

0 .32 889 83 668642 64 6E+02 

LAN 

0 . 21508 6274 6053 12 0E+03 

AP 

0 . 663614 0824 14 7962E+02 

MA 

0 .1332440193797 609E+03 

EA 

0 . 1389650806767234E+03 

P 

0 .204583942699 5563E+01 

SLR 

0 . 7992415410708262E+04 

PR 

0 . 6937305768487B44E+04 

APR 

0 . 94260423 05 06012 9E+04 

PH 

0 . 559 17  076  84  8  7  843  8E+03 

APH 

0 . 3 0479 073 05 060 12 9E+ 04 

C3 

-0 . 2435936693446952E+02 

TA 

0 .1444016201663212E+03 

RA 

0 . 61643 8 63 2 2 942 16 7E+ 02 

DEC 

-0 . 16 12 6 04 0 8 6662 16 8E+0 2 

VPA 

0 . 8423117145710 005E+02 

AZ 

0 . 119 0 6 02 0 65 2 9 04 3 1E+ 03 

RMAG 

0 . 9120312967737114E+04 

VMAG 

0 . 622 018670722573 8E+01 

CENTRAL  BODY  IS  EARTH 


*** INITIAL 
{INERTIAL  SYSTEM) 


CONDITIONS*** 


"NORAD"  TRUE  OF  REF.  --  EARTH  EQUATOR 


X 

-0 . 2158590909184290E+04 

Y 

-0 . 6 6079170 044 6276 9B+04 

Z 

0 .3943505275263057E+01 

VX 

0 .6544250260129505E+01 

VY 

-0 . 1923 6275518 5 9 584E+01 

VZ 

0.4407890047245081E+01 

SMA 

0 . 81810 2 7 89 17 666 8 7E+04 

ECC 

0. 152263999999999 8E+00 

INC 

0 . 3288339999999999E+02 

LAN 

0.2518591999999999E+03 

AP 

0. 108368000000000 0E +02 

MA 

0 . 3521514999999999E+03 

EA 

0 . 3 5074 9 02 115 9513 9E+ 03 

P 

0 . 204559707754 0084E+01 

SLR 

0 . 7991356276595909E+04 

PR 

0 . 6935351860854726E+04 

APR 

0 . 942 67 03 92267 8 64 9E+04 

PH 

0 . 55 72168 608547 254E+03 

APH 

0 . 3048568922678649E+04 

C3 

-0 . 2436129085937650E+02 

TA 

0 . 3492230657582 670E+03 

RA 

0 . 2 5 19 0 9473 9 04 2 65 6E+0 3 

DEC 

0.3250298192318978E-01 

VPA 

0 .9141873484191045E+02 

AZ 

0.5711661425979332E+02 

RMAG 

0 . 69515535963 06017E+04 

VMAG 

0 . 81213 94 52 8 9 02 080E+ 01 

***  SECTIONING  SUMMARY  *** 

NUMBER  OF 

SECTIONS 

SCHEDULED  -  1 

NUMBER  OF 

SECTIONS 

COMPLETED  -  1 

CENTRAL  BODY  TIME  OF  CROSSING 

TIME  INTO  FLIGHT 

CAUSE  OF 

CROSSING 

H 

M 

S 

H  M  S 

EARTH  SEPT  5 

7 

35 

14.000 

240  0  0.326 

SPECIFIED  TIME  OF  FLIGHT  REACHED 

***  FILE  GENERATION  SUMMARY  *** 

AN  ORB-1  FILE  HAS  BEEN  GENERATED 

START  TIME  OF  THE  FILE  -  940826073514.000 

END  TIME  OF  THE  FILE  «  940905073514.000 


***  NUMBER  OF  GTDS  ERRORS  ENCOUNTERED-  0** 
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TRUTH  DATA  GENERATION 


COWELL  DC  &  EPHEM  DECK  (FIT  TO  SGP4 -BASED  EPHEMERIS)  TO  GENERATE  TRUTH 


CONTROL 

DC 

EPOCH 

940826.0 

ELEMENT1 

1 

2 

1 

8183.787 

ELEMENT2 

251.2 

OBS INPUT 

15 

940826073514.0 

ORBTYPE 

2 

1 

1 

60.0 

DMO  PT 

OBSDEV 

21 

22 

23 

500.0 

OBSDEV 

24 

25 

26 

50.0 

END 

DCOPT 

PRINTOUT 

1 

4 

CONVERG 

30 

1 

l.D-3 

END 

OGOPT 

DRAG 

1 

1.0 

ATMOSDEN 

1 

S CP ARAM 

l.D-6 

SOLRAD 

1 

1.0 

POT FI ELD 

1 

2 

MAXDEGEQ 

1 

21.0 

MAXORDEQ 

1 

21.0 

STATEPAR 

1 

STATETAB 

1 

2 

3 

4.0 

DRAG PAR 

1 

SOLRDPAR 

1 

END 

FIN 

CONTROL 

EPHEM 

OUTPUT 

1 

2 

1 

940905.0 

ORBTYPE 

2 

1 

1 

60.0 

OGOPT 

DRAG 

1 

1.0 

DRAGPAR 

0 

ATMOSDEN 

1 

SCPARAM 

l.D-6 

SOLRAD 

1 

1.0 

SOLRDPAR 

0 

POTFIELD 

1 

2 

MAXDEGEQ 

1 

21.0 

MAXORDEQ 

1 

21.0 

OUTOPT 

21 

940826073514.0 

END 

FIN 

VANGARD2 

073513.6735 

0.153069  32.7 

11.1  352.022 

940831073514.0 


500.0  500.0 

50.0  50.0 


100. DO 


5.0  6.0 


OUTPUT  VANGARD2 

073514.0  86400.0 


100. DO 


940905073514.0  450.0 


DATA: 


59001A 


59001A 
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TRUTH  DATA  GENERATION 


COWELL  EPHEMERIS  OUTPUT  SUMMARY: 


ITERATION  NUMBER  6 

********************************************* 


SATELLITE 

EPOCH 

COORD.  SYSTEM 
CENTRAL  BODY 


* 

POSITION  ERROR  RMS  (M) 

960.56962 

* 

VANGARD2  59001 

* 

CURRENT  WEIGHTED  RMS 

1.1041600 

* 

START  TIME  940826 

74244.000 

* 

PREDICTED  WEIGHTED  RMS 

1.1066222 

* 

940826  73513.674 

* 

PREVIOUS  WEIGHTED  RMS 

1.1072310 

* 

END 

TIME 

940831 

73514.000 

* 

SMALLEST  WEIGHTED  RMS 

1.1072310 

* 

MEAN  OF  1950.0 

* 

RELATIVE  CHANGE  IN  RMS 

2 . 7735445 6E- 03 

* 

NO. 

OBS. 

AVAILABLE 

5760 

* 

PENALTY 

0.00000000 

* 

EARTH 

* 

DC  HAS  CONVERGED 

* 

+ 

NO. 

OBS. 

INCLUDED 

5760 

*********************************************** 

NO. 

OBS. 

ACCEPTED  ! 

5757  99  PCT 

OBSERVATION  SUMMARY  BY  TYPE 


TYPE 

X 

Y 

Z 

XDOT 

YDOT 

ZDOT 

TOTAL  NO. 

960 

960 

960 

960 

960 

960 

NO.  ACCEPTED 

960  (100%) 

960  (100%) 

960  (100%) 

960  (100%) 

960  (100%) 

957  (  99%) 

WEIGHTED  RMS 

0.5797 

1.127 

1.444 

0.3973 

0.7494 

1.706 

MEAN  RESDUAL 

-119.6 

472.0 

327.2 

6.770 

9.954 

54.11 

STANDARD  DEV 

264.1 

307.7 

643.5 

18.68 

36.12 

65.94 

OBSERVATION  SUMMARY 

BY  20  DEGREE 

TRUE  ANOMALY 

INTERVAL 

OINTERVAL 

00-20 

20-40 

40-60 

60-80 

80-100 

100-120 

120-140 

140-160 

160-180 

TOTAL  NO. 

240 

234 

264 

282 

312 

348 

378 

408 

432 

NO.  ACCEPTED 

240  (100%) 

234  (100%) 

264  (100%) 

282  (100%) 

312  (100%) 

348  (100%) 

378  (100%) 

408  (100%) 

432 

(100%) 

OINTERVAL  180-200  200-220  220-240  240-260  260-280  280-300  300-320  320-340  340-360 

TOTAL  NO.  426  408  384  342  300  276  252  246  228 

NO.  ACCEPTED  426  (100%)  408  (100%)  384  (100%)  339  (  99%)  300  (100%)  276  (100%)  252  (100%)  246  (100%)  228 

(100%) 


THERE  ARE  8  SOLVE  PARAMETERS 


NUMBER  OF  DYNAMIC  SOLVE  PARAMETERS 
X  Y  2 


vx 


SOLRAD 


NUMBER  OF  STATION  LOCATIONS  BEING  SOLVED  FOR 
NUMBER  OF  BIASES  BEING  SOLVED  FOR  «  0 


***  FILE  GENERATION  SUMMARY  *** 
NO  FILES  HAVE  BEEN  GENERATED 

*  NUMBER  OF  GTDS  ERRORS  ENCOUNTERED  TO  DATE  «  6  * 
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NEEC  OPTIMIZED  INPUT  DECK 


CONTROL 

DC 

VANGARD2 

59001A 

EPOCH 

940826.0 

073513.6735 

ELEMENT1 

3 

6 

1 

8177.913 

0.1522640 

32.8834 

ELEMENTS 

251.8592 

10.8368 

352.1515 

OBS INPUT 

15 

940826073514.0 

940829073514.0 

ORBTYPE 

5 

1 

1 

43200.0 

1.0 

DMOPT 

OBSDEV 

21 

22 

23 

100. 

100. 

100. 

OBSDEV 

24 

25 

26 

10. 

10. 

10. 

END 

DCOPT 

PRINTOUT 

1 

4 

CONVERG 

30 

1 

l.D-4 

END 

OGOPT 

DRAG 

1 

1.0 

ATMOS DEN 

1 

SPDRAG 

0 

SCPARAM 

l.D-6 

100. DO 

SOLRAD 

1 

1.0 

SPSRP 

0 

SPGRVFRC 

1 

1 

3 

3.0 

1.0 

3.0 

SPZONALS 

6 

5 

13 

SPMDAILY 

6 

€ 

4 

POTFIELD 

1 

2 

MAXDEGEQ 

1 

21. 

MAXORDEQ 

1 

21. 

RESONPRD 

86400.0 

STATEPAR 

3 

STATETAB 

1 

2 

3 

4.0 

5.0 

6.0 

DRAGPAR 

1 

DRAGPAR2 

1 

1 

SOLRDPAR 

1 

SSTESTFL 

1 

2 

0 

0.0 

SSTAPGFL 

1 

0 

0 

1.0 

0.0 

1.0 

END 

FIN 

CONTROL 

EPHEM 

OUTPUT 

VANGARD2 

59001A 

OUTPUT 

1 

2 

1 

940905. 

073514.0 

86400.0 

ORBTYPE 

5 

1 

1 

43200.0 

1.0 

OGOPT 

DRAGPAR 

0 

DRAG 

1 

1.0 

' 

ATMOS DEN 

1 

SPDRAG 

0 

SCPARAM 

l.D-6 

100. DO 

SOLRAD 

1 

1.0 

SPSRP 

0 

SOLRDPAR 

0 

SPGRVFRC 

1 

1 

3 

3.0 

1.0 

3.0 

SPZONALS 

6 

5 

13 

SPMDAILY 

6 

6 

4 

POTFIELD 

1 

2 

MAXDEGEQ 

1 

21. 

MAXORDEQ 

1 

21. 

RESONPRD 

86400.0 

OUTOPT 

21 

940826073514.0 

940905073514.0 

450.0 

END 

FIN 
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SUMMARY  OP  NEEC  OPTIMIZED  DSST  FIT  TO  COWELL  TRUTH; 


ITERATION  NUMBER  6 


SATELLITE 

EPOCH 

COORD.  SYSTEM 
CENTRAL  BODY 


VANGARD2  59001 
940826  73513.674 

MEAN  OF  1950.0 
EARTH 


+**********************t*************+******** 


POSITION  ERROR  RMS  (M) 
CURRENT  WEIGHTED  RMS 
PREDICTED  WEIGHTED  RMS 
PREVIOUS  WEIGHTED  RMS 
SMALLEST  WEIGHTED  RMS 
RELATIVE  CHANGE  IN  RMS 
PENALTY 

DC  HAS  CONVERGED 


69. 419240 
0.37449436 
0.37647487 
0.37650718 
0.37650718 
5 . 34604546E-03 
0.00000000 


START  TIME  940826 
END  TIME  940829 
NO.  OBS.  AVAILABLE 
NO.  OBS.  INCLUDED 


74244.000 

73514.000 

3456 

3456 


**********************************************’ 


NO.  OBS.  ACCEPTED  3443  99  PCT 


OBSERVATION  SUMMARY  BY  TYPE 


TYPE 

X 

Y 

Z 

XDOT 

YDOT 

ZDOT 

TOTAL  NO. 

576 

576 

576 

576 

576 

576 

NO.  ACCEPTED 

568  ( 

98%) 

572  {  99*) 

576  (100*) 

576  (100*) 

575  (  99*) 

576  (100*) 

WEIGHTED  RMS 

0.4409 

0.3934 

0.3666 

0.3263 

0.3209 

0.3134 

MEAN  RESDUAL- 

0.8460 

-1.781 

2.432  . 

-4 . 2900E-02 

3 . 5866E-02  * 

-0.2093 

STANDARD  DEV 

44.08 

39.30 

36.58 

3.263 

3.209 

3.127 

OBSERVATION  SUMMARY  BY  20  DEGREE 

TRUE  ANOMALY 

'  INTERVAL 

0 INTERVAL 

00-20 

20-40 

40-60 

60-80 

80-100 

100-120 

120-140 

140-160 

160-180 

TOTAL  NO. 

150 

144 

156 

168 

198 

204 

222 

246 

252 

NO.  ACCEPTED 

149  { 

99*) 

144  (100*) 

156  (100*) 

168  (100*) 

198  (100*) 

204  (100*) 

221  (  99*) 

245  (  99*) 

251  ( 

99*) 

0 INTERVAL 

180-200 

200-220 

220-240 

240-260 

260-280 

280-300 

300-320 

320-340 

340-360 

TOTAL  NO. 

252 

240 

234 

204 

180 

168 

156 

150 

132 

NO.  ACCEPTED 

248  ( 

98%) 

239  (  99*) 

231  (  98*) 

203  (  99*) 

180  (100*) 

168  (100*) 

156  (100*) 

150  (100*) 

132 

(100%) 


THERE  ARE  8  SOLVE  PARAMETERS 

NUMBER  OF  DYNAMIC  SOLVE  PARAMETERS  -  8 

A  H  K  P  Q  LAM  CSUBD  SOLRAD 

NUMBER  OF  STATION  LOCATIONS  BEING  SOLVED  FOR  -  0 

NUMBER  OF  BIASES  BEING  SOLVED  FOR  -  0 

***  FILE  GENERATION  SUMMARY  *** 

NO  FILES  HAVE  BEEN  GENERATED 

*  NUMBER  OF  GTDS  ERRORS  ENCOUNTERED  TO  DATE  -  0  * 
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COMPARE  DECK  FOR  NEEC: 


CONTROL 

COMPARE 

VANGARD2 

59001A 

COMPOPT 

CMPEPHEM 

END 

FIN 

1102102  940829073514.0 

940901073514,0 

7.5 

NEEC  COMPARE  SUMMARY: 


EPHEMERIS  COMPARISON  SUMMARY  REPORT 
Number  of  Points  Compared  577 


MINIMUM  POSITION  DIFFERENCE 


MAXIMUM  POSITION  DIFFERENCE 


YYMMDD 

HHMMSS.SSS 

(km) 

YYMMDD 

HHMMSS.SSS 

(km) 

RADIAL 

940829 

155014.000 

4.082393E-06 

940901 

43514.000 

1.103580E-01 

CROSS  TRACK 

940829 

84244.000 

6 . 979701E-05 

940831 

130514.000 

2 . 553895E-01 

ALONG  TRACK 

940829 

95014.000 

9 . 159481E-06 

940831 

194244.000 

4 . 969624E-01 

TOTAL 

940831 

164244.000 

1 . 110800E-02 

940831 

194244.000 

5 . 081243E-01 

MINIMUM  VELOCITY  DIFFERENCE 

MAXIMUM  VELOCITY 

DIFFERENCE 

YYMMDD 

HHMMSS.SSS 

(km/sec) 

YYMMDD 

HHMMSS.SSS 

(km/sec) 

RADIAL 

940829 

115014.000 

3.964064E-08 

940901 

42744.000 

3 . 010421E-04 

CROSS  TRACK 

940830 

42014.000 

1.722035E-07 

940901 

25014.000 

2 . 334248E-04 

ALONG  TRACK 

940830 

145014.000 

3 . 795B44E-07 

940831 

190514.000 

1.350620E-04 

TOTAL 

940829 

150514.000 

3 . 683908E-06 

940831 

195744.000 

3 . 353047E-04 

POSITION  RMS 

VELOCITY  RMS 

(km) 

(km/sec) 

RADIAL 

5 . 1535E-02 

1 

.  2705E-04 

CROSS  TRACK 
ALONG  TRACK 
TOTAL 


1.0831E-01 
1 . 8288E-01 
2 . 1871E-01 
GTDS  COMPARE 


9 . 0419E-05 
5 . 1923E-05 
1.6436E-04 


NORMAL  COMPLETION  OF  JOB 


16 
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DISTRIBUTION  LIST 


AUL/LSE 

Bldg  1405  -  600  Chennault  Circle 
Maxwell  AFB,  AL  361 12-6424 

DTIC/OCC 
Cameron  Station 
Alexandria,  VA  22304-6145 

AFSAA/SAI 

1580  Air  Force  Pentagon 

Washington,  DC  20330-1580 

PL/SUL 

Kirtland  AFB,  NM  871 17-5776 
PL/HO 

Kirtland  AFB,  NM  87117-5776 

Official  Record  Copy 

PL/VTPT/Brian  Whitney 

Dr.  R.  V.  Wick 
PL/VT 

Kirtland,  AFB,  NM  87117-5776 

Professor  Beny  Neta 
Code  MA/Nd 
Naval  Postgraduate  School 
Monterey,  CA  93943 

Dr.  Paul  Cefola 

The  Charles  Stark  Draper  Laboratory 
555  Technology  Square 
Cambridge,  MA  02139 

Dr.  Ron  Proulx 

The  Charles  Stark  Draper  Laboratory 
555  Technology  Square 
Cambridge,  MA  02139 


1  cy 


2  cys 


1  cy 


2  cys 


1  cy 


2  cys 
1  cy 


1  cy 

2  cy 


1  cy 
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Mr.  Rick  Metzinger 
The  Charles  Stark  Draper  Laboratory 
555  Technology  Square 
Cambridge,  MA  02139 

Mr.  Wayne  McClain 
The  Charles  Stark  Draper  Laboratory 
555  Technology  Square 
Cambridge,  MA  02139 

Mr.  Leo  Early 
Microcosm  Inc.,  Suite  230 
Rolling  Hills  Office  Plaza 
2601  Airport  Drive 
Torrance,  CA  90505 

Dr.  Shannon  Coffey 
Code  8233 

Naval  Research  Laboratory 
Washington,  DC  20375-5355 

Mr.  Steve  Casali 

Kaman  Sciences  Corporation 

P.  O.  Box  7463 

Colorado  Springs,  CO  80933-7463 

Dr.  Chia-Chun  (George)  Chao 
P.  O.  Box  92957-M4/946 
Los  Angeles,  CA  90009-2957 

Dr.  Stephen  H.  Knowles 
U.S.  Naval  Space  Command 
Code  63T 

Dahlgren,  VA  22448-5170 


Dr.  Paul  Schumacher 
U.S.  Naval  Space  Command 
Code  63T 

Dahlgren,  V A  22448-5170 


1  cy 
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2cy 


1  cy 


1  cy 


1  cy 


James  R.  Wright 
305  Exton  Commons 
Exton,  PA  19341 

Dr.  Dave  Carter 

The  Charles  Stark  Draper  Laboratory 
555  Technology  Square 
Cambridge,  MA  02139 

W.  Spencer  Campbell 
The  Aerospace  Corp. 

2350  E.  El  Segundo 
El  Segundo,  CA  90245-4691 

Dr  Ahmed  H.  Salama 
Jet  Propulsion  Lab 
4800  Oak  Grove  Drive 
Pasadena,  CA  91109-8099 


